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PREFACE 


The last major Transonic Symposium was held at NASA Ames Research Center 
in February 1981. Since then, significant advances have been made in computer 
hardware, theoretical and computational methods, applications, experimental 
facilities, and testing techniques. Although much research remains to be 
done, these advances now provide us with capabilities in the transonic regime 
which we hardly envisioned at that time. In order to assess the state of the 
art in transonic flow disciplines and to glimpse at future directions, 
the NASA Langley Research Center held a Transonic Symposium on April 19-21, 
1988. Emphasis was placed on steady, three-dimensional external, transonic 
flow and its simulation, both numerically and experimentally. 

Papers were presented by researchers from NASA, industry, and 
universities. The symposium included technical sessions on wind tunnel and 
flight experiments; computational fluid dynamics applications (industry over- 
views and configuration analysis design); inviscid methods and grid generation 
viscous methods and boundary-layer stability; and wind tunnel techniques 
and wall interference. 

The proceedings are published in two volumes as follows because of the 
range of classifications: 

Volume I, Unclassified ( Part^land21^ 

Volume II, Unci assified,J|^0^g^0H^ 

A list of attendees is included in this document. 

General Chairman: P. J. Bobbitt 

Organizing Committee: M. D. Salas 

L. E. Putnam 
J. T. Foughner 


Technical Committee Chairmen: 

Theory and Computational J. L. Thomas 

Applications E. G. Waggoner 

Experiments L. W. McKinney 
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ABSTRACT 

The National Transonic Facility (NTF) has been operational in a combined 
checkout and test mode for about 3 years. During this time there have been 
many challenges associated with movement of mechanical components, operation 
of instrumentation systems, and drying of insulation in the cryogenic environ- 
ment. Most of these challenges have been met to date along with completion of 
a basic flow calibration and aerodynamic tests of a number of configurations. 
This paper reviews some of the major challenges resulting from the cryogenic 
environment with regard to hardware systems and data quality. Reynolds number, 
effects on several configurations are also discussed. 



INTRODUCTION 


The National Transonic Facility (NTF), which was constructed by NASA with 
a goal of meeting the national needs for High Reynolds Number Testing, has been 
operational in a checkout and test mode for about 3 years. The order of mag- 
nitude increase in Reynolds number over existing transonic wind tunnels pro- 
vided by the NTF, figure 1, is the result of operating at cryogenic tempera- 
ture and stagnation pressures to 8.8 atmospheres. Although the cryogenic 
temperatures provide some significant and well documented benefits from a 
Reynolds number standpoint, the harsh environment also provides equally sig- 
nificant challenges for reliable operation of large mechanical systems and 
instrumentation. 

The approach followed during the 3 years of operating in a combined 
checkout and test mode had some obvious advantages for a facility like the 
NTF where there is not a significant experience base. Known problems can be 
solved while identifying and solving those problems that will only show up by 
using the tunnel in a testing mode. The end result is a fully operational 
facility at an earlier date. However, there are also some disadvantages. 

Most significant among them is that testing during this time period is at a 
much reduced level of efficiency. For a facility like NTF where a very high 
level of efficiency is important, it is difficult for both operators and ob- 
servers to maintain perspective during this test period. 

Over the past 3 years a host of operating problems resulting from the 
cryogenic environment have been identified and solved. These ranged from 
making mechanical /el ectrical systems functional to eliminating temperature- 
induced vibration, to minimizing the effect of moisture outgassing from the 
thermal insulation. Additionally, a preliminary flow calibration has been 
completed, and a series of aerodynamic tests has demonstrated data quality 
and provided Reynolds number effects on several configurations. Currently, 
a major effort is under way, through the summer of 1988, which is devoted to 
improving operating efficiency with a goal of being ready to efficiently 
support both research and development testing requirements by the fall. 

This paper will review some of the more significant efforts during this 
time period and summarize the NTF status concerning hardware and instrumenta- 
tion systems, operating constraints imposed by the cryogenic environment, data 
quality, and some Reynolds number data. 
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TESTING AND CHECKOUT EXPERIENCE 

The testing and checkout experience is summarized in figure 2. The ini- 
tial start-up experience (prior to 1985) is reported in references 1 to 3. As 
stated previously, the testing that has been accomplished to date has a two- 
fold purpose of providing aerodynamic data and exposing testing problems 
associated with the tunnel and instrumentation systems. Nine of the configu- 
rations tested are shown in figure 3 which collectively utilized the maximum 
capability of the NTF at both ambient and cryogenic conditions over the Mach 
number range. Two other aircraft configurations have also been tested but 
will not be discussed here. The primary model used for checkout of both tun- 
nel and Instrumentation systems was the Pathfinder I which has a high aspect 
ratio wing with a supercritical airfoil (reference 4). This model was first 
installed in the tunnel during the first quarter of 1985. 

The major areas receiving attention during the checkout are listed at the 
bottom of figure 2 and include model access, process controls, moisture in the 
tunnel, model vibration, and tunnel /test instrumentation. All of these prob- 
lems were worked simultaneously as indicated by the figure and were phased in 
with the testing schedule as appropriate. The tunnel was unavailable for 
cryogenic operation during most of the first half of 1987 due to a failure of 
an expansion joint in the liquid nitrogen supply system. 

Model Access. - Access to the model requires the movement of large compo- 
nents - wTtfrTrTTfie - tunnel , (figures 4, 5, and 6), over the temperature range 
from ambient to fully cryogenic. The details of this system are defined in 
references 3 and 5. The test section plenum is isolated from the rest of the 
tunnel circuit by large isolation valves. The “process of putting these valves 
in place involves unlocking and translating a large section of the contraction 
cone and the high speed diffuser away from the plenum; this process uses dual 
electrical driven actuators with a seven-foot stroke which must operate in 
phase on each component. Additionally, the ability to make up limit switches 
or components to fairly close tolerance where the components may be exposed to 
large temperature excursions is required. With the plenum vented to atmo- 
spheric pressure, the 9- by 12-foot doors in the pressure shell are opened and 
the test section sidewalls are dropped so that access housings may be inserted 
to encapsulate the model as shown in figures 5 and 6. The’ reliable movement 
of these large components at cryogenic temperatures has required several modi- 
fications to the basic actuation concept. These modifications have resulted 
from operational experience and were implemented over the past 3 years. A 
final series of modifications to these components is being made during the 
current enhancement period which should make it a fully operational system. 

Process Controls . - The primary controls for the tunnel are closed loop 
and provide fast response interactive control of pressure, temperature, and 
Mach number and control of model attitude (see figure 7). The controls for 
the test section variables (tunnel walls and re-entry flaps) are also closed 
loop but have slower response requirements. ^ 


A detailed description of the process controls is provided in reference 6 
and only a brief summary will be presented here. The approach used in design 
of the pressure, temperature, and Mach number controls involved using a mathe- 
matical model of the process to determine circuit response characteristics in 
order to establish the design criteria for the control hardware and the ini- 
tial control laws. This approach as indicated in figure 8 required measure- 
ment of actual tunnel response characteristics for verification of the control 
laws and update of the math model. This effort was a first priority in the 
early tunnel operation and indicated that the accuracy required could not be 
obtained due to insufficient system resolution and excessive instrumentation 
noise. To correct this problem the microprocessors, servo control valves, and 
instrument! on were upgraded. The upgraded stagnation pressure and Mach number 
control systems were completed in the first half of last year. The Mach num- 
ber system provides a rapid response and will control around set point to 
within ±0.001. The last system to become fully operational was the tempera- 
ture control. Initially, this system was based on a measurement of the liquid 
nitrogen flow rate entering the tunnel. Minimal success was obtained in 
making this measurement and an alternative approach which calculates the flow 
rate has been recently implemented and provides adequate performance. 

Moisture Contamination . - During the early operation of the NTF at Cryo- 
genic temperatures, a coating was observed on several models that had a frost- 
like appearance. Quantifiable measurements of the extent of the coating were 
not obtainable. However, it was clear that frost-like crystals were forming 
and that reflected light was making it visible with a television system. An 
extensive study was undertaken to determine the contaminating substance and 
its source as well as its possible effect on the aerodynamic data. This study 
identified the contaminant to be water. 

Further evaluation of the incoming liquid nitrogen, various tunnel purge 
technique, and samples of the thermal insulation system identified the insula- 
tion as the water source. The closed-cell polyisocyanhrate foam used for the 
thermal insulation has been shown to have approximately 1.5 percent water 
by weight. (As a reference, dry wood has approximately 6-percent water by 
weight.) The problem, as shown in figure 9, is that at the cryogenic tempera- 
tures only a very small quantity of water is required to saturate the free- 
stream gas. For example, at -70°F only about 0.2 pounds of water is required 
to saturate all of the gas in the tunnel at atmospheric pressure, while at 
+70°F about 200 lbs of water is required. Two questions naturally arise at 
this point. How to eliminate the moisture, and what is its effect on the 
aerodynamic data measurements? Extensive investigations have been conducted 
in both areas. With regard to the question of eliminating the moisture, the 
studies have indicated that the simplest solution appears to be drying the 
tunnel and then maintaining a dry environment. Drying has been accomplished 
in the past by maintaining the tunnel at the maximum warm temperature (approx- 
imately 150°F) for several hours in dry nitrogen gas with a periodic purge. 

The diffusion of water from the insulation in the NTF is a maximum at this 
condition as indicated in figure 10. The figure also illustrates the reduc- 
tion in diffusion with decreasing temperature which is, of course, a favorable 
effect. Since in the nitrogen operating mode the tunnel is continuously 
vented, i.e., purged to maintain constant stagnation pressure, test conditions 
with very small amounts of condensation can be obtained. Aerodynamic studies 
conducted under these conditions, which will be discussed later, have indi- 
cated that the resulting effects on the aerodynamic data are not measurable. 


5 


Therefore, the plannned approach when the tunnel is brought back into opera- 
tion at the conclusion of the current enhancement period will be to dry it for 
several days and keep it closed to atmospheric air. It should be pointed out 
that the drying process takes place in static as well as operating conditions; 
therefore, it continues to dry over weekends and nonworking shifts. 

The two cryogenic tunnels in Europe with internal insulation, KKK at Koln 
and T2 at Toulouse, have experienced similar moisture concerns and have also 
reached the conclusion that they can achieve conditions of dryness where 
effects on the aerodynamic data are not measurable. We are convinced from 
experience to date, with the NTF and other studies, that procedures can be 
followed that will eliminate moisture contamination as a concern regarding 
aerodynamic data quality. 

Model Vibration. - Significant vibration of the model/balance combination 
In the lateral plane has been encountered at some conditions since the initial 
operation of the NTF. After some period of operation the vibration was found 
to be more severe when the structure was cold. The model pitch system is 
illustrated in figure 11 and is composed of an arc sector driven by a hydrau- 
lic cylinder. Restraint is provided by a series of bearing pads located at 
both the top and bottom of the sector. The loads are transmitted through the 
pads to the internal tunnel structure. There is also a fairing on the down- 
stream part of the sector that is fixed to the tunnel structure and provides a 
cavity for instrumentation leads. The attachment of the fixed fairing to the 
arc sector is a slip joint which allows the sector to move independently of 
the fairing. The bearing pads, while providing restraint, also have clearance 
to allow for thermally induced movement of the internal structure. 

The vibration problem has been investigated both experimentally and ana- 
lytically. In the experimental investigation, the Pathfinder model was used 
as a test vehicle. Both it and the model support system were extensively 
Instrumented as follows: 

a. Six component force balance and 3-axis accelerometer package in the 
Pathfinder model 

b. Pressure tranducers in the fixed fairing and test section walls 

c. Accelerometers on the fixed fairing, bearing pads, and the surrounding 
tunnel support structure 

d. Strain gages and thermocouples on the tunnel support structure 

The analytical Investigation involved detailed calculations of the 
dynamic structural respose and of the unsteady aerodynamic characteristics of 
the model support system. 

Experimental observations were made of the dynamic structural response 
characteristics, Mach number, dynamic pressure, test temperature, and coldsoak 
time. The installation of vortex generators and splitter plate on the fixed 
fairing afterbody helped to identify unsteady flow at the rear of the arc sec- 
tor as one of the sources of dynamic excitation. Eventually, the primary fac- 
tor governing the dynamic response of the model was found to be the clearance 
tolerance between the bearing pads and the surface of the arc sector. A 
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procedure was evolved for setting the clearance adjustment which reduced the 
model dynamic response to acceptably low levels independent of temperature 
cycling. Some of the test results are illustrated in figure 12 which shows 
the dynamic yawing moment as a function of stagnation temperature for several 
of the test configurations of the arc sector. Although the low level shown in 
the figure for the last case is completely satisfactory, it may be sensitive 
to adjustment with time due to temperature cycling. 

Instrumentation. - For several years prior to initial operation of the 
NTF an extensive research and development program was undertaken by the 
Instrument Research Division at Langley to develop instrumentation systems for 
basic measurements of forces and moments, pressures, and angle of attack that 
would function with both reliability and accuracy in the cryogenic environ- 
ment. The results of this program were the development of strain gage bal- 
ances that were not temperature controlled and heated instrument packages for 
pressure and angle of attack measurements. These instrument systems were 
developed in cryogenic chambers and verified to the extent possible, in the 
0.3-Meter Transonic Cryogenic Tunnel and indicated good performance and sound- 
ness of the basic concepts. Upon application of these instruments to models 
in the NTF, there were several system problems that had to be resolved. The 
most significant of these were the effect of leads required to support the 
pressure instrumentation system on the balance axial force component, and the 
calibration system for the electronic scanning pressure (ESP) system. Recent 
tests, both in the NTF and in the cryogenic checkout chamber, have indicated 
satisfactory resoul ti on of both of these problems. 

As will be discussed in a later section, a strong concern from the outset 
has been data quality. In this regard, an extensive effort has been carried 
out to eliminate data scatter due to electrical noise and extraneous signals. 
Care was taken from the beginning to provide a "clean" instrumentation ground, 
but as is often the case, extraneous signals can creep into the system. A 
significant effort has been made over the past year to track down and elimi- 
nate noise sources that were infiltrating the data system. The results of 
this effort are indicated in figure 13. The figure illustrates a reduction in 
data scatter by a factor of 5 to a level less than 10 microvolts. This level 
corresponds to a balance error of less than 0.1 percent of full scale. 


OPERATING CONSTRAINTS 

The NTF has the general appearance of, and Is often talked about, as a 
typical continuous operating fan driven wind tunnel that has the potential to 
mass produce data on a continuous basis. Although In principle the potential 
exists, it is somewhat misleading to think about a large cryogenic tunnel in 
this context. In reality the operation must be viewed much as a blow- down 
tunnel with efforts directed toward minimizing run time and maximizing data 
gathering rates. The two main constraints with regard to run time are liquid 
nitrogen logistics and costs, and model access. It was shown in the early 
1970's that the application of cryogenics to wind tunnel testing provided the 
most cost effective approach to high Reynolds number testing. This does not 
mean, however, that testing at high Reynolds numbers will carry the same cost 
and degree of complexity as testing in the many ambient temperature and pres- 
sure tunnels around the country. This is illustrated in figure 14 where the 
energy cost ratio (liquid nitrogen is a form of energy) is shown as a function 
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of Reynolds number ratio for the NTF compared to the Langley 16-Foot Transonic 
Tunnel. An order of magnitude increase in Reynolds number has an energy cost 
ratio of approximatley 44. In the cases where high Reynolds number is re- 
quired for research or development testing, the cost is cheap compared to 
flight test. But because of this cost difference, test programs and objec- 
tives must be carefully defined and supported by adequate precursor studies at 
low Reynolds numbers. Therefore, it is not well suited for researchers to 
exercise broad freedom or for indiscriminate development programs. 

LN 2 Supply System. - Liquid nitrogen is supplied to the NTF by a commer- 
cial "air "s'eparati on pi ant located adjacent to the Langley property and con- 
nected to the NTF site by a pipeline. The capability of this system is shown 
in figure 15. The challenge is to optimize the interface of the plant which 
operates continuously 24 hours per day at a 300 ton per day rate with the 
Intermittent operation of the NTF, which can use LN 2 to a maximum rate of 
30 tons per minute. This requires, of course, storage tanks to serve as a 
buffer. The current system has a 2100 ton storage tank at the plant and a 
700 ton tank located at the NTF with capability to transfer approximately one 
tank per 24 hour day. The system can sustain a use rate of 2100 tons per week 
and if the tunnel has not been using nitrogen for a while, can build up to a 
maximum quantity available of 4800 tons for a week. This tends to optimize on 
two week test programs using liquid nitrogen assuming all tanks are full at 
the start. The maximum transfer rate of the pipeline is 656 tons or approxi- 
mately one NTF tankful per day; therefore, the maximum use rate is one tank- 
ful per day. As shown in figure 16, if operation is started on Monday morn- 
ing with the equivalent of four NTF tankfuls available and used at a rate of 
1 tank per day, with 0.5 tanks being resupplied by the plant, and no use over 
the weekend, the total supply is exhausted by the end of the second week. The 
options at that point are either to operate on one-half tank per day (the 
plant output) which is not practical, or to allow 8 days for the nitrogen sup- 
ply to be replenished and use the tunnel in the air mode or some other capac- 
ity during this period. This use scenario can be changed, of course, by 
Increasing the plant capacity and/or storage and transfer rate. A decision to 
increase the plant capacity must take into consideration long term use rates. 
The most immediate benefit will be provided by an Increase In storage capacity 
and associated transfer rates. Current plans are to triple the NTF onsite 
storage capacity in FY 1990. This will tend to optimize at about a three week 
test program in nitrogen. 

The most important aspect of the tunnel operation from cost considera- 
tions is the speed of changing test conditions and data acquisition. This is 
driven primarily by the degree and quality of automation. Figure 17 illus- 
trates the impact of time per data point for a typical test condition of 
M = 0.8 and temperature = -250°F. This illustration includes time to change 
angle of attack and Mach number and adjust temperature and pressure. When the 
NTF was first brought on-line and manual control was used, times averaged 
greater than 45 seconds per point. After the current enhancement period it is 
expected to be under 30 seconds per point with a goal of 7 to 10 seconds per 
point with further refinements in the control and operating system. The re- 
duction in electrical signal noise discussed earlier has a direct impact on 
this in that it reduces the number of data samples required to be averaged to 
obtain high quality data. 
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DATA QUALITY 


In this section the status of flow calibrations and efforts to assess 
data quality will be discussed. The discussion will include both steady and 
dynamic aspects of the flow calibration, the ability to measure Mach numbers 
and angle of attack, importance of wall boundary effects, and the effect of 
moisture contamination or "frost" on the data. 

Flow Calibration. - As reported previously, references 3 to 10, the MTF 
has been operated throughout the operating envelope as shown in the upper left 
of figure 18. The initial calibration looked at the distributions of tempera- 
ture, pressure, Mach number, and flow angularity on the tunnel centerline as 
measured by the model upright and inverted. More recent efforts concerned the 
details of the Mach number calibration and the effect of temperature on flow 
angle. 

The ability to vary the test section wall angle provides the capability 
to maintain zero Mach number gradient through the test section, thus eliminat- 
ing model buoyancy effects in the empty tunnel for all test conditions. A 
typical variation of wall angle with Mach number is shown in the left of fig- 
ure 19. This particular case is for a warm temperature of 100°F and varies 
from approximately 0.2° convergence at M = 0.2 to 0.4° divergence at M * 1.15. 
A typical correction to Mach number as calculated from the plenum static or 
reference pressure is also shown. In general the quality of the steady flow 
appears to be excellent and sufficient adjustments are available on the test 
section geometry, i.e., wall angle and re-entry flap settings, to eliminate 
the existence of any Mach number gradients in the empty tunnel. 

In order to investigate the effect of tunnel structural deformation due 
to temperature gradients in the structure on flow angle with temperature, fre- 
quent measurements of flow angle have been made during the last several test 
programs. The data presented in figure 20 are typical of results obtained 
from these measurements. With the exception of one point the effective flow 
angle is always equal to or less than 0.02° which is approaching the measure- 
ment accuracy. This effect will continue to be monitored until sufficient 
history is developed to have a firm basis for determining the required fre- 
quency of flow angle measurements. 

NTF Dynamic Flow Quality. - In 1980 while the NTF was still under con : 
struction, the plans' for flow quality measurement were described extensively 
in reference 11. The measurements were to consist primarily of fluctuating 
pressure and velocity measurements. Since that time, hot-film probes and 
fluctuating pressure gages have been operated in the test section at the 
locations indicated in figure 21, and further tests are planned as also shown. 


Some results of the measurements made to date are shown in figures 22 and 
23. These results have been excerpted from work by W. B. Igoe on a proposed 
doctoral dissertation to be submitted to the George Washington University. 
Figure 22 shows the root-mean- square fluctuating static pressure measured in 
the NTF as a function of Mach number for a unit Reynolds number of 6 million 
per foot. The fluctuating pressures were measured in air on the test section 
sidewall and have been divided by free- stream static pressure. Fluctuating 
static pressures have been measured in the free stream In a number of other 
large transonic wind tunnels using a 10° cone on the centerline (see 
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reference 12 for example). Some of these results have been included for 
comparison in figure 22 for a Reynolds number range of 1.3 to 5 million per 
foot. 


The results of hot- film probe measurements at low Mach number in the test 
section are shown in figure 23. These measurements were made in the free 
stream in air with the slots closed. Results are shown for the streamwise 
measurements at total pressures of 1 and 8 atm. Although there is consider- 
able scatter at a pressure of 1 atm, the measurement levels are about 0.1 per- 

cent, which was the target level for the NTF. Streamwise hot-wire probe mea- 
surements made in the Langley Low Turbulence Pressure Tunnel (LTPT) at total 
pressures of 1 and 10 atm (reference 13) are shown for comparison. 

Further measurements in the NTF are planned using fluctuating static 

pressure probes in the test section free stream, and hot-film probes in the 
test section, settling chamber, and in the vicinity of the cooling coil and 
screens. The measurements will be made over the full operating range of the 
NTF. By the time the measurements are completed, the dynamic flow quality of 
the NTF will have been fully documented. 

Mach Number and An gle of Attack Sensitivities. - Prior to lookina at 
either Reynolds number effects or the effects of frost on data quality, it was 
desirable to obtain an assessment of the sensitivity of the model to the basic 
test parameters, Mach number and angle of attack, and some qualitative assess- 
ment of our ability to measure them. As stated in the introduction, the 
primary model used for assessing data quality was the Pathfinder I Model. 

This model was built early in the program for the purpose of developing model 
design and fabrication technology as well as providing a research model. It 
has an aspect ratio 10 wing with a supercritical airfoil section typical of 
supercritical airfoil design technology of the mid 1970's. As a result, shock 
movement is very sensitive to small changes around the design point in Mach 
number and angle of attack. This made it a good model for use in assessing 
our ability to measure Mach number and angle of attack in the NTF. Figures 24 
to 26 illustrate the results obtained from this assessment. The figures show 
pressure distributions for an inboard and an outboard wing station at a Mach 
number around 0.82 with transition fixed at 10 percent of chord. The data of 
figure 24 shows that a Mach number increment of 0.0038 results in a shock 
movement of about 6 percent of chord for the outboard wing station. In light 
of the high degree of sensitivity to Mach number, a series of runs were com- 
pared where all variables except Mach number were held constant and Mach 
number varied in very small Increments of 0.0001 to 0.0006. Wing pressure 
distributions for these cases are shown in figure 25. The data show an 
orderly progression of shock movement and the well defined curves suggest that 
both wing pressures and Mach number are being measured with a high degree of 
accuracy. 

The sensitivity of this model to angle of attack is illustrated in fig- 
ure 26. In this case, the angle of attack increment was 0.054° with a Mach 
number difference of 0.0013. Again, the shock movement is about 5 to 5 per- 
cent. However, about one half of the movement can be attributed to Mach 
number effects. These data support the point that the highest quality instru- 
mentation is an absolute requirement for using the NTF to understand incremen- 
tal effects of Reynolds number and compressibility. Further they support the 
conclusion that a high degree of accuracy is currently available in the Mach 
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number, angle of attack, and model pressure measuring systems. The data also 
underscore the importance of being able to accurately assess wall boundary 
effects. Recent research on this subject by W. B. Kemp (ref. 14), P. A. 

Newman and associates (ref. 15), is described by Dr. Newman in a separate 
session of this symposium. Sophisticated techniques have been developed which 
utilize measured tunnel wall static pressures to calculate model induced vari- 
ations of Mach number and upwash through the test section. Figure 27 shows 
typical contours, in the region of the model, of wall induced Mach number 
corrections. For this size model at the conditions illustrated the correc- 
tions are relatively small, aM = 0.001, but as illustrated in the previous 
figures, corrections of this magnitude are significant if high quality data 
are to be obtained. 

Moisture Contamination "Frost" Effects. - Having established confidence 
that small incremental effects of the basic test parameters, Mach number, and 
angle of attack could be both controlled and measured, an investigation to 
assess the possible effects of frost on the data was undertaken. Care was 
also taken to insure that comparable test parameters were obtained where the 
only significant variable was that in one case frost was visible on the model 
and in the other case it was not visible. Wing pressure distributions from 
these two cases are compared in figure 28. Although care was exercised in 
setting the test parameters it should be noted that the Mach number is dif- 
ferent by 0.0004 and the angle of attack by 0.01°. These differences are 
believed, based on previous discussion, to account for the small difference in 
shock location shown in the data of figure 28 for the outboard wing station. 
After accounting for the difference in shock location, there is a small dif- 
ference in pressure level ahead of the shock that may be a small effect of 
frost. In general the two cases are in very close agreement and provide con- 
fidence that when planned tunnel drying procedures are used and the inside of 
the tunnel kept closed to atmospheric air, frost on the models will not be a 
problem with regard to data quality. 


REYNOLDS NUMBER EFFECTS 

During this initial checkout and testing phase several Reynolds number 
sensitivity studies have been conducted using both the conventional air mode 
and the cryogenic capability 

Air Operation. - The operating envelope in air is shown in figure 29. It 
is restricted above a Mach number of 0.4 by the drive power in the variable 
speed motors. A maximum Reynolds number capability of about 20 million per 
foot at a Mach number of 0.38 is available. This Reynolds number is the 
maximum available in this speed range and tunnel size in the United States. 

One of the attractive features of this capability is that a constant dynamic 
pressure line tends to be close to the maximum drive power boundary. There- 
fore models designed for testing at high subsonic or transonic speeds at 1 to 
2 atmospheres can also be tested at high Reynolds numbers at the lower Mach 
numbers at the same dynamic pressure and model loads. One case where this was 
done was the EA-6B wing modification program. The objectives and overall 
results from this program are described in reference 16. In summary, it 
involved modifying the wing leading-edge slat and trailing-edge flap airfoil 
sections with a major goal of improving maximum lift at loiter and maneuver 
conditions. A photograph of the EA-6B configuration mounted in the NTF is 
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shown in figure 30. The effect of Reynolds number on the lift characteristics 
at a Mach number of 0.30 for the basic and modified configuration are shown in 
figure 31. The 1.4 million Reynolds number case corresponds to testing at 1 
atmosphere stagnation pressure. Two points are significant. First, as would 
be expected, Reynolds number effects on CL max are large for all configura- 
tions. Second, modifications to the leading-edge slat showed no benefit at 
the low Reynolds number condition; however, at the higher Reynolds number the 
slat benefit was approximately equal to that of the flap. These data are an 
example of the potential pitfalls of relying on low Reynolds number data for 
configuration refinement. The effect of Reynolds number on lateral and direc- 
tional stability as measured by Cig and Cng are shown in figures 32 and 33. 

The effects are not as dramatic as those shown for C|_ max , but in general a 
stabilizing increment in lateral stability was obtained with increasing Rey- 
nolds number. The exception to this was the basic configuration above a = 16° 
and the modified wing configuration with glove and vertical tail extension 
above a equal about 14°. A more stabilizing effect in Cng was obtained for 
all configurations with increasing Reynolds number. 

Cryogenic Operations . - Results from three of the configurations tested 
during this period will be briefly discussed here. They include the Path- 
finder I, a Lockheed high-wing transport which was a Lockheed wing tested on 
the Pathfinder I fuselage, and a rather large submarine model. 

The Pathfinder I (figure 34) was tested over a range of conditions; how- 
ever, most of the test was aimed at evaluating instrumentation and moisture 
concerns. Most of the high Reynolds number data was obtained in the early 
test program and in retrospect may have been contaminated with moisture 
effects and instrumentation errors. Therefore, the data shown in figure 35 
are for more intermediate Reynolds numbers (RN = 5 and 18 million) which were 
known to be free of instrumentation errors. These data are for a Mach number 
of 0.82, a constant angle of attack and essentially constant lift coefficient 
with transition fixed at approximately 10-percent chord. These conditions are 
the design point for the wing. The natural transition point of the 18 million 
Reynolds number case was estimated to be essentially at the trip location. 

The effects of Reynolds number are relatively small with only a slight aft 
movement of the shock indicated at the outboard wing station. 

The data for the Lockheed configuration shown in figure 36 show a much 
more pronounced effect of Reynolds number. A photograph of the model is 
presented in figure 37. For the case with transition fixed, the strip was 
located at a constant 1-inch aft of the wing leading edge. The difference in 
shock location between transition fixed and free at a Reynolds number of 5 
million is about 18 percent of chord. A significant difference was obtained 
by increasing the Reynolds number to 30 million. A further increase to 40 
million produced a negligible effect. These data clearly support the well- 
known need to be able to monitor transition location as a function of Reynolds 
number if effects at Reynolds numbers less than full scale are to be 
interpreted. 

The submarine test in NTF was unique in that full scale Reynolds number 
was obtained in a wind tunnel for the first time. A photograph of the model 
and Reynolds number velocity envelope is presentd in figure 38. Details of 
the test and data have been omitted for classified security reasons. In the 
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test both static and dynamic measurements of pressures were obtained through 
the boundary layer in the plane of the propeller for evaluation and develop- 
ment of scaling laws. As can be seen on the right of figure 38, data in the 
past using a 6 percent scale model have been obtained at Reynolds numbers 
about an order of magnitude below full scale. The test in the NTF extended 
these data well into the region of full-scale submarine operation. 


SUMMARY 

The National Transonic Facility has been operational in a checkout and 
test mode for the past 3 years. During this time there have been many chal- 
lenges associated with testing in a large cryogenic wind tunnel. For the most 
part they centered around the effect of large temperature excursions on the 
mechanical movement of large components, the reliable performance of instru- 
mentation systems, and an unexpected moisture problem with dry insulation. 

Most of these challenges have been met, and it is expected that the rest of 
them will be met during a major effort that is ongoing through the summer of 
1988 to Improve operational efficiency. Also, during the past 3 years a 
preliminary flow calibration has been completed and a data quality evaluation 
conducted along with high Reynolds number aerodynamic tests of several config- 
urations. Tests were also conducted that provided major inputs to several 
programs that were not discussed for classification reasons. There is still a 
requirement for a major effort to develop and implement flow visualization and 
diagnostic techniques for maximum utilization of the facility. The current 
NASA facility revitalization program includes funding for these activities in 

FY'89 and 90. However, from a basic facility standpoint, we believe that it 

will be ready to efficiently support research and development requirements by 

the fall of 1988. 
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Figure 1. National Transonic Facility perspective and Reynolds number 
operating envelope. 
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Figure 3. Photographs of models tested In the NTF. 
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Figure 5. Schematic of model access system illustrating access tubes in the 
inserted and retracted positions. 
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Figure 8. Chronology of events for control system checkout and upgrade. 
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P = 1 ATM. 



Figure 9. Water vapor required to saturate free-stream gas as a function of 
stagnation temperature. 
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Figure 10. 


Rate of water transfer from insulation to free-stream gas as a 
function of stagnation temperature. 
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Figure 15. Characteristics of the liquid nitrogen supply system. 
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Figure 16. NTF maximum utilization based on LN£ supply system. 
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Figure 17. The effect of data sampling rate on liquid nitrogen costs 
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Figure 18. Typical calibration results. 
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Figure 19. Variations of test section wall angle required for zero Mach 

number gradient at the model location and correction to reference 
Mach number with Mach number. Stagnation temperature = 1CIU°F. 
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Figure 20. Typical variation of tunnel flow angle with temperature as 
measured by model upright and inverted. M = 0.82. 
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Figure 21 . Location of dynamic flow quality measurements. 
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Figure 22. Variation of fluctuating static pressure in the test section with 
Mach number for several major wind tunnels. 
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FREESTREAM TURBULENCE 



M 

Figure 23. Comparison of turbulence measurements in the test section for the 
NTF and LTPT. 



Figure 24. The effect of Mach number on wing pressure distributions for the 
Pathfinder I model. RN = 6.5 x 10 ; transition fixed at 0.10c. 
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Figure 30. EA-6B model mounted in NTF. 
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Figure 31. The effect of flap/slat modifications on the low speed lift 
characteristics for the E.A-6B aircraft model . M = 0.30. 
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Figure 32. The effect of Reynolds number on rolling moment due to sideslip 
for the EA-6B. M = 0.30. 
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Figure 33. The effect of Reynolds number on yawing moment due to sideslip for 
the EA-6B. M = 0.30. 



Figure 34. Photograph of the Pathfinder I model mounted in the NTF. 
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Effect of Reynolds number on wing 
Pathfinder I model. M = 0.82; a = 
0.10c. 


pressur 

2.57®; 





distributions for the 
transition fixed at 



inmiiRiiiiiii ilium mu iiiimiii muni i i i r m m iiiiiiiiini hi d i i mn nnmniR 




LOCKHEED HIGH WING TRANSPORT CONFIGURATION 

I : 



Figure 36. Photograph of Pathfinder I fuselage with the Lockheed transport 
wind installed in the NTF. 
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Figure 37, Effect of Reynolds number on wing pressure distributions for the 
Lockheed high wing transport configuration. M = 0.80; C M = 
Constant. 
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Summary 

Aerodynamic forces and moments for a slender wing-body configuration are summarized from an 
investigation in the Langley National Transonic Facility (NTF). The results include both longitudinal 
and lateral-directional aerodynamic properties as well as sideslip derivatives. Results have been selected 
to emphasize Reynolds number effects at transonic speeds although some lower speed results are also 
presented for context. The data indicate nominal Reynolds number effects on the longitudinal aerody- 
namic coefficients and more pronounced effects for the lateral-directional aerodynamic coefficients. The 
Reynolds number sensitivities for the lateral-directional aerodynamic coefficients were limited to high 

angles of attack. 


Introduction 


Recent interest has developed in advanced aerospace vehicles which are capable of very high speed 
flight. Examples of such vehicles include a variety of advanced transport concepts designed for supersonic 
cruise as well as transatmospheric vehicles such as the proposed X-30. These vehicles all tend to be 
slender due to high speed considerations, although they stiff eihbrace a wide range of configurational 
concepts (i.e., wing-bodies, waveriders, accelerators, etc.). The aerodynamic challenges for such vehicles 
are by no means limited to high speed concerns such as cruise design or aerothermal heating. Most 
aerodynamic subdisciplines (e.g., stability and control, propulsion integration, transonic flow, high angle 
of attack, etc.) present unique and often conflicting challenges for these vehicles. Extending the current 
aerodynamic data base for such a broad range of concepts and issues would constitute a vast research 
endeavor and possibly require more time than is practical. However, focused investigations for selected 
configurations could provide insight to certain fundamental aerodynamic issues in a timely manner. 

The present investigation is directed toward transonic Reynolds number effects for a slender wmg- 
body configuration of the accelerator class. Some discussion of lower speed and lower Reynolds number 
data is also provided for perspective. The accelerator class of configuration tends toward body-dominant 
conical geometries with slender wings. As a consequence, the wing and body related aerodynamics are 
very closely coupled. Some prominent aerodynamic features for this class of configuration include conical- 
like shock structures and boundary layer flows and, at high angles of attack, forebody separated flows 

along with wing (leading edge) vortex flows. , , . 

This research is part of a broader experimental program at NASA Langley. The purpose of this 
program is to (i) design a force-and-moment wind-tunnel model with suitable configuration parametncs 
which is based upon one of the configurational concepts and (ii) examine selected aerodynamic phenom- 
ena over an appreciable range of Reynolds numbers and Mach numbers. The status of this program will 
be briefly addressed. 
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Symbols 


b wing span 

Cq drag coefficient, Drag/ qooS r ef 

Cd,o drag coefficient at zero lift 

Cl lift coefficient, Lift/ q^Sref 

Ci body-axis rolling-moment coefficient, Rolling Moment/ qooS r efh 

beta derivative of body-axis rolling-moment coefficient 
C m pitching-moment coefficient, Pitching Moment/ qooS re fC 

Cn normal-force coefficient, Normal Force/ qooSref 

C n body-axis yawing-moment coefficient, Yawing Moment/ q^Sref b 

C n/} beta derivative of body- axis yawing-moment coefficient 
C mean aerodynamic chord of reference wing planform 

i total body length 

Moo freestream Mach number 

qoo freestream dynamic pressure 

R Reynolds number based on i 

r n nose radius 

S re f area of reference wing planform, extended to model centerline 
a angle of attack, degrees 

ft angle of sideslip, degrees 

$b frustum angle, degrees 

6 C cone angle, degrees 

A leading-edge sweep angle, degrees 

Abbreviations 

LTPT Low Turbulence Pressure Tunnel 

NTF National Transonic Facility 

UPWT Unitary Plan Wind Tunnel 


Configuration and Test Program 

Basic geometric features of the configuration are presented in figure 1. The fuselage was comprised 
of a cone/cylinder/frustum with a cone half angle of 5 degrees, aboattail angle of 9 degrees, and an 
overall length of three feet. The maximum fuselage diameter was 12.87 percent of the body length 
and the sharp nose radius was approximately 0.014 percent of the body length. The delta wing was of 
unit aspect ratio (75.96 degrees leading-edge sweep) with a symmetric 4 percent thick diamond airfoil 
section and a span of 30 percent body length. The leading and trailing edges were sharp. The wing was 
mounted with zero incidence such that the trailing edge fell at 92 percent of the body length. Moments 
were referenced about the quarter chord point of the mean aerodynamic chord for the wing planform 
extended to the plane of symmetry; this occurred at 62 percent of the body length. The vertical tail had 
a leading-edge sweep of 70 degrees, a trailing-edge sweep of approximately -2 degrees, and a symmetric 
4 percent thick diamond airfoil section. Additional details of the model geometry have been reported 
by Fox et al. (reference 1). A photograph of the model mounted in NTF is presented in figure 2. 

The overall range of test conditions for the NTF experiment are summarized in figure 3. Reynolds 
numbers are based upon the reference body length of 3 feet. The tests were conducted for Mach numbers 
ranging from 0.3 to 1.15 and Reynolds numbers ranging from 18 million to 180 million. The maximum 
Reynolds number data were obtained at M^, = 0.6. Test conditions were accomplished with total 
pressures nominally ranging from 2.0 to 7.3 atmospheres and total temperatures nominally ranging from 
120 down to -225 degrees Fahrenheit, The model was sting mounted on an internal six-component force 
balance. The support mechanism included a roll coupling so that pitch and roll could be combined to 
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achieve angle of attack and sideslip. 

A more detailed description of the test program is presented in figure 4 along with the NTF tunnel 
envelope as reported by Fuller (reference 2). The test was structured such that (i) Reynolds number 
effects could be studied at a subsonic and a transonic freestream Mach number and (ii) Mach number 
effects could be studied at a low and a high Reynolds number. Both longitudinal and lateral-directional 
aerodynamic properties were investigated up to an angle of attack of approximately 20 degrees. Sideslip 
derivatives were computed from data taken at +4 and -4 degrees of sideslip. These data were only 
obtained at freestream conditions corresponding to the “corners” of the test matrix shown in figure 4. 
Results for the present paper are focused on the Reynolds number data taken at a freestream Mach 
number of 0.9. 

The data were obtained in NTF with the test section floor and ceiling slotted and the side walls 
solid. The measurements have been compensated for temperature effects, and conventional corrections 
have been applied to the data for the effects of deflection due to load, flow angularity, and base pressure. 
These corrections were, in general, small. No buoyancy corrections have been applied to the data. 
However, these effects were also found to be small. Tests in NTF occurred in early February, 1988. 

The test program for this wind-tunnel model encompasses additional facilities to NTF as shown 
in figure 5. In particular, the model design permits supersonic testing in the Langley Unitary Plan 
Wind Tunnel (UPWT) as well as low-speed Reynolds-number testing in the Langley Low Turbulence 
Pressure Tunnel (LTPT). Included in figure 5 is the tunnel envelope for UPWT as reported by Jackson 
et al. (reference 3), the tunnel envelope for LTPT as reported by McGhee et al. (reference 4), and an 
indication of the freestream conditions at which testing has been completed. Thus far, data have been 
obtained for Mach numbers ranging from 0.2 to 4.5 and Reynolds numbers ranging from 1 million to 
180 million; these results have been obtained with the same wind-tunnel model. Preliminary supersonic 
results from the UPWT investigation may be found in the paper by Covell et al. (reference 5). Results 
from the LTPT experiment have been reported by Fox et al. (reference 1) as well as by Luckring et al. 
(reference 6). 

Both the UPWT and the LTPT tests addressed a substantially broader range of configuration 
parametrics than was investigated in NTF. The configuration variables for the LTPT investigation 
included fuselage nose bluntness, vertical tails, and canards. The UPWT investigation included these 
same variables as well as wing incidence, longitudinal wing position, and wingtip-mounted vertical fins. 

The current test program includes plans for further testing in the UPWT to obtain data at low 
supersonic speeds. In addition, a set of nominally half-scale models have been fabricated for testing at 
hypersonic speeds. 


Results and Discussion 

The general effects of Reynolds number on longitudinal aerodynamic properties are summarized in 
figure 6 for a freestream Mach number of 0.9. As would be expected, Reynolds number had minimal 
effects on the lift and pitching moment data. The lift-curve slope evidences a break at approximately 
4 degrees angle of attack beyond which nonlinear lift effects are observed. The pitching moment data 
show a nose-down break at a comparable angle of attack. These effects are primarily associated with the 
separation- induced leading-edge vortex flow from the wing. The data of figure 6c show a reduction in 
the zero-lift drag coefficient of approximately 25 counts due to an increase in Reynolds number from 24 
to 45 million. The shape of the drag polar was unaffected by this increase in Reynolds number. Further 
increases in Reynolds number had little effect on the drag. 

The results of figure 6c include wave drag increments as indicated by the data presented in figure 7. 
Here the drag coefficient is presented for several freestream Mach numbers ranging from 0.6 to 1.15 at a 
fixed Reynolds number of 90 million. At a freestream Mach number of 0.9, the zero-lift drag coefficient 
has roughly doubled as compared to the results for a freestream Mach number of 0.6; this increment is 
primarily associated with wave drag. Additional discussion of the zero-lift drag rise will be included in 
the section regarding theoretical estimates. In general, the Reynolds number effects for the longitudinal 
forces and moments were nominal. 



Contrary to the longitudinal results, Reynolds number has a more pronounced effect on the lateral- 
directional aerodynamic properties; this effect occurs at high angles of attack. An example is presented 
in figure 8 for the variation of yawing moment with angle of attack at zero sideslip. These data were 
obtained at a freestream Mach number of 0.3 over a range of Reynolds numbers in the LTPT investigation 
reported by Fox et al. (reference 1). All lateral-directional properties in this paper are presented in the 
body axis coordinate system. 

The yawing moment is essentially zero up to a critical angle of attack of approximately 12 degrees. 
Beyond this angle of attack, nonzero values of the yawing moment develop due to asymmetric forebody 
separation and demonstrate a strong sensitivity to Reynolds number. However, the onset angle of attack 
for the asymmetric loads shows little effect of Reynolds number. The initial buildup of yawing moment 
(2 to 3 degrees beyond the onset angle) also shows little effect of Reynolds number. The critical angle of 
12 degrees is on the order of twice the cone semiapex angle, as would be expected from previous forebody 
research such as has been reported by Keener and Chapman (reference 7). These yawing moment trends 
are representative of the other lateral-directional aerodynamic coefficients. These data, along with the 
other results reported by Fox et al. (reference 1), served as precursor information for the high Reynolds 
number investigation in NTF, 

The model configuration for the data of figure 8 differs from the configuration for the NTF tests 
in two respects. First, the sharp nose used for the NTF experiments was a replacement for the one 
utilized for the LTPT test which had become damaged. The second difference is that the vertical tail 
was removed for the data presented in figure 8. 

Reynolds number effects for the current investigation are first addressed by presenting results over 
a range of freestream Mach numbers at both a low and a high Reynolds number test condition, figure 9. 
Before addressing the Reynolds and Mach number effects, it should be noted that the yawing moment 
has the opposite sign at high angles of attack as compared to the results from the LTPT investigation (cf, 
figure 8). This indicates that the flow asymmetry has occurred in the opposite sense. This can be caused 
by either (i) minor differences in the geometry of the nose or (ii) minor differences in flow angularity 
between the tunnels. However, for each test the asymmetry tended to occur either with one sense or the 
other throughout the test; it was very repeatable. 

At a Reynolds number of 24 million (figure 9a) the data show minimal Mach number effects for the 
angle-of-attack range investigated. A lack of sensitivity to Mach number was also observed by Fox et al. 
(reference 1) at a Reynolds number of 9 million for Mach numbers ranging from 0.2 to 0.375. However, 
at a Reynolds number of 90 million (figure 9b) the data do evidence compressibility effects for angles of 
attack in excess of approximately 16 degrees. 

The results presented in figure 9 also demonstrate significant Reynolds number effects at high angles 
of attack. The nonlinear reversal in yawing moment which occurred at a Reynolds number of 24 million 
did not occur at a Reynolds number of 90 million within the angle of attack range investigated. The data 
presented in figure 10 indicate that this change in thebigh angle of attack yawing moment is generally 
associated with high Reynolds number flow. At a freestream Mach number of 0.6 (figure 10a) the data 
for the two lower Reynolds numbers both show the yawing moment reversal whereas the data for the two 
higher Reynolds numbers do not evidence this effect. The transonic case (figure 10b) shows a similar 
trend. In addition, the high Reynolds number yawing moments do not appreciably change beyond 16 
degrees angle of attack. This effect was not observed at Moo = 0.6. It is difficult to determine from the 
data specific Reynolds numbers at which the changes occur. 

The data of figure 10 show limited Reynolds number effects in the 10 to 16 degree angle of attack 
range. This differs from the results presented in figure 8 where Reynolds number sensitivities were 
manifested at only 2 to 3 degrees angle of attack beyond the onset angle of attack for flow asymmetry. 
Therefore, it appears that the angle of attack at which Reynolds number effects become evident in the 
lateral directional coefficients increases as the Reynolds number itself increases. Confirmation of this 
observation will require further testing. 

Sideslip derivative data were obtained at nominally the limiting freestream conditions of the test 
matrix shown in figure 4. The results presented in figure 11 show compressibility effects on the lateral- 
directional stability derivatives at a low and a high Reynolds number. As was observed for the yawing 
moment data of figure 9, the low Reynolds number data (figure 11a) show virtually no compressibility 



effect whereas at the high Reynolds number condition (figure lib) significant compressibility effects are 
indicated for high angles of attack. The results presented in figure 12 indicate that Reynolds number 
effects were limited to high angles of attack and were most prevalent at low speeds. The data of figures 
11 and 12 show that neither Mach number nor Reynolds number had any significant effects on the 
lateral-directional stability derivatives below approximately 14 degrees angle of attack. 

Theoretical Estimates 

A preliminary theoretical analysis of the longitudinal forces and moments was conducted to provide 
design loads as well as to provide some insight to the longitudinal aerodynamic phenomena. Calculations 
were performed with the vortex lattice program of Margason and Lamar (reference 8) as extended by 
Lamar and Gloss (reference 9) to account for separation-induced vortex lift effects by the leading-edge 
suction analogy of Polhamus (reference 10). This method was selected because it has proved over many 
years to provide reasonable estimates of longitudinal forces and moments for a wide range of applications 
as reported by Lamar and Luckring (reference 11), for example. The method was also chosen because 
(i) it tends to provide conservative load estimates (i.e., errors result in over predictions of the loads) 
and (ii) it is a very rapid method to utilize. These attributes are principally due to Polhamus’ suction 
analogy concept which allows nonlinear integral properties associated with leading-edge vortex flows to 
be extracted from a simple linear theory computation. 

Theoretical estimates for the effects of compressibility are presented in figure 13. The normal force 
results are for a fixed angle of attack of 10 degrees whereas the pitching moment results are for a fixed 
lift coefficient of 0.3. Differences between the attached flow theory and the vortex flow theory are due to 
the vortex lift increment predicted for the wing by the suction analogy. Although the trend with Mach 
number is reasonably well predicted by the theory, the magnitudes of normal force and pitching moment 
are not. The differences between the vortex-flow theory and the experiment are larger than would be 
expected from prior experience; they are primarily due to a poor representation of the fuselage in the 
computation as a flat plate. This approach neglects the nonlinear interaction of the leading-edge vortex 
with the thick body. 

A surface grid representation of the configuration (without tail) is presented in figure 14 which 
illustrates the relative size of the body to the wing. Near the forward portion of the wing the body 
thickness will tend to crowd the leading-edge vortex off of the wing. This effect reduces the vortex lift 
increment which also results in a negative pitching moment increment for the assumed moment reference 
point. Methods which properly account for the vortex-body interaction have been shown to accurately 
predict force and moment properties for configurations similar to the one of the present investigation. 
An example has been given by Luckring and Thomas (reference 12) for the wing-body configuration 
tested by Stahl et al. (reference 13). 

Computations for the zero-lift drag rise have also been performed using the analysis system reported 
by Middleton et al. (reference 14). Calculations are presented in figure 15 along with experimental results 
at a Reynolds number of 90 million. The theoretical drag is comprised of a skin friction increment based 
upon the method of Sommer and Short (reference 15) along with a standard supersonic wave drag 
increment; form drag effects were not included in these estimates. 

The computed friction drag provides a reasonable estimate from which the transonic drag rise is 
evident. The experimental drag coefficient at a freestream Mach number of 0.3 is less accurate than 
the other data shown on the figure due to the reduced loads at this freestream condition. This relative 
difference in accuracy is conjectured to be a leading cause for the seemingly high experimental value 
of C D|0 at this Mach number. The supersonic drag estimate is higher than the experimental value 
by approximately 60 counts. A comparable drag increment between theory and experiment was found 
by Compton (reference 16) for the boattail drag of a geometrically similar afterbody when suitably 
normalized. 
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Concluding Remarks 


Selected results have been presented from an experimental investigation in the Langley National 
Transonic Facility (NTF) of a slender wing-body configuration. The tests were conducted at Reynolds 
numbers ranging from 18 million to 180 million based on total model length and at Mach numbers 
ranging from 0.3 to 1.15. The configuration is similar to the accelerator class of vehicles which have 
been considered (along with other configurational concepts) for future high-speed aerospace vehicles. 

Experimental results for the effects of Mach number and Reynolds number on the longitudinal forces 
and moments were found to be nominal. However, the effects of Mach number and Reynolds number on 
the lateral-directional forces and moments were more pronounced. These effects only occurred at high 
angles of attack. Yawing moments became less nonlinear at the high Reynolds number test conditions. 
Compressibility was found to have a larger effect at high Reynolds numbers than was observed at low 
Reynolds numbers. In addition, the angle of attack at which Reynolds number effects became evident 
seems to have increased as Reynolds number itself increases. 

Simple theoretical methods based upon linear theory were found to provide less accurate estimates 
of the longitudinal forces and moments than is usually achievd. This was due to the lack of representing 
the nonlinear wing-fuselage interaction effects as regards the leading edge vortex flow. Approximate 
estimates of the zero-lift drag coefficient were obtained at subcritical and supersonic conditions using 
conventional methodology. 
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Figure 2.- Model mounted in NTF. 
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Figure 3.- Range of test conditions. 
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Figure 4.- Test program for NTF. 
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Figure 6.- Reynolds number effect on longitudinal aerodynamic properties. Mqq = 0.9. 
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INTRODUCTION 

On December 17, 1937, B. Melvill Jones presented the first Wright 
Brothers' lecture at Columbia University in New York (ref. 1). His lecture, 
which was entitled "Flight Experiments on the Boundary Layer," dealt 
specifically with the first British flight observations of transition of the 
boundary layer from laminar to turbulent flow. These data, Jones conclude^, 
showed that it is possible to retain a laminar layer over at least one-third of 
the whole wing surface even when the chord Reynolds number is as high as 8 j 
millions. In the 50 years since this presentation, much flight research has 
been performed to explore the potential of laminar-flow control for drag 
reduction. Both passive control and active control by suction (designated as 
natural laminar flow and laminar-flow control, respectively) have been 
researched and impressive results achieved. The successes of the early natural 
laminar-flow (NLF) flight testing were remarkable, with the achievement of an 
extent of laminar flow and transition Reynolds numbers which were not to be 
exceeded in flight for over 40 years. Nevertheless, mid-century manufacturing 
capabilities were such that insufficiently smooth or wave-free wing surfaces 
led to failure of attempts to transfer this technology to practice. The 
experience with laminar-flow control (LFC) nearly paralleled that of NLF. LFC 
was recognized as a potentially more powerful means for achieving extensive 
laminar boundary layer flow, although admittedly more complex from the systems 
standpoint. LFC flight research began in the 1940's and peaked in the I960' s 
with the USAF/Northrop X- 21 program, the most ambitious LFC flight test to 
date, which attempted to achieve full chord and full span laminar flow on a 
swept wing (ref. 2). Two WB-66 aircraft were fitted with new, full chord 
suction controlled laminar-flow wings and flight tested over 3 years. The main 
result of the program was that laminar-flow control was observed to be 
aerodynamically achievable, but surface qua! ity and structural complexity still 
appeared formidable barriers to LFC applications. 

For a period of almost 10 years, research in NLF and LFC was dormant. The 
energy crisis of the early 1970's revived interest in the technology and flight 
testing resumed. Today, the prospects for a practical technology are brighter 
than ever. We have a greater understanding of the phenomena involved and new, 
less-complex systems concepts are evolving. Critically important to this new 
outlook is the fact that our manufacturing capabilities have dramatically 
advanced since the 1960 ' s and the needed wing surface quality appears within 
our reach. 
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In this paper, the flight testing conducted over the past 1C years in the 
NASA LFC program (ref. 3) will be reviewed. The LFC program has been directed 
towards the most challenging technology application, the high subsonic speed 
transport. To place these recent experiences in perspective, earlier important 
flight tests will first be reviewed to recall the lessons learned at that time. 
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ACRONYMS 


AFWAL 

Air Force Wright Aeronautical Laboratory 

TACT 

Transonic aircraft technology 

NACA 

National Advisory Committee on Aeronautics 

NPL 

National Physics Laboratory 

LFC 

Laminar-flow control 

mac 

Mean aerodynamic chord 

LTPT 

Low-Turbulence Pressure Tunnel 

Max 

Maximum 

PW 

Pratt and Whitney 

NLF 

Natural laminar flow 

LE 

Leading edge 

OASPL 

Overall sound pressure level 

rpm 

Revolutions per minute 

dB 

Decibel 

N-OASPL 

Normalized overall sound pressure level 

Hz 

Hertz 

HLFC 

Hybrid laminar-flow control 


SYMBOLS 

a 

Speed of sound 

S 

Pressure coefficient 

c 

Mean chord 

c 

Chord 

R c 

Reynolds number based on chord 

r t 

Transition Reynolds number 

R x 

Reynolds number based on x 

U 

Free- stream velocity 

M 

Mach number 
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X 


Distance along chord in streamwise direction 


l 

C l 

L 

A 

M 

N 

N 

N 


FAN 


TS 


CF 
P 

S 

t 

GREEK 

P 

P 

A 

V 

P 


Liter 

Lift coefficient 

Lift 

Area 

Fan Mach number 
Engine rpm 

N-factor for Tollmien-Schlichting (TS) calculation 
N-factor for crossflow (CF) calculation 
Pressure 

Surface distance along normal chord 
Time 

Density 

Laminar coefficient of viscosity, or microns 
Sweep angle 
Span station 
Sideslip 
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EARLY LAMINAR- FLOW FLIGHT RESEARCH 


A brief review of some of the most significant past efforts on laminar 
flow is beneficial in understanding the needs for further research. Examples 
of some of the most pertinent flight tests will be presented to highlight the 
knowledge gained. In such a review, wind tunnel tests cannot be ignored 
because of the often elucidative impact they had. 

The earliest known attempts to attain extensive regions of laminar flow in 
flight were made in the late 1930's and early 1940's. Both NLF with favorable 
chordwise pressure gradients and active LFC with boundary layer suction were 
investigated. 

The B- 18 flight test by the NACA in 1939 (ref. 4) was a major milestone in 
the development of NLF. Therein, an attempt was made to prolong the run of 
laminar flow to higher Reynolds number than had previously been achieved by 
flight testing a 17-foot chord, 10-foot span wooden NLF glove on the wing of 
the test aircraft (Figure 1). An exceptional effort was made to. evaluate 
surface quality effects by working the wing to previously unattained smoothness 
and fairness. The flight test clearly displayed the importance of surface 
discontinuities and finish. In fact, the adverse effect of surface 
disturbances (surface waves, two-dimensional type steps, and three-dimensional 
type roughness) was the most pervasive factor observed in the early tests and 
continued to be the principal cause of limited laminar flow in most future 
flight and wind-tunnel tests. Although the severity of the surface 
disturbances was always shown to be aggravated by increased unit Reynolds 
number, it was not until considerable research (made possible by the 
development of the Langley Low-Turbulence Pressure Tunnel , LTPT) had been 
completed in the late 1940' s and 1950' s that an understanding of this 
phenomenon was developed, resulting in a quantitative ability to predict the 
magnitude of permissible three-dimensional disturbances (ref. 5). It was not 
until the 1960 ' s that a quantitative ability to predict permissible waviness 
and two-dimensional type disturbances was developed. 

The maximum transition Reynolds number, based on free-stream conditions and 
distance to the position of transition, attained in the B-18 flight tests was 
about 11.3 million with laminar flow to 42.5-percent chord on an NACA 35-215 
section with only a 3-percent chord loss of laminar flow due to engine and 
propeller noise. This NLF transition Reynolds number was not to be exceeded in 
flight for over 40 years, until the NASA F-ll 1/TACT NLF glove flight tests to 
be discussed later in this paper. The B-18 flight test was very encouraging in 
its time, because it indicated that the flight environment was possibly more 
benign for laminar flow than wind tunnels which until then had achieved laminar 
flow only at lower Reynolds number. The wind-tunnel tests were highly 
compromised because the higher unit Reynolds numbers of the tunnels exacerbated 
the roughness problem. Later wind-tunnel tests, in the quiet LTPT (ref. 6), 
suggested that the B-18 maximum length of NLF was constrained by the glove 
dimensions or achievable aircraft unit Reynolds number, and even higher 
transition Reynolds numbers might be obtained in flight. Indeed, the 
wind-tunnel experiments performed by Braslow (refs. 6 and 7) showed natural 
transition Reynolds numbers of 14 to 16 million for 6-series airfoils. 
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During World War II, several military aircraft were built with NACA 
6-series airfoils, which were designed to achieve extensive natural laminar 
flow. Perhaps the most notable of these airplanes was the P-51 Mustang. But it 
is doubtful that much laminar flow was achieved on these aircraft because 
attention was not given to the surface quality that was required to maintain 
laminar flow. These aircraft flew in a harsh environment for obtaining laminar 
flow (i.e., at high speeds and low altitudes such that the unit Reynolds number 
was high) which placed stringent demands on surface smoothness and fairness. 

But after the war, attempts were made to see if NLF technology could be reduced 
to practice. The flight tests of King Cobra and Hurricane aircraft reported in 
references 8 through 10 are examples of such efforts (see figures 2 and 3). 

The King Cobra used production wing surfaces that were highly polished and 
filled to reduce waviness. The Hurricane employed an NLF section in a special, 
"low-drag construction" wing thought to be suitable for the maintenance of 
laminar flow. With highly polished surfaces good NLF performance was achieved 
on these aircraft, but underlying concerns with the practicality of the wing 
surface tolerances and maintenance defeated these efforts. Now, some 40 years 
later, the general aviation industry is just beginning to explore the use of 
NLF on aircraft for which the Reynolds number capability was more than 
demonstrated by the early NLF flight testing. Many general aviation aircraft 
now fly at higher altitudes, where unit Reynolds numbers are lower, and recent 
advancements in wing fabrication techniques now offer the possibility of 
routinely producing small aircraft with sufficient surface smoothness and 
fairness. 

Active laminar -flow control with boundary-layer suction has been used in 
attempts to extend the laminar flow into the region of adverse chordwise 
pressure gradient, which is not possible to any appreciable degree with NLF. 
Suction through porous materials, multiple slots, and perforations were tried 
with various degrees of success (e.g. refs. 11-15). Three Vampire aircraft 
(figure 4 and refs. 12 through 14) with a number of suction surface 
configurations (continuous porous, perforated, and porous strips) and an F-94 
aircraft (figure 5 and ref. 15) with suction through multiple slots were flight 
tested in the mid-1950's. The F-94 tests were very encouraging. With 69 slots 
between 41-to 95-percent chord, full -chord laminar flow to length Reynolds 
numbers of 36.4 million was obtained on the F-94. The addition of slots and 
suction in the favorable gradient (x/c less than 41 percent) was found to 
significantly broaden the lift-coefficient range for low drag with laminar-flow 
achievement. Laminar flow was lost behind shock waves on the F-94 when the 
aircraft speed was increased to the point where the local Mach number on the 
airfoil surface exceeded about 1.09. An important observation of the F-94 
flight program was that the remains from insect impacts at low altitudes became 
subcritical at high speeds and altitudes above 20,000 feet for the boundary- 
layer flow of the unswept, F-94 wing. Unfortunately, this experience did not 
prevail in later flight test of swept wings, for which smaller critical 
roughness height has been observed in the regions of boundary-layer crossflow 
(ref. 16). The Vampire aircraft tests experienced unusual surface roughness 
difficulties. Continuous suction (from 6 to 98-percent chord) through a porous 
panel cloth (covered with nylon) or through 0.007 inch diameter perforations 
proved nearly as successful as the slotted F-94 surface, but each of these 
surfaces was thought to be impractical to manufacture and maintain. Two, "more 
practical" surfaces were tested, but with poor results. One incorporated 
porous strips of suction with sintered metal inserts; the other had 
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perforations, 0.020 inch in diameter (the smallest holes then thought to be 
practical for manufacturing) . The metal inserts caused surface discontinuities 
under flight loads, and the larger perforations caused transition by 
introducing unstable secondary flow in the boundary layer. 

In the early 1960's, the most ambitious LFC flight program to date was 
im n ky the Northrop Company. Under U.S. Air Force sponsorship, two 

WB-66 aircraft were modified with slotted suction wings and designated X-21 
experimental aircraft (figure 6 and ref. 2). At the end of the program, 
full-chord laminar flow with suction was routinely obtained at Reynolds number 
of about 20 to 25 million. This was only after a long and difficult effort to 
improve performance through the systematic isolation and solution of problems, 
many due to wing sweep. 

The most troublesome phenomenon encountered with the X-21 involved 
leading-edge turbulence contamination, a problem unique to swept wings. On the 
X-21, and at about the same time on a swept slotted-suction wing mounted 
vertically on the fuselage of a Lancaster bomber (figure 7 and ref. 14 and 17), 
the significance of this problem became apparent. Although previous 
small-scale wind tunnel and flight experimentation by the British (refs. 18 and 
^ ki 1 ndi cated the existence of the spanwise turbulence contamination 
problem; its significance had gone unrecognized until the large scale flight 
tests. Subsequent flight and wind-tunnel tests indicated that leading-edge 
scale was a predominant factor and that proper treatment of the inboard wing 
leading edge could prevent turbulence contamination of the swept wing from 
disturbances that propagate down the wing leading edge along the attachment 
line (e.g. ref. 16). Although this phenomenon is now understood, it requires 
careful attention in the design of large LFC aircraft. 

Another adverse effect of wing sweep on the ability to attain laminar flow 
had been found earlier during flights by the British with an AW52 airplane in 
1951 (ref. 18) with a natural laminar-flow airfoil. A series of tests were 
performed where transition was shown to occur very close to the leading edge as 
a result of the formation of streamwise vortices in the laminar boundary layer. 
Later, an inability to obtain laminar flow in the last 20-percent chord of a 
Vampire trainer aircraft (ref. 14) was attributed to the forward sweep of the 
trailing edge. This sweep-induced boundary-layer instability was caused by the 
IdrgG crossflows rosulting from strong, fstvorsblo or ddvorsG chordwisG pressure 
gradients on swept wings. Research prior to the X-21 program (ref. 20) showed 
that the proper application of suction is effective in controlling this 
crossflow instability; a result borne out in the X-21 flight tests. 

Although structural flaws in the X-21 wing design produced surface waves 
and discontinuities that required liberal use of filler for smoothness, 
extensive laminar flow was routinely obtained at cruise altitudes of 40,000 ft. 
at Mach 0.75. A composite of the best wing surface laminar-flow performance is 
shown in figure 8 with remarkably good results obtained on the upper and, to a 
lesser degree, the lower wing surfaces. Flexing of the wing in flight 
continually deteriorated the surface quality due to the filler loss, but a 
series of 12 flights showed good repeatability even with major surface 
discrepancies on the last flight, figure 9. Nonetheless, the X-21 wing 
structure was just not good enough to provide the surface quality needed for a 
convincing demonstration that LFC was ready for application. The poor laminar- 
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flow performance at lower altitudes and higher chord Reynolds number was 
undoubtedly due to the aggravated effects of poor surface quality at higher 
unit Reynolds numbers (figure 10). Still, maximum-length laminar-flow Reynolds 
numbers up to 45.7 million were observed in some areas. 

Another phenomenon adverse to achieving laminar flow was also realized and 
investigated during the X-21 program. It was noted that flight through visible 
cirrus clouds, and sometimes very light haze, caused loss of laminar flow. At 
cruise altitudes, cirrus clouds are composed mainly of ice crystals; 
entrainment of the crystals in the boundary layer produced local turbulence 
leading to the Toss of laminar flow (figure 11). Turbulent vortices shed by 
ice particles in the boundary layer were thought to trigger transition for 
certain combinations of particle size, concentration, and residence time in the 
boundary layer. At the termination of the X-21 program, concerns about this 
phenomenon and other unanswered issues on the operation of LFC aircraft were 
high. 

In summary, when interest in laminar-flow technology was rekindled by the 
energy crisis in the early 1970's, the fundamental aerody n amic concepts of both 
passive and active laminar-flow control had been well established, verified in 
wind-tunnel tests and demonstrated in various flight tests. The aerodynamics 
of the technology appeared to be well in hand. Laminar flow to Reynolds 
numbers up to 16 million had been observed on two-dimensional NLF sections, and 
it was not clear that an upper bound on the transition Reynolds number had been 
reached. Suction control had been demonstrated for boundary layers in adverse 
pressure gradients and on swept wings at Reynolds numbers well above 16 
million; specifically, full chord laminar flow to about 36 million chord 
Reynolds for the former and 46 million for the latter. Yet, no practical 
application had been made with any suction method. The ability to manufacture 
and maintain aircraft surfaces with admissible tolerances^ considerably smaller 
than required for turbulent aircraft, and at acceptable cost was still viewed 
as a formidable challenge. Neither suction slots nor perforations could be 
manufactured economically within required tolerances, and the latter were 
believed to generate disturbances that adversely affected the ability to attain 
large length Reynolds numbers. Criteria for the proper design of slots were 
greatly improved during the X-21 flight program. With respect to perforated 
surfaces, early research indicated the need for hole diameters smaller than 
could be practically fabricated at that time. Porous surfaces with the 
required structural characteristics and aerodynamic smoothness were not 
available. 

Over the past decade, NASA and the aircraft industry have launched 
programs to continue the development of this technology and to provide the 
information needed for objective decisions on its application to new aircraft. 
The flight tests reviewed in this paper have been an integral part of those 
efforts. 
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F— 111/TACT NLF GLOVE FLIGHT EXPERIMENT 


The NACA 6-series airfoils were originally developed for low-drag, NLF 
applications. In actuality, these airfoils were used on many of the early jet 
aircraft because they had very good performance as turbulent airfoils. 

However, modern supercritical airfoil technology has since led to improved 
airfoils with greatly enhanced turbulent performance (i.e., drag divergence 
Mach number, thickness ratio and lift coefficient capability). For this 
reason, in the late 1970's, the Boeing Company designed a new supercritical, 

NLF airfoil in a NASA contract study to evaluate NLF for transport aircraft 
applications (ref. 21). With the Boeing airfoil as a starting point, a new 
supercritical NLF airfoil was designed at the Langley Research Center and 
flight tested at the Ames-Dryden Flight Research Facility on the F— 1 1 1/TACT 
aircraft (refs. 22-24). The objective of the flight test program was to 
investigate natural laminar flow at transonic speeds. 

A supercritical, natural laminar-flow airfoil glove was installed on the 
right wing panel of the F- 1 11/TACT aircraft (figure 12). The glove was made of 
fiberglass skins with an inner core of polyurethane foam and bonded to the 
metal wing skin. For symmetry, an uninstrumented glove was also installed on 
the left wing panel. The glove had a 6-foot span, a 10-foot chord, and was 
finished to "sail-plane" quality. The glove airfoil design pressure 
distribution (figure 13) had a favorable gradient that extended to about 
70-percent chord on the upper surface (dCp/dx/c = -0.4) and to about 50-percent 
chord on the lower surface (dCp/dx/c = -0.8). The airfoil design lift 
coefficient was 0.5 at a Mach number of 0.77 and a Reynolds number of 25 
million. On the upper surface at this condition, supersonic flow extended from 
about 20-percent chord to 70-percent chord where the favorable gradient 
terminated in a weak shock. The glove was installed on the airplane to achieve 
the design pressure distribution at 10 degrees of leading-edge sweep (figure 
12); however, wind-tunnel tests had indicated that the pressure distributions 
at the higher sweep angles (up to 26 degrees) were acceptable for obtaining 
transition data at these conditions (i.e., no leading-edge peaks or premature 
adverse gradients). In hindsight, the low design sweep angle of the glove was 
very conservative, but at that time, some studies had been very pessimistic 
regarding the amount of laminar flow that could be obtained at even moderate 
sweep angles and Reynolds numbers approaching 30 million (ref. 21). 

Results from the flight-test program (ref. 24) indicate that the maximum 
extent of laminar flow was about 55-percent chord on the upper surface at 10 
degrees of sweep for a chord Reynolds number of 28 million. However, as the 
wing sweep was increased to 26 degrees, the transition location moved forward 
to the 10 to 20-percent chord range (figure 12). On the lower surface at 28 
million chord Reynolds number, the maximum extent of laminar flow was about 
50-percent chord (the start of the adverse gradient) and this was achieved to 
sweep angles as high as 15 degrees. 

The wind-tunnel pressure distributions on the glove upper surface were 
much smoother than those obtained in flight (figure 13), particularly at the 
higher sweep angles. Although the majority of the wind-tunnel results have not 
been published, stability analyses are presented for five cases in reference 
25. Based upon the wind-tunnel pressure distributions, these analyses 


67 


predicted transition locations significantly further aft than those measured in 
flight on the upper surface. The irregularities in the flight upper-surface 
pressure distributions, which led to premature transition, were apparently 
caused by shocks propagating onto the glove from the inboard wing and not by 
surface waves in the glove skin. In retrospect, the 6-ft. span of the glove 
was too small to isolate the glove from the flow on the remainder of the basic 
F- 1 1 1/TACT wing over a broad range of conditions. Even for the design point at 
10 degrees of sweep, there was a weak shock wave on the glove near 55-percent 
chord that limited the extent of laminar flow to this point instead of further 
aft near the pressure minimum (figure 13). Since the lower-surface flow was 
subcritical, the lower-surface flight pressures were much smoother than those 
obtained on the upper surface, and in several cases laminar flow was obtained 
to the pressure minimum (approximately 50-percent chord). However, the steeper 
favorable gradient on the lower surface (figure 13) was not suitable for 
achieving large runs of laminar flow at the higher sweep angles because of 
increased crossflow instability. 

The F- 1 1 1/TACT NLF experiment was brief, and consequently the transition 
data were very limited. However, the results were very encouraging. The 
maximum transition Reynolds numbers of about 15 million on the upper surface 
for 10 degrees of sweep, and 14 million on the lower surface at 15 degrees of 
sweep were significantly higher than values obtained in previous NLF flight 
tests. The closest comparable flight test had been conducted over 40 years 
earlier on the B- 18 bomber previously discussed. During that test, a maximum 
transition Reynolds number of about 11.3 million was obtained. 


F-14 VARIABLE SHEEP TRANSITION FLIGHT EXPERIMENT 


Since the F- 1 1 1/TACT NLF glove pressure distributions had not been 
designed to minimize crossflow at the higher sweep angles, and since the 
maximum extent of laminar flow on both the upper and lower surface was 
determined by adverse pressure gradients, even larger transition Reynolds 
numbers at moderate sweep angles seemed possible. In addition, the techniques 
for fabricating and bonding large and very smooth foam and fiberglass test 
surfaces or gloves to metal wings had been developed and proven acceptable for 
flight testing. Consequently, the F-lll/TACT NLF experiment paved the way for 
a follow-on program that could provide a much broader transition data base. 

The F-14 Variable Sweep Transition Flight Experiment was initiated in 1984 
(refs. 26 and 27)) with flight tests being completed in 1987. These tests were 
conducted with an F-14 (figure 14) on loan to NASA from the Navy. Obtaining 
transition data was the primary objective of the program - not airfoil design 
verification. Therefore, only the upper surface of the wing was gloved in 
order to provide a laminar-flow test surface. The gloves extended from about 
10-percent chord on the lower surface to about 60-percent chord on the upper 
surface (spoiler hinge line) and covered the majority of the variable-sweep 
outer panel (figure 15). Four rows of flush static pressure orifices and three 
arrays of hot-films were distributed along the span for determination of the 
local wing pressure distributions and transition locations. These data and the 
other associated flight parameters were monitored in real time on the ground 
during all the testing. 
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Two gloves were flight-tested during the program: one was a "clean-up" or 

smoothing of the basic F- 14 wing (modified NACA 6-series airfoils), while the 
second involved significant contour modifications to the basic F-14 wing. The 
second glove, designed at NASA Langley (ref. 27), provided more moderate 
favorable pressure gradients than the "clean-up" glove, and achieved more of a 
two-dimensional type flow (straighter isobars) over a larger part of the span. 
Both gloves were constructed of fiberglass skins with an inner core of 
polyurethane foam (ref. 28). Measurements taken on the gloves with a 
mechanical deflection gauge having support feet two inches apart indicated wave 
amplitudes no larger than 0.002 in. Representative pressure distributions at 
several Mach numbers are presented in figures 16 and 17 for both gloves. 

The Langley glove design provided a wide variety of pressure distributions with 
different favorable gradients to about 50-percent chord over a broad Mach 
number range (0.6 to 0.8). 

The transition location at 0.7 Mach number on the Langley-designed glove 
is presented in figure 18 as a function of wing sweep for altitudes of 20,000 
ft. and 30,000 ft. Transition locations for the "clean-up" glove are presented 
in reference 27 and for the F-ll 1/TACT NLF glove in figure 12 of the present 
paper. However, to compare various transition or laminar-flow experiments, 
transition Reynolds number is a more appropriate parameter for comparison than 
just transition location. Therefore, the maximum transition Reynolds number 
observed in several of the more significant flight and wind-tunnel experiments 
are presented in figure 19. For this figure, transition Reynolds number is 
based on free-stream conditions, rather than local conditions. In addition to 
the F- 111 and F-14 experiments, included in figure 19 are results from several 
natural laminar flow tests: the B- 18 flight test (ref. 4); the King Cobra 

flight test (refs. 8 and 9); a T-33 flight test (unpublished data from the 
Boeing Company); the 757 NLF glove flight test conducted by the Boeing Company 
(refs. 29 through 31); and low-speed wind-tunnel tests conducted in the 12-Foot 
Tunnel at NASA-Ames Research Center (ref. 32) and the LTPT at the Langley 
Research Center (refs. 6 and 7). Prior to the F— 111 and F-14 flight tests, the 
highest NLF transition Reynolds numbers for airfoils or wings had been obtained 
in the LTPT at Langley and the 12-Foot Tunnel at Ames. These are very quiet 
tunnels and only until recently have airplanes (i.e., jet-powered aircraft) 
been able to match the Reynolds number capability of these facilities. More 
importantly, very few large aircraft have had the capability of providing large 
runs of laminar flow. 

As previously discussed, results from the F-lll/TACT NLF Glove Experiment 
had exceeded the prior maximum values for natural laminar-flow transition 
Reynolds numbers that had been obtained in flight on the B-18 and King Cobra. 
Results obtained during the F-14 VSTFE indicate maximum transition Reynolds 
number values exceeding F-lll/TACT and wind-tunnel values up to 30 degrees of 
sweep. For the F-14 VSTFE, a maximum transition Reynolds number of about 17.6 
million was obtained at 15 degrees of sweep, 13.5 million at 20 degrees, and 12 
million at 25 degrees. Beyond 25 degrees of sweep, maximum transition Reynolds 
number decreased rapidly to about 5 million at 35 degrees of sweep. It should 
be pointed out that for all the maximum transition Reynolds number cases on 
both the F-lll and F-14, the amount of laminar flow was limited by either 
adverse pressure gradient or shock wave location. This suggests that even 
higher transition Reynolds numbers are possible in flight. 
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In comparison to the NLF tests, as would be expected, maximum transition 
Reynolds numbers for most suction or laminar-flow control experiments are much 
higher. As previously discussed, transition Reynolds numbers of about 30 to 36 
million were obtained in flight on the Vampire and F-94, and a value of about 
46 million was obtained in a small area of the X-21 wing. However, with 
suction only in the leading-edge region of swept wings, the transition Reynolds 
number for natural laminar-flow designs can be significantly increased — 
possibly doubled. This concept, called hybrid laminar.flow control (HLFC), is 
discussed later in the paper. 


757 WING NOISE SURVEY AND NLF GLOVE FLIGHT TEST 


In 1985, under a NASA contract, the Boeing Company performed a flight test 
to measure the acoustic environment in cruise on the wing of a 757 aircraft 
with a view towards the determination of the potential effects of the acoustic 
environment on boundary-layer transition (refs. 29, 30, and 31). Prior to this 
flight test, there were no extensive measurements of the noise environment on 
the wing of a commercial transport with wing-mounted, high-bypass-ratio 
turbofan engines. Engine noise concerns had led to conservatism in LFC 
aircraft design studies, with designs restricted to aft engine placement with a 
potentially severe adverse impact on performance and a degradation of LFC fuel 
savings potential. A major part of the 757 flight test was an attempt to 
achieve a limited amount of laminar flow over the wing and measure the impact 
of the engine noise intensity on the extent of laminar flow. Although the 
primary goals differed, an interesting parallel exists between the 757 and the 
B-18 tested some 45 years earlier. As with the B-18, the 757 experiments 
yielded encouraging results with regard to engine noise effects on laminar 
fTow. 


Boeing removed a leading-edge slat on the 757 wing just outboard of the 
starboard engine and installed a 10-foot span, NLF glove constructed of dense 
foam with a fiberglass epoxy overlay to produce a smooth, nearly wave-free 
surfa ce ( fig ure 20). The glove was designed to achieve laminar flow on both 
the upper and lower surface, with 20-to 30-percent chord laminar flow expected 
without adverse engine noise effects. This anticipated result was made 
possible in part by unsweeping the wing to 21 degrees at the glove location and 
by the favorable pressure distribution over the wing. A single microphone was 
installed on the glove leading edge and eight others were installed on each of 
the wing surfaces’ (upper and lower) — three on the glove and five distributed 
over the remaining wing surface (figure 20). Hot films were used to detect the 
transition front on the upper and lower surfaces. Measurements were made over 
cruise altitudes of 25 to 41 thousand feet at Mach numbers of 0.63 to 0.83. 

The starboard engine was throttled from maximum continuous thrust to idle at 
cruise speeds and altitudes. 

Typical measurements of the overall sound pressure level (0ASPL) at the 
design cruise condition are shown in figure 21. Generally, the 0ASPL is lower 
on the upper surface of the wing with measurements ranging from Ul-to 131- 
decibels. On the lower surface, the measured levels ranged from 121- to 136-dB. 
Collectively, the acoustic data presents a rather confusing picture, 
undoubtedly due to the broad range of flight conditions and changing phenonmena 
that influence the OASPL's on the wing. Additional uncertainty is due to the 


acoustic instrumentation because the magnitude of probe interference and self 
noise generation are difficult to assess. However, some important observations 
that have strong implications in laminar-flow wing designs were made. 

Some order in the acoustic data is achieved by normalization of the 
OASPL's with the ambient pressure over the altitude range of the test 
conditions. These data are shown for two microphone locations in figure 22. 
Normalized OASPL's are presented for a microphone on the glove and one aft of 
the glove for various flight Mach numbers and engine power settings; the latter 
reflected by fan exhaust Mach number variations. Flight Mach number and engine 
power setting effects are measureably different on the upper and lower wing 
surface. The lower surface normalized OASPL's show a strong dependence on 
engine power setting with about a 20 dB. increase occurring from engine idle 
(fan Mach number equal to 0.7) up to maximum continuous thrust (fan Mach number 
equal to 1.28) when the lower surface acoustic characteristics seem dominated 
by engine noise. Engine power setting has little influence on the OASPL s on 
the upper surface, but significant variations occur with flight Mach number. 

The wing appears to effectively shield the upper surface from radiated engine 
noise and the dominant noise sources are presumed to be of aerodynamic origin. 
The data present strong evidence that the wing upper-surface flow field has a 
major influence on the radiated acoustic field, particularly at higher cruise 
speeds when shock waves occur on the wing. The supercritical flow over the 
upper surface inhibits forward radiation of sound from downstream sources, 
aerodynamic or engine related. 

Attempts have been made to analyze the 757 acoustic data and make 
comparisons with theoretical predictions (ref. 31). A procedure developed by 
the Lockheed Georgia Company under NASA contract has been used (ref. 33). To 
the authors' knowledge this is the only code available to make near field noise 
predictions that include all the potentially relevant noise sources at flight 
cruise conditions. However, the theory lacks inclusion of the important 
effects of scattering, refraction, and reflection of sound fields due to the 
airframe or flow fields about it. For this reason, predictions for only the 
lower surface OASPL's and spectra have been made. We will not discuss these 
results in any depth herein. Generally, the results indicate that our ability 
to predict the acoustic environment at high cruise speeds and altitudes is 
poor. Theory suggests that the lower wing surface noise should be dominated by 
the fan exhaust broad band shock noise at cruise thrust conditions, which is 
consistent with the observed correlation with the fan Mach number; but the 
predicted levels of OASPL are 10 to 40 dB. too high. Trail ing-edge noise is 
predicted to be an important aerodynamic noise source, particularly at aft wing 
locations; the data doesn't confirm this. Convective and dynamic amplification 
effects have large impacts upon the predictions. These effects or the 
methodology for their implementation are made suspect by the data. Clearly, 
more analyses of these acoustic data are needed to unravel the confused picture 
presented by the data and theory. With further analyses, the broad range of 
conditions for the 757 data could possibly permit useful calibration of the 
Lockheed code. 

The amount of laminar flow obtained on the NLF glove was very encouraging. 
This result indicates that the acoustic environment may be benign enough to 
achieve extensive laminar flow on wings with wing-mounted engines. The results 
are summarized in figure 23 wherein the design condition and conditions of 


71 



maximum extent of laminar flow are shown. A maximum of nearly 30-percent chord 
laminar flow was obtained on both surfaces. At the design condition, best 
results were obtained on the upper surface, athough laminar flow was not 
uniform across the gloved span and was most extensive inboard. The upper 
surface pressures on the glove showed peaks in the outboard region which 
presumably led to earlier boundary-layer transition. Transition was more 
uniform across the lower surface with 26-percent chord laminar flow achieved 
when the aircraft was sideslipped to reduce the leading-edge sweep by 6.8 
degrees. 

On the upper wing surface, the extent of laminar flow was essentially 
unaffected by engine power setting. Since the power setting had no effect on 
the upper-surface noise levels, the unchanging extent of laminar flow is not 
surprising. On the lower surface, however, the noise levels varied over 20 
dB., but almost imperceptibly small (2 to 3-percent chord decreases at most) 
changes in the extent of laminar flow were observed. Over the range of flight 
conditions, boundary-layer stability analysis (ref. 31) identified stationary 
crossflow vortices in the boundary layer to be highly unstable and possibly the 
dominant disturbances leading to transition. The results for the design flight 
condition are typical (figure 24). The crossflow disturbance amplification is 
two to three orders of magnitude greater than the amplification of 
Tollmien-Schl ichting disturbances in both the upper and lower surface boundary 
layers. The small observed effect of variations in engine noise level on the 
transition location on the lower surface may indicate that engine. noise does 
not have a significant effect on crossflow disturbances. If crossflow 
disturbance growth in the leading edge is controlled by suction, laminar flow 
much more extensive than achieved in this flight test could be possible even in 
the presence of engine noise. However, in an HLFC application, the 
Tollmien-Schlichting wave growth may be comparable or greater than the 
crossflow disturbance growth; engine noise might then be expected to limit the 
extent of laminar flow. 


FLIGHT DATA/BOUNDARY-LAYER STABILITY THEORY CORRELATIONS 


Under NASA contract, the Boeing Company has attempted to correlate the 
F- 11 1/TACT and 757 data using linear boundary- layer stability theory (refs. 24 
and 31). Their approach attempts to account for an interaction of crossflow 
disturbances and the Tollmien-Schl ichting disturbances as predicted by linear 
stability theory by cross plotting the amplification N factors for these 
disturbances at transition. Stationary crossflow vortices and oblique 
Tollmien-Schlichting waves (inclined to the streamline at the angle of greatest 
amplification) are considered. A total of 21 flight-test cases were analyzed 
for the 757 NLF glove and are shown in figure 25. Included are points for the 
F-lll/TACT NLF glove flight test that were reported in reference 24. One point 
from a 20 degree swept NLF glove on a T-33 aircraft (unpublished data by the 
Boeing Company) is also included. The F-lll/TACT and 757 data complement one 
another, since the former is mostly for conditions where Tollmien-Schlichting 
waves were dominant, and in the latter, crossflow disturbances were dominant. 
The band enclosing the data is the Boeing recommended transition criteria. 
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This approach has been critized because only stationary crossflow vortices 
and Tol Imien-Schl ichting waves are examined and their interaction is presumed. 
Stationary crossflow vortices have been shown to not always be the most highly 
amplified crossflow like disturbances in the boundary layer (ref. 34). But 
whether or not a particular type of disturbance is present in the boundary 
layer would seem to be the major question. The existence of stationary 
crossflow vortices is well established with observations in many experiments. 
Some authors have reported observations of non-stationary disturbances in the 
crossflow field (refs. 34 and 35); Kohama (ref. 36) suggests that these 
observances are, in fact, evidence of secondary instabilities produced by 
growth of the primary instabilities, stationary crossflow vortices. The 
existence of an interaction of Tollmien-Schl ichting and crossflow disturbances 
is likewise a controversial issue, and the authors of references 37 and 38 have 
proposed that in transition predictions one need only consider the most 
amplified disturbance, with transition occurring when an N factor of around 7 to 
11 first occurs. More analyses of data are required to determine the best 
approach. 

In stability analyses underway at Langley, the F-14, JetStar, and 757 
transition data are being examined with the maximum amplification option of 
reference 39, but these analyses are also including surface and in-plane 
streamline curvature effects (refs. 34 and 40). The initial efforts have 
concentrated on transition data for conditions where Tollmien-Schl ichting wave 
growth is small and crossflowrl ike disturbances dominate the transition 
process. Stability analyses are performed for both stationary and 
nonstationary crossflow disturbances. 

Illustrated in figures 26 through 29 are results for two typical flight 
conditions on the "clean-up" glove of the F-14 aircraft. These flight 
conditions produce strong, favorable pressure gradients (figures 26 and 28) 
that lead to little or no Tol lmien-Schl ichting wave growth and dominance of 
crossflow-like disturbances. The crossflow Reynolds number development 
(figures 26 and 28) is indicative of strong crossflow vortices, and transition 
occurs for both conditions when this parameter exceeds 400, a value somewhat 
higher than the 175 to 300 range observed at low speeds. Previous analyses 
(ref. 41) indicate that the effects of compressibility on crossflow-like 
disturbances are small; comparison of the N factors for stationary and 
nonstationary crossflow vortices, with and without compressibl i 1 ty effects, 
confirms this, as shown in figures 27 and 29. Also, stationary crossflow 
vortices are not the most highly amplified disturbances for these conditions, 
but indeed, nonstationary vortices with frequencies of about 2000-to 3000-Hertz 
are more highly amplified. In the absence of significant compressibility 
effects, the incompressible code developed by Malik and Poll (ref. 34) has been 
used to examine surface and inplane curvature effects on the disturbance 
development. These effects, shown in figures 27 and 29, are quite significant; 
the N factors of the most highly amplified waves are reduced from about 15 to 
around 10 at the measured transition location. Similar calculations are shown 
for two conditions on the 757 NLF glove (figures 30 and 31). Curvature effects 
are again large, and with these effects included, the N factors at transition 
are around 10. Crossflow Reynolds numbers at transition fall in the range of 
300 to 500. 
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The magnitude of the curvature effects in all the flight data analyzed to 
date gives rise to concerns over any previous attempts to correlate transition 
data with stability codes ignoring curvature effects and questions the 
generality of those correlations. Immediate plans are to begin examination of 
data for flight conditions with stronger Tollmien-Schlichting wave growth to 
explore the possibility of interactions of Tollmien-Schlichting waves and 
crossflow vortices. 


THE NASA LEADING-EDGE FLIGHT TEST 


Earlier in this paper, some of the key laminar- flow flight programs that 
laid the foundation for today's knowledge were briefly reviewed. These flight 
tests removed any doubt that extensive laminar flow could be achieved in 
flight. They did not, however, resolve concerns relative to the practicality 
of producing surfaces sufficiently smooth and wavefree, and of maintaining the 
required surface quality during normal service operations. In the late 1970's, 
with the recent significant progress made in the development of new materials, 
fabrication techniques, analysis methods, and design concepts, a reexamination \ 
of these issues appeared warranted. j 

■ . ■ • ; - I 

Previous experience had shown that the leading-edge region of the swept ; 

wing presented the most difficult aerodynamic problems associated with 
attainment of laminar flow. In addition, the leading edge is subject to 
foreign object damage, insect impingement, rain erosion, icing, and other j 

contaminants. Also, an anti-icing system, an anti-contaminant system, and a I 

suction and perhaps purge system must all be packaged into a relatively small 
leading-edge box volume. Most of these problems are common to all the concepts 
under consideration for the achievement of extensive laminar flow, and solutions 
are needed to establish the practicality of laminar flow for various types of J 
aircraft. j 

In 1980, the NASA Leading-Edge Flight Test (LEFT) program was initiated as j 
a flight validation of two leading-edge systems then under development in NASA 
contract efforts with industry. The flight program objectives were to (1) j 

demonstrate that required leading-edge systems can be packaged into a wing 

leading-edge section of a size representative of a commercial transport 
aircraft, and (2) demonstrate systems performance under operational : 

conditions representative of subsonic commercial transport aircraft. Complete i 
LFC leading-edge systems were installed in the leading-edge box of a, Jet Star 
airplane (figure 32). Description s of the systems illustrated in figures 33 j 

and 34 are provided in references 42 and 43. Two leading-edge test articles j 

were built and flown using a perforated and a slotted suction concept. Each 
spanned about 6-foot of the wing and had the same external contour, ; 

dimensionally about equivalent to the leading-edge box of a DC-9-30 at the mean 
aerodynamic chord. The wing leadingredge sweep was 30 degrees. Different. ! 

systems were used in each test articTe. One used suction through approximately j 
1 million, 0.0025-inch diameter, 0.035-inch spaced, electron-beam perforated j 

holes in a 0.025-inch thick titanium skin to maintain laminar flow on the test ; 
article upper surface. A Krueger-type flap served as a protective shield 
against insect impact on this leading edge. In future applications, the 
Krueger shield could also serve as a high-lift leading-edge device. A 
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freezing-point depressant liquid. Propylene Glycol Methyl Ether (PGME) was 
sprayed on the perforated, wing upper surface from nozzles mounted underneath 
the shield to augment the insect shield protection and to provide an anti-icing 
capability. To prevent clogging of the perforations by the wetting fluid, a 
purging system was included to clear the LFC passages by pressurizing the 
subsurface and thus remove PGME fluid from the LFC ducts and surface. The 
second test article used suction through 27 narrow spanwise slots (about 
0.004-inch wide) on both upper and lower titanium surfaces. This test article 
contained anti-contaminant and anti-icing systems consisting of PGME fluid 
dispensed through dual purpose slots in the leading edge. Purge was also 
provided for this leading edge. 

After an initial flight test program to optimize the system's performance, 
the LEFT systems were flight tested in a simulated airline service in different 
geographical areas, seasons, and weather conditions in the United States 
(figure 35). During the simulated service, one-to four-flights per day were 
made from three "home base" airports (Atlanta, Pittsburgh, and Cleveland). A 
total of 62 flights to 33 airports were made. Flights were made from Atlanta 
in July 1985, Pittsburgh in September 1985, and Cleveland in February 1986. 

The weather experienced thus varied from severe summer to severe winter 
conditions. To realistically simulate typical transport operations, an on-off 
operation of all systems was imposed; no adjustments were made prior to or 
during flights. Transport cruise flight conditions were emphasized, but 
investigations were also made of the ability to attain laminar flow at other 
than cruise conditions. References 44 through 47 provide a summary description 
of the program results. 

The emergence of electron-beam perforated titanium as a practical 
manufacturing surface which meets laminar-flow waviness specifications with 
practical aircraft fabrication methods is considered a major development of the 
LEFT program. The perforated titanium leading-edge presented no difficult 
fabrication problems. This test article yielded clearly superior performance 
(relative to the slotted configuration) and was in virtually the same condition 
when flights ended in October 1987, as when flights began in November 1983. 

Four years of flying resulted in no degradation of laminar-flow performance as 
a result of service, and no evidence of any deterioration in surface quality 
was observed. Essentially, compl ete laminar flow on the test article was 
consistently obtained from 10,000 to 38,000 feet altitude with no need for any 
special maintenance. 

The results ’obtained with the slotted-surface test article, however, were 
not as favorable. Fabrication of this configuration involved some extremely 
difficult problems that led to a suction surface that was only marginally 
acceptable with respect to surface smoothness and waviness. This was reflected 
in consistently poorer laminar-flow flight performance than for the test 
article with the perforated surface. Still, as much as 80 percent of the 
slotted upper surface suction area was observed to be laminar in routine flight 
service. 

Since no attempt was made to obtain laminar flow beyond the front spar, 
the LEFT tests should not be interpreted as showing that perforations are 
aerodynamical ly better than slots. Indeed, the perforated approach should be 
pursued with caution because additional flight testing is required to larger 
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values of length Reynolds number. At higher Reynolds number, the experience of 
the early flight tests with larger holes (i.e., progressive performance 
deterioration with increased Reynolds number) could be repeated. Slots may, 
therefore, be preferred at higher Reynolds number. Accordingly, it is clear 
that more development of fabrication techniques for slotted suction surface 
configurations is required; some initial work in this direction has been 
undertaken (ref. 48). 

The LEFT program relaxed concerns about the operational loss of laminar 
flow when entering clouds or haze. It provided some confirmation of an 
extensive analysis of world-wide cloud-cover (based on 6250 flight hours of 
specially instrumented commercial aircraft) which resulted in an estimate of 6 
percent for the amount of flight time spent in clouds and haze (ref. 49). 

During the simulated service flights, measurements were also taken of the time 
spent in clouds and haze. These LEFT results, based on 6 hours and 52 minutes 
of data taken during 13 flights within the United States, showed that clouds 
and haze were encountered about 7 percent of the time (ref. 47). No effort was 
made to avoid cloud encounters, and a sample of one flight including a cloud 
penetration is shown in figure 36. As expected, laminar flow was lost during 
cloud penetrations, but was regained afterwards. The small percentage of time 
that clouds are encountered indicates that laminar flow loss during cloud 
penetrations in cruise will not appreciably decrease the large economic and 
fuel gains predicted for laminar-flow transport aircraft. However, potential 
cloud encounters en route and flight management to avoid clouds could be 
operational considerations for future aircraft. 

To summarize the LFC systems performance during the simulated service, all 
operational experience was positive. No dispatch delays were encountered due 
to the LFC systems. There was no need to adjust suction system controls 
throughout the test range of cruise altitude, Mach number, and lift 
coefficient. Laminar flow was obtained after exposure to heat, cold, humidity, 
insects, rain, freezing rain, snow, ice, and moderate turbulence. The insect 
alleviation systems were required during descent as well as ascent and were 
effective when used. Perforated test article results indicated that the 
supplemental spray system is not necessary for LFC transport airplanes equipped 
with a properly designed insect shield/high-lift device, although the spray 
system may be necessary for anti-icing purposes. Ground deicing of the LFC 
test articles was no more difficult than normal deicing of commercial 
transports, and snow and ice accumulation was easily eliminated using hand-held 
deicing equipment. The NASA LEFT simulated airline service flights 
demonstrated that effective practical solutions for the problems of suction 
laminar-flow aircraft leading edges are available for commercial transport 
aircraft. 

The LEFT program has been the only laminar-flow flight test with suction 
since the X— 21 ended in 1965. The original intent of the LEFT program was to 
examine systems suitable for future laminar-flow control aircraft, but these 
systems would be equally applicable to hybrid laminar-flow control aircraft 
that use suction only in the leading-edge box. 
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HLFC FLIGHT EXPERIMENT 


One of the most significant developments of the NASA research on laminar 
flow in the last' few years has been the recognition of a hybrid laminar-flow 
control concept that integrates LFC and NLF and avoids the objectionable 
characteristics of each. The leading-edge sweep limitation of NLF is overcome 
by applying suction in the leading-edge box to control crossflow instabilities. 
Wingshaping for favorable pressure gradients to allow NLF over the wing box 
removes the need for inspar LFC suction and greatly reduces the system 
complexity. The possibility of achieving extensive laminar flow on commercial 
or military transport aircraft is offered with a system no more complex than 
that already proven in the Leading-Edge Flight Test Program on the NASA 
JetStar. To explore this possibility, the NASA Langley Research . c ® n *®£» 

Air Force Wright Aeronautical Laboratory, and the Boei "9 Co^ercial A 
Company have initiated a cooperative flight program. A high Reynolds number 
HLFC Flight Experiment will be performed on a 757 aircraft equipped with a 
partial-span HLFC system for the upper surface of the left wing. 

The test aircraft and test region are illustrated in figure 37. ^ 20 "! oot 

span of the wing just outboard of the left engine pylon will be modified. A 
new leading-edge box will be installed with suction achieved through a 
perforated titanium surface. The structural concept will be similar to that 
used on the JetStar for the Leading-Edge Flight Test and will include a 
leading-edge Krueger integrated into the full wing high-lift system and 
designed to also be an insect shield for the wing (figure 38). The leading- 
edge 9 box will be contoured to achieve the desired pressure distribution over 
the test surface (figure 39). Analyses indicate that this can be accomplished 
without changing the inspar contour of the 757. Indeed, measurements of the 
757 production wing surface have shown that only minor shaving or fi 9 
some rivets will be necessary to meet laminar-flow smoothness and fairness 
criteria. The inspar production wing surface will thus serve as the t 
surface downstream of the new leading-edge box. 

The HLFC concept, untried to date in flight or in the wind tunne ^? V'jH ^ 
evaluated to chord Reynolds numbers over 30 million at the cruise conditions o 
modern transport aircraft. An extended flight test program is planne 
calendar year 1990 to achieve operational experience with HLFC J nd y 

evaluate the potential for future applications. Success could !ead to the 
long-awaited transfer of this technology to the drawing board and ultimately to 

practice. 


CONCLUDING REMARKS 


The potential benefits of laminar-flow technology have been so enticing 
that possibly no other technology has received such persistent attention in 
flight research over so long a time. The misgivings of the critics are a g 
with the accomplishments of this research. The aerodynamic issues seem nearly 
resolved, and the manufacturing capabilities of the airframe ee 

have advanced to the point that the aerodynamic criteria for smooth, wave-free 
wing surfaces is a practical production goal. The current NASA, AFWAL, a 
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Boeing HLFC Flight Experiment could provide the verification needed to place 
this technology in practice. Initial applications may provide only modest 
improvements, but with the confidence of success, bolder steps could 
revolutionize aircraft design. 
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(CIRCA 1940) 



• NACA 35215 Airfoil 

• Wooden Glove 

C ± 1/1000 inch Waves) 

• Transition Reynolds No. 

~ 11.3 x lO 6 
at 42.57. Chord 

• Engine/ Propeller Noise 

Reduced Laminar Run 
by 37. Chord 



Figure 1. - The B- 18 NLF Glove Flight Test. 



(CIRCA - 1945) 


• NACA 662X116 section 

• Polished and filled production 

wing surface ( i; 1/1000 inch waves) 

• Transition Reynolds No. —llxlO 6 

• Transition beyond C p m j n , 

• Operational evaluations 

Insects 
Dust 8t rain 
Maintenance 

~ Surface quality deterioration 



Figure 2. - The King Cobra NLF Flight Test. 
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Figure 3. - The Hurricane NLF Flight Test 
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Figure 4. - The Vampire LFC Flight Experiments 
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OF POOR QUALITY 


( 1954 - 1957 ) 


• Modified NACA 65 213 airfoil 

• Slotted suction surfaces 

12 slots (41 to 95% chord) 

69 slots (41 to 95% chord) 

81 slots ( 8 to 95% chord) 

« Full chord laminar flow 
06 ^ Moo ^ 0.7 
12.2x10 6 ^ R c 36.4x10 6 

• For M lqcal 1.09, shocks caused 


loss of laminar flow 




Figure 5. - The F-94 LFC Flight Experiments. 


(1963 - 1965 ) 



Step-by-step growth 
in laminar area 
throughout program 



Figure 6. - The X-21 LFC Flight Program. 







Upper Right Wing Standard Suction 
M 00 =O.75 Alt. = 40,000 ft. 



% Chord 


Figure 9. - Maintenance effects on the repeatability of laminar flow on the X-21. 


(10,000 Ft. Altitude) 



Wing Station 


Figure 10. - X-21 high chord Reynolds number results. 
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■■■ Laminar Flow 

Design Turbulent Flow 



• M 00 = 0.73 Alt. = 41,160 ft. • = 0.75 Alt. = 43,120 ft. 

• Large Particles • No Large Particles 

s?150/i., 20/ •Small Particles 

• Small Particles 20 /z or less, 7000/& 

20 M or less, 7000/JL 

Figure 11. - Effect of ice particle encounters on X-21 results. 
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Figure 12. - The F-lll/TACT NLF Glove Flight Experiment. 
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Upper Surface Lower Surface 

Figure 13. - Design and flight pressure distributions for the F- 111 /TACT NLF glove 
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Figure 14. - F- 14 test-bed aircraft for the Variable Sweep Transition Flight 
Experiment (VSTFE). 
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Glove Geometric Characteristics 



Glove Planform Layout 


Figure 15. - F-14 VSTFE glove details. 


A = 20 ° , 77 = 0.68 



x/c 


Figure 16. - F-14 VSTFE "clean-up" glove pressure distributions 
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Figure 19. - Maximum transition Reynolds number for several natural laminar 
flow experiments. 
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Figure 20. - The 757 Wing Noise Survey and NLF Glove Flight-Test Program 


ORIGINAL PAGE f$ 
OF POOR QUALITY 


94 




Mqq = .8, Alt. = 39,000 ft. 



Figure 21. - Overall sound pressure level distribution on the 757 wing at 
cruise conditions. 


N-OASPL = OASPL - 20L0G P/P 40 K 



Upper Surface 
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Figure 22. - Normalized overall sound pressure level versus fan Mach number fo 
the 757 wing. 










N TS -Tollmien-Schlichting Waves 
N Cp -Stationary Crossflow Vortices 



Figure 25. - Boeing stability analyses correlation of the Tollmien-Schlichting vs. 
crossflow N-factors with flight transition data. 



M co = 0.80 
R c = 21 xIO 6 
A =20° 
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V Transition Location 


Figure 26. - F-14 "clean-up" glove pressure and crossflow Reynolds number distributions 
at a Mach number of .80 and a chord Reynolds number at 21 X 10^. 
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Mqo =0.80 R c = 21 x 10 6 C=7.61 ft., A- 20° 
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Fiqure 27. - Crossflow N-factors (with and without curvature effects) for the F-14 
J "clean-up" glove at a Mach number of .80 and a chord Reynolds number 

of 21 X 10 6 . 
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R c = 24x10® 
A =20° 

C = 7.61 ft. 


V- Transition Location 


Figure 28. 


F-14 "clean-up" glove pressure and crossflow Reynolds number distributions 
at a Mach number of .80 and a chord Reynolds number of 24 X 10 . 


98 
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Figure 29. - Crossflow N-factors (with and without curvature effects) for the F- 14 
"clean-up" glove at a Mach number of .80 and a chord Reynolds number 
of 24 X 10 6 
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Figure 30. - 757 NLF glove pressure and crossflow Reynolds number distributions on 
the upper and lower surface at two flight conditions. 
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(Incompressible Theory w & w/o Curvature) 



x/c 

Figure 31. - 757 NLF glove crossflow N-factors (with and without curvature effects) 
on the upper surface and lower surface at two flight conditions. 



Figure 32. - JetStar test-bed aircraft for the NASA leading-Edge Flight Test Program. 
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• Suction on upper surface only 

• Suction through electron-beam- 

perforated skin 

• Leading-edge shield extended for 

insect protection 

• De-icer insert on shield for ice 

protection 

• Supplementary spray nozzles for 

protection from insects and ice 




Figure 33. - The Leading-Edge Flight Test Program perforated test article. 
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Figure 34. - The Leading-Edge Flight Test Program slotted test article. 
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Figure 36. - Typical flight profile from the LEFT Program simulated airline service. 
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Figure 37. - The 757 test-bed aircraft for the Hybrid Laminar Flow Control (HLFC) 
Flight Experiment . 
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Technical Features : 

• Cruise Conditions 

- M«= 0.8 _ 

- R c 30x10° 

- Alt. = 34K to 42K ft. 

• Laminar Flow, 38 to 

62% chord 

• Suction to Front Spar 

• Microperforated 

Titanium Suction 
Surface 


HLFC Panel 



Krueger Deployed 


• All Metal Construction 

• Operational Systems 

- Leading Edge 
Krueger/Insect 
Shield 

- Anti-icing 


Milestones: 

• Contract Award Nov '87 

• Critical Design Review Dec '88 

• First Flight Feb '90 

• Complete Program Nov '90 


Figure 38. - The HLFC Flight Experiment technical features and milestones. 


Cruise Pressure Distribution 



Figure 39. - The 757 wing leading-edge modification for the HLFC Flight Experiment. 


104 



LAMINAR- FLOW WIND TUNNEL EXPERIMENTS 


N 91-24136 


Sy_ 02 - 

/ 1 96 / 


William D. Harvey, Charles D. Harris, 
William G. Sewall , and John P. Stack 
NASA Langley Research Center 
Hampton, Virginia 


SUMMARY 

Although most of the laminar flow airfoils recently developed at the NASA 
Langley Research Center have been intended for general aviation applications, new 
low-drag airfoils have been designed for transonic speeds and wind tunnel 
performance tested. The objective was to extend the technology of laminar flow to 
higher Mach and Reynolds numbers and to swept leading edge wings representative of 
transport aircraft to achieve lower drag and significantly improved operation 
costs. This research involves stabilizing the laminar boundary layer through 
geometric shaping (Natural Laminar Flow, NLF ) and active control involving the 
removal of a portion of the laminar boundary layer (Laminar-Flow Control, LFC), 
either through discrete slots or perforated surface. Results show that extensive 
regions of laminar flow with large reductions in skin friction drag can be 
maintained through the application of passive NLF boundary-layer control 
technologies to unswept transonic wings. An even greater extent of laminar flow and 
reduction in the total drag level can be obtained on a swept supercritical airfoil 
with active boundary- layer control. 
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INTRODUCTION 


Long range market forecasts, operating costs, and strong competition among the 
airlines have provided incentives for the development of new commercial transport 
aircraft with substantially improved fuel consumption and operating costs (ref 
1). One method of improving fuel consumption is through new laminar-flow 
technologies and advanced theories in aerodynamics (refs. 2—4), whereby drag 
reductions are made possible by extending and maintaining laminar boundary layers on 
wing surfaces, thus reducing skin friction drag. The most direct technique for 
maintaining a laminar boundary layer is through the passive, or natural laminar-flow ! 
(NFL) approach which utilizes favorable pressure gradients to avoid Tollmien- 
Schlicting instability that leads to transition. These favorable pressure gradients 
are provided by judicious shaping of the airfoil geometry (refs. 5-9). Maintaining • 

laminar flow by favorable pressure gradient alone becomes increasingly more 
difficult with increasing Reynolds number and wing sweep at transonic speeds; 
however, and some form of active or laminar-flow control (LFC) may be required 
(refs. 2-4, 10, 11). 

Natural laminar-flow airfoils have been developed and extensively applied in [ 

the past on gliders with very low-drag coefficients realized in flight at low-speeds I 
and Reynolds numbers. More recently, carefully designed NASA laminar-flow airfoils ! 
for general aviation aircraft have been developed and wind tunnel tested (refs. 6-9, ! 

12-15) with significantly improved performance characteristics at higher Reynolds 4 

numbers than those in the past. Several of these new low— drag airfoils have been ; 

applied on different general aviation (GA) aircraft. For example, the NLF(1)-04l6 1 

is in production on the Swearingin SX-300 and Freedom Master FM-2 aircraft and the 
NLF(1)-0215F on the NEKO Lancair 210 and Prescott Pusher 200 aircraft (refs. 6, 

7). A Cessna 210 aircraft (refs. 16, 17) has been equipped with a full-span glove 
of the NLF( 1 ) — 04 1 4F (refs. 8, 12) over the original wing for flight tests which 
verified performance predictions. Also, the NLF(1)-0414 has been slightly modified 
and applied on the prototype CIRRUS VK 30 aircraft for flight certification. 


The above airfoils for general aviation application have achieved drag 
coefficients that are lower than conventional turbulent airfoils used for either GA 
or transport aircraft, as illustrated in figure 1. The faired curves are measured 
drag polars for several different airfoils (with flap deflection) at subsonic and 
transonic speeds for a chord Reynolds number of 6 million. The low-speed turbulent 
LS(1)-0417 airfoil is discussed in reference 18 and the 11-percent thick 
supercritical airfoil in reference 19. The drag polars for the NLF(1)-0215F and 
NLF(1)-0414F airfoils show a low drag coefficient "bucket" over a rather wide range 
of lift coefficients which is primarily due to the reduction of skin friction 
through long runs of laminar flow on both the upper and lower surfaces. These long 
runs of laminar flow were promoted by shaping the pressure distributions to minimize 
the amplification of local disturbances that influence transition on the forward 
portion of the airfoil. 
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Conventional ainfoils for transonic aircraft operate at substantially higher 
Mach numbers, Reynolds numbers, and sweep than general aviation aircraft and are 
usually designed for fully turbulent boundary layers which result in higher drag 
values. From figure 1 it can be seen, for example, that the drag obtained on the 
NLF( 1 )-04l 4F airfoil is about 40- to 50-percent lower than either the turbulent low- 
speed LS ( 1 )— 04 1 7 or supercritical airfoil drag over a range of lift coefficients. 

The challenge then, is to extend the technology of NLF and LFC to higher Mach 
numbers, Reynolds numbers, and sweep with low drag. The performance verification of 
such laminar-flow concepts on airfoils at transonic speeds requires wind tunnels 
with excellent flow quality environments. 

This paper presents experimental wind tunnel results and analysis for two 
advanced NASA laminar-flow airfoils designed for high subsonic and transonic flow 
conditions . 


SYMBOLS 


c 

C, 


Wing or airfoil chord 

^ . section drag force 

Drag coefficient, a 

Q c 

CO 


Lift coefficient, 


sect ion lift force 
q c 


P " P„ 


Pressure coefficient, 


Suction coefficient, 


(pw) v 

(pu) 


frequency 


M 

N 

R 


Mach number 

Natural logarithm of ratio of distance amplitude to its value at the 
neutral stability point 

pU 
r 00 

Reynolds number per unit length, — — 


pU c 
r 00 

R Reynolds number based on chord, 

c y 

t Thickness 

x Chordwise distance measured from leading edge 

a Angle of attack 
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p 


Density 


6 f Flap angle 

A Sweep angle 

y Viscosity 

Subscripts : 
c Chord 

CF Cross flow- 

max Maximum value 

min Minimum value 

N Normal to leading edge 

s Suction 

t Total 

tr Transition 

TS Tollmien-Schlichting 

w Wall or wake 

» Free stream 
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LOW-DRAG AIRFOIL TESTING 


This section presents a review of the mechanisms which cause transition; 
comments on the ability to properly simulate flight conditions in a transonic wind 
tunnel; and discusses recent advances made in diagnostic instrumentation as appli 
to low-drag airfoil testing. 

One of the most important parameters to be determined during performance 
verification of low-drag airfoils is the location of boundary-layer transition. 
Equally important, however, is understanding the mechanisms which cause 
transition. Transition depends on many factors, including Tollmien Schlic 1 g 
instability (TS), cross-flow instability (CP), attachment line contamination (AL) , 

and Gortler vortices (GV). 

The boundary- layer profile in the direction of the external flow is very 
similar to that of a two-dimensional boundary layer and by a process of receptivity, 
smail disturbance waves (TS) are developed. These small disturbances grow initially 
through a linear process, followed by a nonlinear growth and breakdown to 
transition For unswept wings, TS disturbance waves are the dominant cause of 
transition and adverse pressure distributions and increasing Reynolds number have a 
major influence in their growth. 

Wind tunnel experiments on swept wings (ref. 10) have. shown that transition may 
occur near the leading edge due to the instability of cross-flow velocity profiles 
the boundary layer. Cross-flow (CP) velocity profiles exhibit inflection points near 

the surface, and boundary-layer stability theory shows ^ ^oss-noi . 

lead to transition (refs. 20-24). The major parameters that influence 
instability are wing sweep, pressure distribution, and Reynolds numbers. 

Spanwise turbulent contamination, induced by instability of the 
juncture boundary-layer or leading-edge roughness, can propagate along the leading 
edge of swept wings (ref. 25). The occurrence of such a disturbance along the 
leading-edge attachment line (AL) could result in the entire wing surface boundary 
laver becoming fully turbulent for attachment line momentum thickness Reynolds 
numbers 6 exceeding about ,50. Attachment line transition is dependent - «ng sweep 
angle, leading-edge radius, Reynolds number, and flow acceleration in and along the 

attachment line region. 

on the lower surfaces of cambered airfoils, GSrtler vortices (GV) may develop 
as a result of a centrifugal instability of the flow over regions of concave surface 
geometry (refs. 26-29). Such disturbance vortices are counter rotating pairo 
lift off the surface and travel streamwise while remaining in the boundary layer. 
Gortler vortices are normally highly unstable and quickly lead to transition, 
however , V their growth may be damped when traveling from concave to convex curvature 

regions - 

All of these instabilities are very important (refs. 30-32), and particular 
attention must be given to avoiding these phenomena when designing and performanc 
testing swept laminar-flow airfoils. 


LAMINAR-FLOW TESTING IN TRANSONIC WIND TUNNELS 

There are certain limitations in transonic test capabilities which must be 
considered. Since low-drag airfoil testing is fundamentally a boundary ayer 
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stability problem, the ability to properly simulate flight conditions in a wind 
tunnel can be strongly affected by Reynolds number scale (ref. 33) » environmental 
disturbance levels (refs. 34-36), model surface conditions (refs. 37-41), and tunnel 
wall interference (refs. 42-44). Also important are the methods and instrumentation 
, used to measure disturbance levels and boundary-layer transition characteristics. 


DIAGNOSTIC INSTRUMENTATION 

The test results obtained with any model is a direct function of the test 
section flow characteristics. Transonic tunnels are, in general, not well 
calibrated in terms of free-stream dynamic quantities, but pressure, density, and 
temperature fluctuations are very important flow disturbance sources and must be 
known to properly evaluate the turbulence levels (ref. 36). Hot-wire anemometers 
were used to measure__mass flow fluctuations from which components of velocity 
fluctuations (u, v, w) have been derived. The acoustic environment was also 
measured using either pressure transducers or microphones. 

Another type of anemometry consists of individual thin-film sensors, which may 
be attached directly to the model surface or flush mounted (ref. 45). it operates 
on principles similar to the hot-wire anemometer and is used to detect transition. 
The single-type film sensors (ref. 45) have recently been configured in a multiple 
sensor array (ref. 46) to simultaneously detect separation and reattachment 
locations or shock- boundary layer phenomena. More details about advancements in 
instrumentation for transonic testing are presented in reference 47. 

Figure 2 shows the locations of individual thin-film sensors which were 
deposited on the ends of tiny quartz rods and mounted flush with the surface of a 
model. The sensors were maintained at 80 °C above local recovery temperature by a 
constant temperature anemometer. The installation technique for the hot-film 
sensors on this model is given in reference 45. Examples of output RMS voltage 
traces and their corresponding boundary- layer states are also shown in figure 2. 

The amplitude of the fluctuating sensor voltages indicates the relative state of the 
boundary layer based on assigned arbitrary values of intermit’terncy from 1 to 10. 

The low amplitude dynamic signal indicates the lower shear stress that is 
representative o ? a laminar boundary layer (intermittency factor = 1). The higher 
amplitude dynamic signal indicates the higher shear stress that represents a 
turbulent boundary layer (intermittency factor = 10). The transitional boundary 
layer has both the laminar and turbulent type signals (intermittency 
factor - 3 to 8) , and it is the relative amount of each, or intermittency factor, 
which indicates the degree of transition. 

Figure 3 shows an example of the Multielement Dynamic Shear Stress Sensor 
(MEDS ) technique (ref. 47) that was used to detect shock-induced separation on the 
upper surface of an airfoil at transonic speeds. The MEDS 3 consists of a large 
number of closely spaced (0.1-inch) individual Nickel films that were vacuum 
deposited on a thin substrate of thickness 0.000010 inches in a straight line 
array. The substrate is bonded to the model surface with sensor array oriented 
streamwise and leads removed at the model spanwise ends. The figure illustrates 
that shock induced separation can be clearly detected by the phase reversal 
phenomenon observed by P. Stack and S. Mangalam (ref. 46) in low-frequency dynamic 
shear stress signals. When the shear stress fluctuations from adjacent surface 
thin-film sensors are correlated, a zero phase is observed upstream of the shock 
(fig. 3(a)). A -180° phase shift is observed across the shock-induced separation 
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(fig. 3(b)); a +108° phase shift is observed across reattachment (fig. 3(c)); ana 
the signals are again in phase downstream of the shock- induced separation and 
reattachment (fig. 3(d)). 


TRANSONIC FACILITIES 

A brief description of the transonic wind tunnels used to obtain the data 
presented in this paper follows. 

8-Foot Transonic Pressure Tunnel 


The NASA Langley 8-Foot Transonic Pressure Tunnel (8-Ft. TPT) is a closed- 
circuit single-return variable density continuous-flow wind tunnel with a 
contraction ratio of 20:1. The test section walls are slotted (5-percent porosity) 
top and bottom, with solid sidewalls fitted with windows for schlieren flow 
visualization. In 1981, the facility was modified for flow quality improvements and 
reconfigured for testing of a large-chord swept laminar— flow control (LFC) airfoil 
at transonic speeds (refs. 11, 48). A honeycomb and five screens were installed in 
the settling chamber to suppress the disturbance levels in the test section. A 
contoured liner was installed on all four walls of the test section to simulate 
free-air streamlines about an infinite yawed wing. This contoured liner produces a 
contraction ratio of 25:1 and covers the floor and ceiling slots. An adjustable 
sonic throat (0.78 < < 0.84) is located at the end of the test section to block 

upstream propagation of diffuser noise. The honeycomb and screens were installed as 
permanent additions to the facility and the liner as a temporary addition, to be 
removed at the end of the LFC experiment (ref. 48). 

The combination of honeycomb, screens, and choke provide very low disturbance 
levels of p/p and u/u = 0.055 percent, in the test section at M o = 0.8 and 

r oo oo 

R/ft » 3 x 10 6 , based on conventional single hot-wire measurements and data 
reduction techniques. However, results also obtained using a new and more accurate 
three-wire, hot-wire probe technique at the same test conditions indicated a 
disturbance level of u/u,,, = 0.3 percent (ref. 47). 


OHIO STATE UNIVERSITY 6- BY 22-INCH TRANSONIC TUNNEL (6 x 22 TT) 

The Ohio State University 6- by 22-Inch Transonic Tunnel is a closed-circuit 
blowdown type facility (ref. 49). The test section is 6-inches wide and 22-inches 
high, has plane parallel sidewalls, perforated walls on the top and bottom with 
separate, self-aspirated plenums, and an array of bars in the exhaust to choke the 
flow and fix the test section Mach number. 

Air is supplied from a 1500 cubic foot storage tank at 2600 psi through two 
preset control valves, enters the settling chamber through a distribution device and 
a perforated bulkhead, passing through a honeycomb and six screens, and enters the 
test section through a three-stage contraction region. Models are mounted in a 
circular sidewall port , which may be rotated for angle-of-attack from 0° to ±180°. 
Model pressure distributions are acquired by scani valves, and Kulite gages are used 
for diagnostic testing. Wake profiles for the determination of drag are acquired 
by means of a single traversing probe at transonic speeds or by multiple probes at 
low speeds. 
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The maximum stagnation pressure is 4.423 atmospheres with a Mach number range 
from 0.2 to 1.1. The operational unit Reynolds number range is dependent upon Mach 
number. For example, 


M 

R/ft 


0.2 

2 to 

7 x 10 6 

.4 

4 to 

14 x 10 

.8 

7 to 

23 x 10 


RESULTS AND DISCUSSION 
Unswept High-Speed NLF Airfoils 

High subsonic cruise speeds provide a more difficult design challenge than Tow 
subsonic cruise speeds when laminar-flow wings are used. At high speeds, low sweep 
is desirable to minimize cross-flow boundary-layer instabilities and leading-edge 
contamination; high sweep is desirable to minimize wave drag and delay transonic 
drag rise. Consequently, a compromise must be made to prevent boundary-layer 
instabilities from causing transition near the leading edge and, at the same time, 
not reduce the drag-rise Mach number to an intolerably low value. Since thinner 
airfoils rather than reduced sweep are an alternate way of delaying drag rise, 
airfoil/wing thickness will also be a variable. Another negative factor in high- 
speed design is the larger Reynolds numbers that are inherent and the adverse effect 
they have on all instability mechanisms. Careful tailoring of the airfoil contour 
and the associated pressure distributions is an important activity at low speed but 
at high speed it is a "make-or-break" one. Favorable gradients required for laminar 
flow drive the design C^'s down and the wing area up. The need for highly accurate 
inviscid and viscous flow codes as well as boundary-layer stability routines is 
clear. 

Based on the encouraging results obtained by geometric shaping to achieve 
extensive laminar flow on both surfaces of the low-speed NLF(1)-0414F airfoil in the 
wind tunnel (ref. 12) and in-flight (ref. 16), effort has been recently directed 
towards extension of the concepts to higher speed NLF airfoils. Details of one of 
the resulting two-dimensional concepts, HSNLF(1 )-021 3F> are reported in references 
13 and 15. Wing body integration of the airfoil is discussed in reference 13 and 
design of an integrated trail ing-edge flap is presented in reference 14. 

The HSNLFO )-0'21 3F, has been fabricated of metal and tested (ref. 15) in the 
NASA Langley LTPT and 6- by 28-Inch TT to investigate its low-speed high-lift 
characteristics and high-speed drag-rise characteristics. The airfoil was designed 
for a lift coefficient of 0.26, Mach number of 0.70, Reynolds number of R = 1 1 x 
10 6 , and t/c = 0.134. This particular airfoil was designed for application to an 
essentially unswept wing but could be used for a swept wing as well. 

The HSNLF ( 1 )-021 3F airfoil profile and design pressure distribution are shown 
in figure 4. The design shockless "pressure distribution was expected to provide 
laminar flow rearward to about x/ c - 0.55 on the upper and x/c = 0.70 on the lower • 
surfaces up' to R c S 10 x TO 6 . The upper surface forward of 50-percent chord employs 
the NLF(1 )-04l 4F airfoil geometry; the upper surface aft of the 50-percent chord was 
modified from that of the NLF( 1 ) — 04 1 4F to minimize the possibility of turbulent 
separation (refs. 2, 13, 15 ). The bottom side of the nose was slightly modified to 
minimize off-design pressure peaks on the lower surface. At the same time, the 
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pressure distribution at the design Mach number and lift coefficient remains 
essentially subcritical, thereby avoiding recompression shock waves. 

Initially, the HSNLF( 1 ) — 02 1 3F airfoil was tested both at low speeds 
(M £ 0.3) and high speeds (0.35 £ M £ 0.80) and the results are presented in 
reference 15. Near the design lift coefficient, 0.23 £ C £ 0.25, the low-speed 
tests revealed both low-drag coefficients (0.0038) and transition locations with 
individual surface thin-film sensors at x/c = 0.50 on the upper surface and 
x/c = 0.60 on the lower surface. However, at the design Mach and Reynolds numbers, 
the high-speed test showed much higher drag coefficients (0.0078 to 0.0083) at the 
design Mach and Reynolds numbers. This higher drag level was attributed to 
premature transition caused by poor tunnel flow quality in the NASA Langley 
6 x 28 TT. 

To investigate the 1 am i nar,- f 1 ow potential of the HSNLF( 1 )-021 3F airfoil, 
another wind tunnel test was conducted on the same metal model that was used in the 
NASA Langley 6 x 28 TT, with multielement thin-film sensors to detect transition, 
separation, reattachment, and shock location (ref. 46). The 6- by 22- inch transonic 
wind tunnel at the Ohio State University (ref. 49) was used for this experiment. 

Test requirements and facility capabilities limited the maximum chord Reynolds 
numbers to 5.0 million. Limited test results were also obtained on a composite 
model in this same facility. 

Some of the results from this experiment are shown in figures 5 and 6, where 
comparisons have been made between the experimental and theoretical pressure 
distributions. Included with the pressure distributions are the locations for the 
beginning, peak, and end of transition as indicated by the thin-film sensors on the 
model's upper surface that extend from 0.3 < x/c < 1.0. The theory (ref. 50) has no 
laminar boundary layer; the start of the turbulent boundary layer was arbitrarily 
set at 0.40c on the upper and 0.50c on the lower surfaces. 

The shock-free experimental and theoretical pressure distributions presented in 
figure 5, at the design Mach number and lift coefficient agree very well up to the 
recompression zone for both the upper and lower surfaces. The theory then shows a 
stronger pressure recovery than the experiment. Good agreement exists between the 
measured and predicted drag of = 0.0057. No wall interference corrections were 
made to the experimental data. 

The hot-film sensors (refs. 46, 47) on the upper surface indicate the beginning 
of transition at x/c = 0.50, near minimum pressure, followed rapidly by peak 
transition at x/c = 0.545, then fully turbulent flow developing at x/c = 0.593- No 
boundary- layer separation was indicated by these upper surface sensors. It should 
be noted that, while the measured transition location agrees with the design 
objective, the results are limited to a Reynolds number of 4.2 million instead of 
the design 11 million. 

Figure 6 shows another experimental and theoretical pressure distribution 
comparison at an off-design Mach number of 0.77 which results in a strong shock wave 
on the upper surface, as indicated by both experiment and theory (ref. 50). 

However, the experiment shows that the beginning of the shock is approximately 0.10c 
ahead of the theoretical location, along with significantly reduced suction levels, 
followed by an entirely different pressure distribution than shown by theory. It 
should be noted that the theory (ref. 50) has no provision for shock-boundary layer 
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interaction, and there were no wall interference corrections made to the data. Thin- 
film sensors in the region of the shock wave indicated an attached laminar 
separation bubble just downstream of the beginning of transition, which occurred at 
the shock wave itself (fig. 3). Therefore, the laminar boundary layer extends 
rearward to the shock wave which interacts with a laminar separation bubble and 
causes it to shift forward. The experimental pressure distribution following the 
separation bubble, i.e., after the reattachment point, develops a positive gradient. 

Figure 6 not only shows large differences between measured and predicted 
pressure distributions , but also between drag coefficients that are caused by ; 

interaction of the laminar boundary layer and shock wave that are not accountable by 
the theory. For the design case (fig. 5), good agreement exists between theoretical i 

predictions of pressure distribution and drag since there is no shock or boundary- 
layer separation. f 

I 

Figure 7 shows a comparison of the transonic drag-rise with increasing Mach 
number obtained on the same metal model of the HSNLF( 1 )-021 3F airfoil in two wind 
tunnels. Also shown for comparison is a single-drag value at M = .7 obtained on a 
composite model of the HSNLF( 1 )-021 3F in the O.S.U. 6 x 22 TT. The angle of attack 
was set at zero and the resulting lift coefficients ranged from 0.26 to 0.28. From 
figure 7, it appears that measured and predicted low to moderate drag coefficients ‘ 

(fig. 5), ranging from about 0.0057 to 0.0065, can be attained on the metal model up ; 

to a Mach number of 0.73 at R c = 4.2 x 10 6 . Apparently, much lower drag values of j 

about 0.0045 are possible at the design M = 0.7 with the smooth composite model I 

compared to that obtained for the metal model in the same tunnel. With transition ! 

fixed at 0.05 chord, the drag coefficient increases to 0.0083 while the drag-rise 
Mach number is reduced to about 0.715. - 


SWEPT SUPERCRITICAL LFC AIRFOILS 

The concept of combining geometric shaping and boundary-layer control through 
suction on airfoils to achieve very low drag dates back to the late 1930's when the 
feasibility of achieving full-chord laminar flow on conventional swept wings, with 
suction applied through many closely spaced surface slots, was established by 
Pfenninger (refs. 10, 51). Results were obtained on large-chord wing sections 
(modified NACA-66012) of 30° sweep and t/c * 0.12 in three different wind tunnels 
(ref. 10) and flight (refs. 37, 39). Flight tests on the X-21 and F-94 aircraft 
(refs. 52, 53) achieved full-chord laminar flow to R c = 22 x 10 6 at M o = 0.80 and 
R = 36 x 10 s at M = 0.72, respectively. The X-21 wing was swept 33° and the F-94 
wfng was unswept. °°Since this research demonstrated the potential for significant 
drag reduction through application of boundary-layer control, interest was generated 
in evaluating the feasibility of combining boundary- layer control with supercriti cal 
airfoil technology at conditions which are typical of high-performance transports. 
About 1975, Dr. Werner Pfenninger proposed such an experiment and a large chord 
swept supercritical airfoil, with suction capability on both the upper and lower 
surfaces, was designed, constructed, and tested in the NASA Langley 8-Ft. TPT (fig. 
8). This NASA airfoil is designated SCLFC ( 1 )-051 3F. Details of the airfoil and 
suction system design along with the test setup have been reported (ref. 11). One 
of the major objectives of the experiment was to evaluate the feasibility of slotted 
and perforated suction surface concepts for laminarization in a large supercritical 
zone. 
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As mentioned above, requirements for this test included modification of the 
wind tunnel to achieve the desired flow quality and test section wall contouring to 
simulate free air flow about an infinitely yawed model at the design Mach number 
(refs. 44, 48). 

Slotted and Perforated Models 


The LFC airfoil model has an aspect of about 1.0. It consists of a wingbox 
structure to which three upper and three lower surface suction panels and a 
trail ing-edge flap are attached (fig. 9). Figure 10 is a sketch showing the slotted 
suction surface with the internal airflow metering and ducting system. Laminar flow 
control by boundary-layer removal on the slotted configuration was achieved with 
suction through closely-spaced spanwise slots (fig. 10) on the airfoil surface. 

After passing through the slots and small underlying plenums, the air passed through 
appropriately spaced metering holes and was collected by spanwise ducts of constant 
cross section with suction nozzles located at the ends. Air from the nozzles passed 
through model evacuation lines, through airflow control boxes which controlled the 
amount of suction to the individual duct nozzles, through variable nozzles, through 
hoses to a collector manifold and, finally, to a 10,000 CFM compressor with a 4.5:1 
compression ratio which supplied the suction. 

In the upper surface midspan region where full-chord laminar flow was expected, 
the suction slots extended rearward onto the flap to x/c = .96. On the lower 
surface, suction extended rearward to x/c = .84. Outside of the laminar midspan 
region, the upper surface spanwise slots extend rearward only to the flap hinge line 
at x/c = .89. A photograph of the swept LFC airfoil installed in the test section 
of the 8-Ft . TPT is shown in figure 11. Figure 11(a) is a downstream view looking 
at the swept leading-edge upper surface, and figure 11(b) is an upstream view of the 
trailing edge. Suction ducts to control the boundary layer in the model/tunnel 
junctures may be seen in the liner at the floor and ceiling (fig. 11(a)). 

Figure 12 is a sketch of the swept LFC airfoil model illustrating replacement 
of the slotted upper surface with three perforated suction surface panels. The 
perforated panels were designed and fabricated by McDonnell Douglas Aircraft Company 
under contract to the NASA Langley ACEE Project Office (ref. 54). The original 
bottom surface slotted panels were retained for this configuration. Figure 13 is a 
sketch showing the perforated suction surface and internal airflow metering and 
ducting system. This surface was fabricated with fiberglass flutes or hat sections 
with a continuous perforated skin bonded to the outer part of the hat sections 
resulting in alternating suction and nonsuction areas. These areas extend spanwise 
and are constant in width. The active suction and nonsuction area widths are 
nominally 0.60-in. and 0.3~in. wide, respectively. The perforations in the surface 
skin, produced by electron beam drilling, are approximately 0.0025 Inches in 
diameter on 0.025-inch spacing, yielding an open area porosity ratio of 0.8 
percent. Beneath each suction strip are spanwise flutes and ducts that are 
connected through a number of metering holes. Calibrated suction nozzles are 
provided in each duct and are connected via hoses to the suction compressor through 
the same control used for the slotted model. The perforated suction surface extends 
from 0.02 £ x/c £ 0.89 along the entire span. The turbulent wedge regions at the 
ends of the model are very similar to those for the slotted model and the suction 
applied and controlled is also similar. 
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Experimental Results 


Figures 1 4 and 15 show a comparison of the measured and theoretical chordwise 
pressure distributions at midspan and = 0.82 for chord Reynolds numbers of 

10 and 20 x 10 6 . Figure 14 shows the measured pressure distribution for the slotted 
model to be essentially shock free at R c = 10 x 10 6 with full-chord laminar flow (as 
evidenced by the surface thin-film gages used for detecting transition) on both 
upper and lower surfaces (see figs. 16, 17). The higher and less uniform than 
design velocities on the upper surface were attributed to classical problems 
associated with wind tunnel testing — model deformation and wall interference. Local 
Mach numbers measured on the test section wall opposite the upper surface indicate 
that the supersonic bubble on the airfoil was larger than design and extended 
practically to the wall. The larger bubble was believed to be due to the inability 
to account for three-dimensional boundary-layer displacement thickness effects with 
the accuracy desired in the design analysis of the contoured liner wall 
(ref. 44). 

The pressure distribution for the perforated model at R c = 10 x 10 s (fig. 14) 
is very similar to that of the slotted model except for slightly higher velocities 
on the upper surface and with a very weak shock appearing at about x/c = 0.70. Even 
with the higher velocities and weak shock on the perforated upper surface, full- 
chord laminar flow was achieved (fig. 18) for R c = 10 x 10 6 . These higher upper 
surface velocities are related to the more positive 0.71° angle of attack required 
for the perforated model as compared to the 0.51 ° angle of attack for the slotted 
model. The higher angle of attack was required to minimize the strength of the 
shock at the end of the supersonic zone as discussed below. 

Measurements under simulated loading conditions indicated that the perforated 
model deformed more in the chordwise direction than the slotted model. Because of 
the larger deformation, the rear panel of the perforated model was closer to the 
test section wall opposite the upper surface than was the slotted model when both 
were positioned at ot = 0.51 This resulted in a saddle back shaped upper surface 
pressure distribution with a strong shock near x/c = 0.70 and a sonic bubble which 
impinged on the wall. Since the angle of attack was adjusted by rotating the model 
about x/c = 0.24, increasing the angle moved the trailing edge of the model away 
from the wall, weakening the upper surface shock and reducing the extent of the 
sonic bubble. 

In order to reduce the velocities nearer to design, it was necessary to modify 
the perforated model mounting blocks and translate the model 0.25-inch further from 
the wall above the upper surface. With this translation, the angle of attack could 
be reduced to 0.71° without strengthening the upper surface shock. At these 
conditions, the upper surface pressure was very near that for the slotted model and 
the extent of the sonic bubble was about the same for the two models. 

With increases in Reynolds number, transition moved forward gradually on the 
upper surface and rapidly on the lower surface (figs. 16— 18). Figure 15 compares 
the resulting pressure distributions at R c = 20 x 10° and indicates very similar 
distributions for the two models except for the previously mentioned higher 
velocities on the upper surface of the perforated model and the tendency for the 
perforated model to show a developing shock on the lower surface near midchord. The 
laminar boundary layer on the lower surface (both models used the same three slotted 
lower surface panels, fig. 9) was unable to withstand the adverse pressure gradient 
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leading into the trailing edge eusp region above approximately R Q = 14 x 10® where 
transition moved rapidly forward and the flow separated in the rear cusp. With 
decreasing extents of laminar boundary layer and the appearance of separation on the 
lower surface, the local effective area distribution of the test section changes, 
resulting in higher free-stream Mach numbers being required to achieve the design 
plateau pressure distribution as Reynolds number increased. A loss in model lift 
and increase in drag also occur with decreasing extents of laminar flow. 

It should be noted that, in general, as the free-stream Mach number increases, 
the corresponding static pressure decreases rapidly, causing the local differential 
pressure drop across a given suction surface to become large with corresponding 
increased suction requirements. This trend may become a challenge for designers of 
high transonic-supersonic LFC transports. 


Analysis of spanwise pressure distributions and transition patterns at M degign 
showed that the flow was nearly two dimensional over the slotted and perforated 
models but that the leading edge peak pressure coefficient tended to vary along the 
span. In addition, transition tended to move forward with increases in Reynolds 
number in a somewhat non two-dimensional fashion. In general, transition tended to 
be more rearward toward the ceiling where leading-edge pressure peaks were lowest, 
and more forward toward the floor where they were higher. 


The measured and theoretical suction coefficient 


( pw) 
(C Q = (pu) 



00 


distributions 


required to maintain full-chord laminar flow over both the slotted and perforated 
upper surfaces at M design =0.82 and R c » 10 x 10® are shown in figure 19. The 

design suction requirements were based on shock-free flow and minimum level required 
to maintain laminar flow. The measured suction required for full-chord laminar flow 
on the slotted model was higher than that predicted, in some cases by as much as 50 
percent, and generally represents the maximum suction capability of the system. The 
measured suction requirements on the perforated model were much higher than on the 
slotted, particularly over the seven suction ducts on the forward panel 
(x/c < 0.26). The relative magnitude of the measured suction requirements over the 
next five suction ducts on the center panel (0.26 < x/c < 0.59) varies somewhat, but 
over the remaining nine suction ducts on the aft panel (x/c > 0.59) , the available 
suction capability of the perforated panel was much lower than that theoretically 
required or that available on the slotted model. The measured suction distribution 
over the aft panel of the perforated model represents the maximum suction capability 
of the panel as designed and installed. The measured suction distribution over the 
forward and center panels represents the maximum suction available after extensive 
hardware modification has been made to increase their suction capability. These 
modifications include enlarging all metering holes in the flutes (fig. 13) in the 
two panels, enlarging the laminar nozzle extensions in the two panels, and adding a 
second suction nozzle to each laminar duct in the forward panel. 


There was, of course, an infinite combination of individual duct suction levels 
and overall suction distribution possibilities. The distributions shown (fig. 19) 
consist, in general, of the maximum suction capability of each suction duct combined 
with the maximum suction capacity of the compressor system used to provide suction 
to the model. Small local variations may be permissible within these overall 
distributions without an adverse effect on the extent of laminar flow within the 
resolution permitted by the chordwise spacing of thin films. (See fig. 2). 

Reductions in the overall suction level of the distributions, by varying compressor 
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controls in amounts large enough for the sum of the suction drag over the entire 
upper surface to be measurably different, generally resulted in either a detrimental 
effect on the laminar-flow pattern (transition behavior) or an increase in wake 
drag. The higher than theoretically predicted suction requirements on the upper 
surface of both models are attributed to uneven velocities/pressures on the upper 
surface discussed above and the wind tunnel environment. 

The higher suction requirements on the forward region of the perforated model 
compared to the slotted model was attributed to several factors. These factors 
include the less rigid construction of the perforated panels and the resulting 
greater, and more uneven, deformation under load; the waviness of the perforated 
surface, some of which is associated with the corrugated substructure, is also 
greater than on the slotted panels. 

Since the lower slotted suction surface was the same for either the slotted or 
perforated model upper surfaces, an example comparison between the measured and 
predicted suction is shown in figure 20 for M = .82 and R Q = 1 0 x 10 6 . The measured 
higher than predicted suction requirements are attributed to the suction control 
required for the minimization of centrifugal Gortler type boundary- layer 
instabilities and interactions with cross flow in the concave regions of the lower 
surface (fig. 8). 

The contributions to the total section drag coefficient for the slotted LFC 
airfoil at M =0.82 are shown in figure 21 over a chord Reynolds number range of 10 
to 20 million. The division of suction drag contributions between the upper or 
lower surfaces may be made since the suction drag is computed duct-by-duct and 
integrated over each surface independently. The wake drag is separated into upper 
and lower surface components on the basis of the assumption that the wake can be 
divided between the upper- and lower-surface at the point on the wake rake where the 
stagnation pressure loss is the greatest. The data (fig. 21) indicate that the 
larger contribution to the total drag is from the lower surface. With full-chord 
laminar flow over the upper and lower surfaces for R Q < 12 x 10 s , the contribution 

to the total drag was about 1/3 due to wake and 2/3 due to suction drag. The sharp 

rise in wake drag on the lower surface between 14 and 15 million Reynolds number is 
associated with rapid forward movement of transition (fig. 17) and separation of the 

boundary layer in the lower aft cusp. The contribution of the suction drag is about 

40 percent for the upper and 60 percent for the unconventional lower surface 
(fig. 21). 

Wind tunnel tests have been conducted earlier on several swept LFC airfoils at 
low speeds (refs. 10, 51) with extensive laminar flow and low drag. These airfoils 
■had suction applied only on the upper surface. Figure 22 shows the upper-surface 
pressure distributions and minimum total drag measured at the design lift condition 
(C. * .3) on two earlier low-speed LFC airfoils of different sweep as compared to 
the present supercr itical LFC airfoil with C^ = .55. The low-speed LFC designs, 
based on standard NACA airfoil profiles, have favorable pressure gradients over the 
first 50 percent or more of the chord, and much less severe adverse pressure 
gradients aft, than the supercritical design. The higher design C L of the SCLFC(I)- 
051 3F airfoil does not permit a long run of favorable pressure graident . Figure 22 
shows only the upper-surface measured minimum profile drag coefficients 
corresponding to the upper-surface pressure distributions at R Q = 10 x 10 6 . As 
might be expected, the supercritical design has a somewhat larger suction drag 
penalty than the N0RAIR model, which has the greatest extent of favorable pressure 
gradient (fig. 22), and only a slightly larger suction drag penalty than the 
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University of Michigan 5- by 7-Foot Tunnel model, while the wake drag contribution is 
about the same in all cases. An increase in wake drag for M = 0.70 is observed on 
the slotted LFC airfoil (fig. 22) and attributed to the formation of a weak shock 
wave near the leading-edge as the supersonic bubble begins to develop. As the 
bubble develops near the leading edge (0.78 < M o < 0.80), full-chord laminar flow 
still exists, but periodic turbulent bursts occur over the upper surface causing an 
increase in wake drag. It may be concluded that the basic phenomenon of applying 
suction laminarization over an extensive supercritical zone with full-chord laminar 
flow at R q = 1 0 x 10 s is feasible as demonstrated on a swept LFC airfoil at high 
lift conditions. It is clear, however, that the "nonideal" transonic wind tunnel 
test environment (refs. 42, 44) and surface waviness (refs. 40, 41) do cause a loss 
in the amount of laminar flow achieved. Consequently, the data obtained are 
conservative; i.e., one should do better in flight. 

Figure 23 shows the in-flight upper surface wing and glove pressure 
distributions for the X-21 and F-94 aircraft compared with that measured on the 
SCLFCO )-051 3F airfoil. Full-chord laminar flow was achieved for each experiment 
and conditions shown. The present wind tunnel model has not only a larger local 
Mach number- but extent of sonic zone than the flight (refs. 52, 53) results, and also 
has a much higher lift coefficient. The X-21 wing and F-94 glove have favorable 
pressure gradients in the nose region and much less severe adverse gradients in the 
aft region than the supercri tical airfoil. 

Figure 24 shows a comparison of the drag measured on the upper surface of these 
three configurations. With full-chord laminar flow on the present model for 
R c £ 12 x 10 s , the total drag level was about the same as that measured in flight on 
the X-21 or F-94. Detailed suction drag measurements were not determined in flignt 
on the X-21, but an estimate of the total drag was made from reference 55. The 
higher total drag measured for the present LFC airfoil for R c £ 12 x 10 s was 
discussed earlier and is attributed mainly to the surface deformation, waviness, and 
the wind tunnel environment. However, the indicated drag level for the upper 
surface of the SCLFC (1 ) —051 3F airfoil, obtained with full-chord laminar flow, is 
about one-third that of the upper surface of a conventional turbulent airfoil. 

High-Speed HLFC Airfoils 


Another and more practical concept for obtaining long runs of laminar flow at 

high Mach number is a hybrid configuration which would combine suction over forward 

regions of the upper surface with natural laminar-flow (NLF) concepts over rearward 

regions (refs. 2, 3). An attempt was made to simulate such hybrid laminar-flow 

control (HLFC) conditions during the 8-Ft TPT LFC experiment "by turning off suction 

over regions of the upper surface. On the lower surface, suction was maintained for 

0 £ x/c £ 0.25 to control the flow between the lower surface and tunnel wall. This 

simulation of HLFC involved first establishing laminar flow to the most rearward 

location, then progress! vely turning off upper surface suction starting with the 

most rearward suction duct and measuring the extent of laminar flow that was 

maintained downstream of the turn-off points. The variation of transition location 

on the upper surface with chordwise extent of suction is shown in figure 25 for 

= 0.82 and two chord Reynolds numbers. The results indicate that laminar flow 

could be maintained well beyond termination of suction. With no suction, laminar 

flow was present back to about 15“ and 25- percent chord, depending on R Q . At 

R = 10 x 1 0 6 , the measured total drag was C^ - 0.0012 with 0 £ (x/c) .. £ 0.96 

u suction 

and increased to C^ ^ 0.0025 with no suction. The C c j w contribution increases and 
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becomes a larger fraction of C dtotal as suction is reduced to zero. Clearly, this 

simulated HLFC approach must be considered somewhat "nonideal istic" and the results 
conservative since no effort was made to seal or smooth the unsucked regions of the 
airfoil. Furthermore, the pressure distributions (figs. 14 and 15) are not like 
that which would be designed for an HLFC configuration. 

Almost any form of LFC could be considered a hybrid conf iguration, since, from 
a practical standpoint, one would probably not apply suction control over the full 
chord on both surfaces. Furthermore, considerable research is required to establish 
the effectiveness of suction control through shock boundary-layer interaction 
regions that occur at off-design conditions. Thus, the attractiveness of a form of 
HLFC that requires limited suction control in the leading-edge region for cross-flow 
instability followed by a favorable pressure gradient to maintain attached laminar 
flow (as far rearward as possible) is substantial. This concept would provide the 
maximum drag reduction with minimum suction requirements and complexity at design 
and off-design transonic speeds. 

NASA Langley has designed, fabricated, and is currently installing an HLFC 
airfoil in the 8-Ft TPT for testing. The HLFC concept was originated by P. J. 
Bobbitt and the design was carried out by S. H. Goradia and J. C. Ferris and is 
reported in reference 56. It should be noted that constraints were placed on the 
HLFC design that required utilization of the present LFC airfoil and liner 
configurations in order to minimize costs and schedules for fabrication and 
testing. Only the upper mid and aft model surface regions (0.25 £ x/c £ 1.0) were 
altered in order to achieve a somewhat favorable pressure gradient downstream of 
suction turn off. The upper forward perforated and lower three-slotted surface 
panels were retained (figs. 9, 12). A new flap was designed and required to match 
the new upper aft surface camber change. 

Figure 26 shows the HLFC airfoil design upper surface pressure and suction 
distribution for M = 0.81, C, = 0.47, R c - 15 x Ip 6 and sweep angle of 23°. 

Suction control is°°applied from 0.025 £ x/c S 0.26 through the original perforated 
panel for the HLFC concept. The HLFC type pressure distribution yields somewhat 
less lift than that for LFC in an effort to avoid shocks and aft separation. 
Predicted drag coefficients (fig. 27) for the upper surface of the HLFC airfoil are 
expected to be on the order of two times the drag levels for the full-chord laminarr 
flow shockless LFC airfoil for E c < 12 x 10 6 . 

DRAG REDUCTION SUMMARY 

Figure 28 is an attempt to summarize the drag reduction achieved for the 
HSNLF( 1 ) —021 3F and SCLFC( 1 )— 051 3F airfoils designed for transonic operation, 
compared with several other recent NASA laminar-flow airfoils developed for low 
speeds shown as a hatched region in the figure. The results shown are for the 
minimum total drag measured which includes both wake and suction drag. The symbols 
are for drag values measured on these airfoils with extensive laminar flow over the 
upper and lower surfaces compared to fully turbulent flow conditions on airfoils as 
also indicated by hatched band. The values shown correspond to thejninimum drag 
levels measured with a and 6 f = 0° but for different R ? 's. The data (fig. 28) 
indicate that an average drag value of about 0.0085 exists when the flow is 
turbulent and M cos A < 0.75 and about 0.0035 with laminar flow achieved by NLF or 
LFC. This represents a total drag reduction of about 60 percent that has been 
achieved and extended to transonic speeds with the new NASA airfoils. Apparently, 
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the presence of a shock on the unswept HSNLF(1 )-021 3F airfoil with laminar flow 
causes wave drag that influences drag rise for M w > 0.7. 


STABILITY METHODS ANALYSIS 

Empirical relations, derived from earlier low-speed and Reynolds number 
experiments, are available and may be used for predicting transition (refs. 5 ' ’ 

58 ). More sophisticated approaches have been developed to determine the transit 
location on swept and unswept wings for incompressible (refs. 20, 59) and 
compressible (refs. 60, 61) flow. 

The more sophisticated methods used for design and analysis of the airfoils 
herein have been described and presented earlier (refs. 20. 59-61). These methods 
are based on the temporal stability theory for laminar boundary layers where the 
local amplification rates are obtained as solutions of the governing Orr-Sommerfeld 
equation as functions of the frequency, wavelength, and orientation angle o e 
disturbance. If it is assumed that a disturbance has an initial ampiitude A 0 at 
some location X n in the boundary layer, then the amplification ratio A/A Q a y 
downstream point X can be calculated. This amplification ratio can be expressed as 
an exponential function e N where N is referred to as the amplification exponent or 
N-factor. N-factors may be related to the growth of disturbance amplitudes and 
used .based on experience, to predict transition conditions. These analysis metnods 
only apply in linear attached flow regions and calculations must be carried out to 
the extent that the amplification exponent N becomes a maximum .or all 
frequencies. When using the N-factor method for design, one is left with the 
arbitrary or empirical choice of a limiting N-value for transition. When used fo 
analysis, limiting N-values are determined by the experimental transition locati 

Calibration of Stability Codes 

A number of computer codes are available to calculate N-factors for Tollmien- 
Schlichting (TS) waves and cross-flow (CF) vortices. Some examples are he 
incompressible SALLY (ref. 59) and MARIA (ref. 20) codes and the compressible MACK 
(ref. 61) and C0SAL (ref. 60) codes. The MARIA code is the most str ^g h ^°rwar 
and economical to apply but is restricted to CF calculations only. The SALL 
may be used for CF or TS analysis. 

Experimental transition data have been recently obtained and analyzed on NLF 
and LFC airfoils (refs. 3, 62), on bodies of revolution (ref. 63), in wind tunnels 
and on wings or gloves In flight (refs. 3, 62, 64-70). The measured pressure 
profiles and suction distributions from a number of these experiments were used 
input to the above mentioned codes to calculate local disturbance growth rates and 
integrated amplification ratios (N-factors). N-factors calculated at the 
experimently measured transition locations for both TS and CF disturbances were 
compared to determine which instability was the most amplified. 

Figures 29 and 30 are an attempt to summarize the calculated lowest N-factors 
at measured transition for TS and CF, respectively. It should be noted that there 
exists higher ref erencable N-factor values than shown in the figures; only the lower 
values were chosen to establish a low N-factor boundary and Mach number trend. In 
figure 29, calculated N-factors for situations where the TS instabilities dominate 
(CF negligible) is given; figure 30 is for situations where CF dominates (TS 
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negligible). The N-factors shown for the NASA NLF(1)-0414 wind tunnel and Cessna 
210 flight tests and the NASA SCLFC ( 1 )-051 3 wind tunnel tests with slotted suction 
surface were calculated using the envelope method of solution in the SALLY and COSAL 
codes, with identical input data for each point (ref. 62). All other results on 
figures 29 and 30 (refs. 64-70) were calculated using either the SALLY, COSAL, or 
MACK codes. The incompressible TS N-factors (open symbols) approach a constant 
lower limit of about N = 9 — while those N-factor values obtained with the 
compressible codes (solid symbols) decrease over the indicated speed range (solid 
curve) and are about 50 percent lower at cos A “ 0.8. These trends (fig. 29) are 
intended to establish a theory calibration^of the allowable lower limit for 
disturbance amplification. The limited number of incompressible CF N-factors (fig. 
30 ) have a minimum value of 9 while those from the compressible code decrease to 
values of 4 or 5. The results strongly indicate that, on the average, a constant 
(N = 9) but conservative N-factor value can be used in applying the incompressible 
codes. However, knowledge of both the conservative trend and level of the N-factor 
(figs. 29, 30) is required when applying the compressible codes. The consistency of 
this correlation validates and extends the use of both incompressible and 
compressible linear stability theory for conservative design to high speeds. The 
fact that the lower limit of the N-factor trend from the compressible codes decrease 
with Cos A implies that "compressibility effects" are not conservative in 
stability analysis; i.e., transition takes place with less amplification than for 
incompressible flow. 

All of the N-factors calculated and presented in figure 29 for the LFC airfoil 
were based on the beginning of transition location (ref. 62). However, if the 
transition point is selected further downstream of this location, the incompressible 
or compressible boundary— layer codes generally predicted laminar separation leading 
to a very rapid increase in the growth of the N— factor over a very narrow frequency 
range. This also caused arbitrary selection of amplification factors that could be 
on the order of twice or more those values shown in figure 29 for the SCLFCO )— 051 3F 
airfoil. Calculated values of N TS along the slotted LFC airfoil chord were found to 
vary depending upon the chordwise’ extent of applied suction and accuracy with which 
transition location could be determined. For example, with laminar flow over less 
than 30-percent chord on the upper suface, the calculated N T q values were from 8.5 
to 17.5 (ref. 62). However, with laminar flow beyond 30-per’cent chord, the N TO 
values were found to decrease dramatically below these values to about 1 or 2, 

Thus, it becomes very important to account for this chordwise variation in N-factor 
for a given design case. The indicated trends (figs. 29~30) could be attributed to 
the wind tunnel disturbance environment; however, flight N-factors are shown that 
agree with the wind tunnel data. 

Figure 31 is an attempt to evaluate the influence of transonic wind tunnel 
turbulence level on the calculated N-factors at transition on models. Shown for 
comparison are calculated values of N TS on the slotted SCLFCO )-051 3F airfoil in the 
Langley 8-Ft TPT (ref. 62) and two bodies of revolution in the NASA Ames 12-Ft PWT 
(ref. 63). Also shown is the empirical relation for the variation of N T q with 
turbulence level generated by MACK (ref. 71) for u/u > 0.1 percent. With the 
exception of a single data point (N TS = 10.5), all of the calculated values tend to 
fall within a horizontal band between N TS values of about 5 to 8 indicating little 
or no influence of tunnel turbulence level as suggested by references 63 and 71 . In 
fact, the LFC airfoil results shown for u/u =0.2 percent (M = 0.7) are for 
essentially full-chord laminar flow on both surfaces up to R°° < 22 x 10 s . For 

u/u = 0.05 percent (M^ = 0.82) full-chord laminar flow existed'for R = 10 7 and 
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N T q * 4.5 and with 0.2 < (x/c) fcr < 0.3 the N TS = 7.5 for R c = 2 x 10 7 . Therefore, 
the present limited results (figs. 29, 30 ) suggest the need for further calculations 
and analysis in order to complete the calibration of compressible boundary- layer 
stability theories. Note that the u/u values just quoted are for single hot-wire 
measurements. 


CONCLUDING REMARKS 

Significant amounts of laminar flow have been achieved through passive ( NLF ) or 
active (LFC) methods at high speeds. This has been demonstrated on recently 
developed low-drag airfoils which use boundary-layer suction control or favorable 
pressure gradient. Wind tunnel tests of these concepts, from subsonic to transonic 
speeds, have shown significant improvements in lift-to-drag ratio over previous 
airfoils. 

Natural laminar flow was maintained over 50-percent chord and shock-free 
conditions on an unswept airfoil designed for high speed and moderate lift. A 
reduction in the drag of about 40 to 50 percent over turbulent designs was 
realized. Suction laminarization , through either a slotted or perforated surface, 
over a large supercritical zone has been shown to be feasible to high-chord Reynolds 
numbers even under "nonideal" surface and test environment conditions on a swept LFC 
airfoil at high-lift coefficients. With essentially shock-free flow at 
M = 0.82, full-chord laminar flow on the upper and lower surface was achieved for 
r” < 12 x 1 0 6 . As Reynolds number increased above R c > 12 x 1 0 6 , transition moves 
gradually forward on the upper surface but laminar flow was maintained over most of 
the supersonic zone on the slotted and perforated surfaces for M* =0.82 . The 
extent of suction laminarization was found to be slightly less for the perforated 
surface than for the slotted surface and attributed to a less rigid model and a more 
wavy surface. A reduction in the drag level of about 80 percent for the upper 
surface only and about 60 percent for both upper and lower surfaces was achieved 
with full-chord LFC. 

Calculations using measured wind tunnel and flight data show differences in the 
trends and the lower bounds of the amplification ratios as calculated by 
incompressible and compressible linear boundary-layer stability theory. Lower 
N-factor levels have been obtained from the compressible stability theory than the 
incompressible; consequently, compressible calculations are not conservative, as 
previously thought. 
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Comparison of drag polars on airfoils at low and high speeds with and 
without flap deflection. 
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An example of flush-mounted surface hot-film gage locations, RMS output 
signals, and laminar pattern. 


129 




PHASE PLOTS 


(x/c) i to (x/c) 2 



(a) Upstream of the shock, (b & c) Across shock-induced separation and reattachment, (d) Downstream of reattachment 
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x/c 

Figure 3* Example of multielement dynamic shear stress sensor detection of shock- 
induced separation by phase reversal. 



Figure 4. Design theoretical pressure distribution and geometry for HSNLF( 1 )-021 3 
airfoil. = .7, = .26, R c = 1 1 x 10^, a « 0°. 


130 



Figure 5. Comparison of experimental and theoretical pressure distributions for 
HSNLFO ) -021 3 airfoil. 


M = 0.77, R c = 4.3 X 1 0 6 , C l= 0.26 



Figure 6. Comparison of experimental and theoretical pressure distributions for 
HSNLFO ) —021 3 airfoil at off-design conditions. 
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Figure 7. 
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TRANSITION 
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Effect of Mach number on drag coefficient for the HSNLF( 1 )-021 3 airfoil 
in two different wind tunnels at zero angle of attack. 



Large chord = 7.07 ft 
Upper and lower surface suction 


FREE-STREAM DESIGN CONDITIONS 


M c0 = 0.82 
A = 23° 

C = 7.07 ft 
R c = 20 X 1 0 6 
C| =0.47 


. Schematic of swept SCLFC(1)-051 
conditions. 


DESIGN CONDITIONS NORMAL 
TO LEADING EDGE 


1^= 0.755 
cjg= 6.508 ft 
R c N= 16.9 x 10® 
(t/c) N = 13.0 percent 
C| N = 0.55 


3F airfoil, sonic zones, and design 
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UPPER PANELS 



LOWER PANELS 

Figure 9. Schematic of SCLFCO ) -051 3F airfoil wingbox and suction panel assembly 



Figure 10. Schematic section of slotted suction surface panel concept. 









(a) View looking downstream at model leading edge and upper surface. 

Figure 11. Photograph of swept supercritical LFC airfoil model and contoured wall 
liner in the Langley 8-Foot Transonic Pressure Tunnel. 
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(b) View looking upstream at model trailing edge and liner step. 
Figure 11. Concluded. 
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Figure 1*J. Comparison of experimental and theoretical pressure distributions for 
SCLFC ( 1 ) -051 3F airfoil with slotted and perforated surface. 

R c > 10 x 10 6 . 
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Figure 15. Comparison of experimental and theoretical pressure distributions for 
SCLFC( 1 )-051 3F airfoil with slotted and perforated surfaces. 

R c = 20 x 10 6 . 
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Figure 16. Measured transition location on upper surface of slotted SCLFC(1 ) -05 1 3F 

airfoil. M = 0.82. 
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Figure 17. Measured transition location on lower surface of slotted SCLFC (1 )-051 3F 
airfoil. M = 0.82. 



Figure 18. Measured transition location on upper surface of perforated SCLFC(I)- 

0513F airfoil. M = 0.82. 
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Figure 


Figure 


Final Configurations, 10 Million Reynolds Number, M®- 0.82 



19. Comparison of experimental and theoretical suction distributions on 

upper surface of slotted and perforated SCLFCO )-051 3F airfoil models. 

M = 0.82, FU = 10 x 10 6 . 



20. Comparison of experimental and theoretical suction distributions on the 
lower surface of slotted SCLFC( 1 ) -051 3 airfoil model. 

M b = 0.82, R c = 10 x 10 6 . 
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Figure 21. Measured variation of drag with chord Reynolds number on SCLFC(1 ) -05 1 3F 

airfoil with slotted suction surface. M = 0.82. 
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Figure 22. Comparison of measured pressure distributions and total drag for upper 
surface only on LFC wind tunnel models over speed range. 
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Figure 23. 
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Figure 24. 



Comparison of measured local Mach number distributions on upper surface 
of wind tunnel and flight LFC wings at transonic speeds. 
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Comparison of drag coefficients for upper surface of swept and unswept 
transonic suction wings. 
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Figure 25. Variation of measured transition for simulated hybrid LFC on swept 
SCLFCO ) -051 3F airfoil with slotted upper surface. 



Figure 26. Design theoretical pressure and suction distributions on upper surface 
of hybrid LFC airfoil. M = 0.811, C, = 0.H7, R - 15 x 10 6 , A = 23°. 
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Figure 28. 


A = 23 , C |_= .5, Mco= .8 



Measured and projected effect of suction extent on drag for LFC and HLFC 
airfoil upper surface. 
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Comparison of measured minimum drag on NASA advanced low-drag airfoils, 
with and without suction, at transonic speeds with previous results 
(hatched areas). 
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Figure 30. Minimum allowable crossflow disturbance amplification N-factors at 
measured transition. 
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Figure 31. Variation of compressible N-f actors at transition with turbulence level. 
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INTRODUCTION 

The X-29A Advanced Technology Demonstrator, figure 1, has been designed and developed to 
demonstrate the advantages of forward wing sweep along with other advanced technologies. Grum- 
man Aerospace Corporation designed and fabricated the flight vehicle under a contract sponsored by 
the Defense Advanced Research Projects Agency and managed by the Air Force. A number of ad- 
vanced technologies have been incorporated into the design. These include, in addition to the obvious 
forward-swept wing, a closely coupled canard, thin supercritical wing sections, discrete variable cam- 
ber, aeroelastic tailoring, three-surface control, relaxed static stability, and a triple-channel digital flight 
control system. 

An extensive flight-test program has been defined and is underway at NASA Ames/Dryden Flight 
Research Facility to acquire high-quality data which will be used to evaluate and correlate the original 
design analyses and ground-test results with flight results. One of the objectives of the flight program 
is to isolate the benefits of specific technologies. However, when one considers the complex interactions 
occurring in the flight environment and the difficulty of extracting specific parametric variations from 
the available data, the task takes on herculean proportions. 

A cooperative program has been defined between the Fluid and Flight Mechanics Branch at NASA 
Ames/Dryden and the High-Reynolds-Number Aerodynamics Branch at NASA Langley to support and 
complement specific flight test objectives. In addition to the flight test elements, the program consists of 
wind tunnel experiments and computational fluid dynamics to enhance the understanding of interacting 
technologies and to isolate individual benefits. The following sections will present brief overviews of the 
flight and ground experimental elements and a more complete discussion of the computational support. 

This paper will mainly address issues and questions associated with the forward swept wmg and 
closely coupled canard. The primary focus will be on research questions which must be addressed to 
obtain high quality ground and flight test data. These data will be used in conjunction with computa- 
tional predictions to complement the analyses required to comprehensively understand the interacting 
technologies. 


BENEFITS OF THE X-29A WING AND CANARD 

Before the three phases of the program are discussed, it is appropriate to include a brief discussion 
of the benefits associated with the forward-swept wing and the canard. For transonic flight, wing 
sweep is used to delay the onset of drag rise. Uhad, et al. (ref. l) demonstrated experimentally 
that a forward-swept wing offers potentially lower wing profile drag than an equivalent aft-swept wmg 
for a transonic maneuvering design point. This yields higher sustained lift coefficients at transonic 
maneuvering conditions. If one compares aft- and forward-swept tapered wings, the local sweep angles 
for the forward-swept wing increase as one progresses from the wing leading edge to the trailing edge, 
while the opposite is true for an aft-swept wing. This results in lower wave drag losses for similar flow 
conditions and shock locations. Moore <md Frei (ref. 2) similarly showed that the trade-off between 
forward- and aft-swept wings could favor the forward-swept wing. Their comparison configurations 
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held the wing area, aspect ratio, taper ratio, and shock sweep constant between the forward- and aft- 
swept wings. The forward swept wing had significant reductions in the wing bending moment at the 
pivot resulting from an inboard shift in the wing center of pressure. If one maintains the same shock 
sweep for the comparisons, the forward swept wing has a lower leading-edge sweep. This allows specific 
supercritical airfoils to be used which have demonstrated reductions in profile drag with decreasing 
leading-edge sweep. If one imposes a constraint of constant bending moment at the pivot and removes 
the constraint of constant aspect ratio, the forward-swept wing yields a higher aspect ratio in comparison 
to the aft-swept wing. This, of course, is beneficial in reducing induced drag. 

Improved handling characteristics for forward swept wings at higher angles of attack were demon- 
strated in an experimental wind tunnel program (ref. 3). In general, flow separation on a forward-swept 
wing begins inboard and progresses toward the wing tip. In contrast, the tendency for an aft-swept wing 
is for the wing tip to separate first with the separation progressing inboard. The X-29A wing shows a 
mild break in the lift curve slope at low speeds and moderate angles of attack (10-12 degrees) resulting 
from inboard separation. However, the lift and pitching moment data are well behaved to high angles of 
attack with no pitch-up tendencies. Data also showed the X-29 full-span ailerons can maintain adequate 
roll control to very high angles of attack. 

In addition to the benefits derived form the forward-swept wing, the close-coupled canard yields 
distinct advantages for the configuration. The canard and wing were designed in conjunction to yield 
an approximately elliptical span loading distribution on the wing at the maneuver design point. The 
canard is in the strong upwash field of the wing which results in an effective moment arm at subsonic 
speeds which is twice that of the geometric moment arm. Another advantage of the close-coupled canard 
results from the effect of the canard downwash on wing-root separation. At higher angles of attack, the 
lift generated by the canard to trim the configuration results in a downwash over the inboard part of the 
wing. This allows the inboard flow to stay attached to higher angles of attack than would be possible 
without the canard influence. The effect of the canard on the wing loading is presented in figure 2. These 
data are from the l/8-scale model test at NASA Ames (ref. 4). The strong influence of the canard on 
wing span loading is evident as the canard deflection angle is changed at a given angle of attack. The 
loading is markedly shifted outboard as the canard deflection is increased. 


FLIGHT TEST PROGRAM 

As a subset of the flight-test objectives for the complete X-29A flight program, specific objectives 
have been identified which are associated with the forward-swept wing and close-coupled canard. To meet 
these objectives, extensive data must be gathered and analyzed over a wide range of flight conditions. 
The aircraft has been extensively instrumented to gather these data, and software has been developed 
to assist in data manipulation and analysis. Figure 3 indicates some of the instrumentation which will 
be used. Static pressure orifices are arrange in five rows on the wing (r; = 0.2, 0.31, 0.49, 0.70, and 
0.91) and two rows on the canard (rj = 0.28 and 0.65). These pressure data are extremely informative 
relative to inferences concerning wing-canard interactions and stall onset and progression. In addition, 
in-flight flow visualization will also be used to qualitatively identify flow characteristics through the use 
of flow cones. These are similar to tufts having a small plastic cone attached to the free end to enhance 
visibility and increase durability. 

One of the potential problems associated with a forward-swept wing is aeroelastic divergence. In 
order to build a wing using conventional materials and manufacturing techniques strong enough to 
counteract this divergence, a significant weight penalty is incurred. Krone (ref. 5) presented the concept 
of advanced aeroelastic ally tailored composites, which are used in the X-29A vehicle, to overcome this 
problem. The concept uses a buildup of specifically oriented plies of composite laminate to form the wing 
skins. The orientation and characteristics of the plies can be designed to have the strength necessary to 
overcome divergence. An additional advantage is that the technique can be used to tailor the aeroelastic 
twist distribution at maneuvering conditions. 

A wing-deflection measurement system is installed on the right wing of the flight-test aircraft. Two 
fuselage-mounted receivers with different focal lengths are used to monitor 12 targets mounted on the 


wing surface, figure 3. These data can be reduced to determine discrete deflections at each of the targets 
or integrated to determine the average twist angle and section translation for the various span locations. 
These data will have several uses including validation of structural design methodology. The data 
are critical in order to provide valid, reliable, transonic computational predictions of the configuration 
flow field for the flight vehicle. The supercritical wing is an efficient but sensitive performer at transonic 
speeds. Small changes in flow conditions or geometry (such as aeroelastic deflections) can have significant 
effects on the flow over the wing surface including shock location, shock strength, flow separation, and 
buffet onset. 

A number of challenges have been encountered as the flight-test program has progressed. Some of 
these difficulties have centered around data analysis. The aircraft’s negative static margin and the use 
of canard, wing flaps, and strake flap simultaneously, figure 1, for pitch control make it exceedingly 
difficult to separate the effect of specific variables. In normal flight the deflections of these control 
surfaces are controlled by the Automatic Camber Control (ACC) system. Hence, the ability to isolate 
specific stability derivatives is quite limited. However, the Manual Camber Control (MCC) system was 
implemented to limit the number of parameters allowed to vary for control purposes. This allows the 
flap to be fixed at a discrete setting during a given portion of the flight, thereby reducing the number 
of control surfaces whose deflections are varying. An example of different inferences derived from data 
for similar flight maneuvers using the two different modes (ACC and MCC) is presented by Waggoner, 
Jennett, and Bates in reference 6. 


NTF MODEL AND WIND TUNNEL TEST 

A 1/16-scale model of the X-29A, figure 4, is being designed and fabricated for testing in the National 
Transonic Facility (NTF) at NASA Langley. The NTF was conceived to provide high-Reynolds-number 
test capability for aerodynamic research and development testing of commercial and military aircraft 
configurations. The NTF is a closed-circuit, single-return, fan-driven, wind tunnel capable of continuous 
operation over a Mach number range of 0.2 to 1.2. It is capable of operating at very low temperatures 
(-320°F) by injecting nitrogen at cryogenic temperatures into the stream. By operating at elevated 
pressures and cryogenic temperatures, a maximum chord Reynolds number of 120 million at Mach 
number of 1.0 for a chord of 0.82 feet is achievable (ref. 7). 

The X-29A test program is intended to provide significant data which will complement the flight-test 
program. In addition to providing basic data for flight-to-tunnel and code-to-tunnel correlation studies, 
Reynolds number effects on a number of sensitive flow areas will be investigated. Over the complete NTF 
testing schedule for the X-29A, it is anticipated that Reynolds number effects will be determined on wing 
surface pressures and wing-canard interactions; supercritical design; shock/boundary layer interactions; 
performance, stability and control; and high- angle-of- attack characteristics. 


1/16 SCALE MODEL 

The NTF model has been fabricated by Grumman Aerospace Corporation under contract with 
NASA Langley. The model can be used to simulate aerodynamic control surface deflections through the 
use of variable incidence canards and interchangeable sets of wing flaperon, rudder and strake flap parts. 
The design of the model allows for the simulation of four wing flap settings, seven strake flap settings 
and four rudder settings. In addition, six positions are available for the variable incidence canard. The 
external lines of the X-29A are modeled in significant detail. Flow through nacelles will afford the first 
simulated inlets for cryogenic operation in the NTF. The flow through ducts will simulate the mass flow 
rate for a condition near the cruise design point. Although the provisions to simulate various mass flow 
rates are not anticipated, the calibration of the inlets for various Reynolds numbers presents a significant 
challenge. 

Data acquisition planned for the test program includes aerodynamic force and moment data, wing 
pressure distributions and wing deformation characteristics. An extensive array of instrumentation will 
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support this task. Static pressure orifices are located at span locations corresponding to wing stations 50 
and 114. The orifice locations were chosen to correspond closely to the orifice locations on the number- 
one flight-test aircraft. A 48-port electronic sensing pressure (ESP) module will be mounted in the model 
nose, as well as an accelerometer-type angle of attack indicator. 

Two, six-component, 1.75-inch diameter, internal strain-gage balances have been constructed. Each 
is intended to be used over a different load range. The normal force design loads of the two balances 
are 1500 and 2500 pounds. It was desirable to have static pressure measurements on the canards. 
However, because of the size of the canard (1/8-inch thick at mid-span) and small diameter of the shank 
attachment to the fuselage, it was not feasible. In lieu of this, a three-component strain-gage balance 
is built into the left canard to measure shear, bending and torsion. A wing-root bending- moment gage 
allows for real-time evaluation of model dynamics. 

Deflection of the wings during the wind tunnel test will be measured by use of a video model 
deformation system. These measurements will be used to determine the magnitude of differences in 
wing shape between the model and flight vehicle under similar flow conditions. 

The fabrication of the model was somewhat unique. The left wing and fuselage half and right wing 
and fuselage half were each machined as one piece out of a high-nickel-content steel suitable for cryogenic 
models (Vascomax 200). The two halves were then welded together yielding a one-piece wing fuselage 
model. The interchangeable flaps also presented a serious design problem. The flap elements are quite 
thin and were designed to attach to the wing through a tongue-in- groove arrangement. This design was 
necessary to ensure compatibility with the severe thermal test environment. The thin tongues on the 
flaps and the precision machine work for the grooves presented a challenge for the model craftsmen. 

NTF MODEL TEST PROGRAM 

The model design point simulates flight Reynolds number at M = 0.9 and 30,000 ft. with the 
NTF operating at minimum tunnel dynamic pressure. Simulated 8-g flight at these conditions yields a 
Cn • q limit of 2,700 psf for the 1/16-scale model. A comparison of the tunnel-test-to-flight envelope 
is presented in figure 5. Superimposed on the envelope are the sting divergence limit boundary for the 
low-angle of attack support system and the envelope of Reynolds numbers attainable with this model in 
the NASA- Ames 11-Foot Tunnel. Because of the versatility of the NTF, note that even with the sting 
divergence limit imposed, much of the airplane flight envelope below M = 1.2 can be simulated with the 
1/16-scale model in the NTF. At Mach numbers less than 0.7, the entire angle-of-attack range can be 
tested at flight Reynolds numbers for the entire flight envelope. 

The model will be tested over a Mach number range from 0.3 to 1.2 at Reynolds numbers from 
2 to 35 million. Initially, the angle of attack will be limited to 20° ; however, with the high- alpha sting, 
angles of attack up to 60° can be achieved. 

COMPUTATIONAL SUPPORT OF THE X-29A EXPERIMENT 

A comprehensive computational program has been identified to support both flight and ground 
testing of the X-29A. The effort up to the present time has concentrated on the applicability of potential 
flow methods. The following sections will briefly describe the codes which have been used, describe the 
code calibration effort and present selected results. 

COMPUTATIONAL TECHNIQUES 

Three primary computational techniques have been identified to support the initial computational 
phase of the program. A three-dimensional, small-disturbance, transonic analysis code for wing-fuselage- 
canard combinations (CANTATA), a three-dimensional full-potential, transonic analysis code for wing- 
fuselage combinations (TAWFIVE), and an aerodynamic/structural analysis system (TAPS) will be 
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employed to exploit the advantages of each technique. Salient features of the computational techniques 
are discussed below. 

CANTATA Analysis Code - The Canard/Tail Transonic Analysis code (CANTATA, ref. 8) is 
characterized by a unique grid-embedding technique which provides excellent resolution for either a 
wing-fuselage-canard or wing-fuselage- tail configuration. The code solves for the flow field about the 
configuration of interest in the direct mode employing an ADI scheme. Design via numerical optimization 
is available in the code but is not anticipated to be used in the current application. Using finite difference 
approximations, a modified small-disturbance potential-flow equation is iteratively solved in a system of 
multiple embedded grids. The modifications to the classical small-disturbance equation are in the form 
of extra terms, which, when added to the equation, provide more accurate resolution of shock waves with 
large sweep angles and a better approximation of the critical velocity where the full- potential equation 
changes from elliptic to hyperbolic type. 

Viscous effects are approximated by coupling a modified Bradshaw boundary-layer computation to 
the finite-difference potential-flow solution. The modified method provides a technique to extend a two- 
dimensional bound ary- layer calculation to account for first order viscous effects (ref. 9). The viscous 
effects are incorporated in the solution by adding the boundary-layer displacement slopes to the wing 
surface slopes. This modifies the wing surface to an equivalent “fluid” wing shape which is then analyzed 
by the potential flow code. 

The basic concepts inherent to the code have been extensively validated by several researchers 
(refs. 10-12), These studies have, in general, been for aft-swept wings and for single lifting surface 
configurations. 

TAWFIVE Analysis Code - A computer code for the Transonic Analysis of a Wing and Fuselage 
with Interacted Viscous Effects, references 13 and 14, is also used in this study. The code utilizes the 
interaction of three-dimensional inviscid and viscous flow solvers to obtain transonic flow-field solutions 
about wing-fuselage combinations. The outer inviscid flow field is solved using a conservative, finite 
volume, full potential method based on FLO-30 by Caughey and Jameson (ref. 15), No modifications 
were made to the internal grid- generation algorithm in FLO-30, which is a body-fitted, sheared, parabolic 
coordinate system. 

A three-dimensional boundary layer for the wing is computed using a compressible integral method. 
The code has the capability of computing laminar or turbulent boundary layer with the methods of Stock 
(ref. 16) or Smith (ref. 17), respectively. An important feature of the code is Street’s treatment of the 
wake (ref. 13). The wake model used in FLO-30 was replaced with a model which satisfies flow tangency 
on the wake displacement body and the pressure jump condition resulting from wake curvature. These 
changes in the code can make significant differences in results obtained on various configurations. 

Transonic Aeroelastic Programs System - The Transonic Aeroelastic Program System (TAPS) de- 
veloped by Campbell (ref. 18) is a method which allows the effects of static aeroelastic wing deflections 
to be included in steady transonic aerodynamic calculations. The method interacts a three-dimensional 
transonic computer code with viscous effects and a linear finite element structural analysis code to cal- 
culate wing pressures and deflections. The nonlinear nature of the transonic flow makes it necessary 
to couple the aerodynamic and structures codes in an iterative manner. TAPS has been arranged in a 
modular fashion so that different aerodynamic or structural programs may be used with a minimum of 
coding changes required. 

For the current study the TAPS program has been utilized to help understand the correlations of 
computations and wind tunnel data with data obtained from the flight test. The flight vehicle has been 
analyzed in the TAPS program to predict the wing deflections under flight loads. These calculations 
will be described in detail in a subsequent section. 


CODE CALIBRATION EFFORT 

The first phase of the computational effort will involve calibration of the transonic aerodynamic 
codes for application to the X-29A configuration. Although these codes have been applied to numerous 
configurations, the X-29 A is such a unique configuration that a calibration phase is required. This utilizes 


both code-to-code and code-to-experiment comparisons. In general, for transonic analysis codes in 
widespread use for analysis and design, there exists an inverse relationship between the physics included 
in the governing flow equations and the complexity of geometry for which one is realistically capable of 
obtaining a solution. The two transonic potential flow codes discussed previously illustrate this point. 
The TAWFIVE code solves the full-potential flow equation coupled with a three-dimensional boundary- 
layer analysis for a wing-body configuration. The CANTATA code, on the other hand, solves a simplified 
flow equation (extended small disturbance) on a more complex (wing-fuselage-canard) geometry. 

The X-29A has a very complex geometry from a computational standpoint. The body, figure 6, 
was extremely difficult to model, particularly in the wing-fuselage juncture region. The discontinuous 
trailing edge resulting from the body strake geometry offered another challenge. The effects of the inlet 
on the flow over the canard and the dihedral effects for the wings were unknown. In addition, there is 
little experience to draw from in solving flow fields including viscous effects about swept-forward wing 
configurations. Hence, there were a number of compelling reasons to undertake the code calibration 
phase of this effort. 

The first task addressed was to generate a detailed computational model of the fuselage, figure 6. 
The objective was to make as few compromises as possible to the geometry while generating a model 
for which computational mappings could be generated by the codes. The wing and canard geometry 
has also been modeled; however, this has been a much more straightforward task. The compromises to 
the wing geometry consisted of a modified trailing edge in the vicinity of the body stake for each of the 
codes and modeling a wing without any dihedral in the CANTATA code because of its Cartesian grid 
structure. 

The actual comparisons on the configuration have been performed in a systematic manner. The 
approach will be described here and selected comparisons from the various steps will be presented in 
the following section. The initial computations were performed on the wing-fuselage combination using 
the TAWFIVE and CANTATA codes. These computations were also compared to available wind tunnel 
data on the wing-fuselage configuration. The differences between the computations and experimental 
data were evaluated in light of previous experience. Code-to-code differences were analyzed in light of 
known differences (full-potential vs. small-disturbance, conservative vs. non-conservative, etc.) between 
the codes on computational results. Comparisons have also been made between computations on the 
fuselage-canard geometry. Each of these sets of calculations and comparisons will help to understand 
the wing-canard interactions. These analyses will also allow the wing-fuselage-canard configuration to 
be analyzed with significantly greater confidence than if the full configuration were analyzed initially. 

Computations on the wing-fuselage-canard configuration in the CANTATA code have been com- 
pared between available wind tunnel data on the l/8-scale model and flight data. The available ex- 
perimental data was surveyed to determine if comparable test points existed between the wind tunnel 
and flight tests. Comparisons of these data were accomplished at appropriate conditions. Differences 
between these data were attributed to either aeroelasticity or Reynolds number. Calculations with the 
TAPS program were accomplished to investigate aeroelastic effects. These data were compared with 
deflection measurements from the flight system. 

The use of any computational technique benefits from the accuracy inherent in the method. How- 
ever, it is also realized that a code can be used effectively if one understands the computational limitations 
and the effects of any inaccuracies which may be present. This becomes increasingly important as one 
approaches the boundaries of practical applicability for a given code. The approach undertaken in this 
study has identified some practical limitations of the computations for this complex configuration. 


SELECTED COMPARISONS 

Initially, computations were performed on the X-29A wing-fuselage combination and compared 
with available wind tunnel data from the l/8-scale model tested at NAS A- Ames (ref, 4). Results are 
presented for a subsonic Mach number (M = 0.6) and a high- transonic Mach number (M = 0.9). The 
computational results are inviscid with no aeroelastic effects included. The transonic analyses were 
performed at the experimental angles of attack and Mach numbers. These results demonstrate the code- 
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to-code differences and cover the range of applicability for the codes up to flow conditions where the 
flow shows sign of significant separation on the wing. 

Results are compared at M = 0.6, a = 7.7°, in figure 7. The comparisons are quite good at 
this condition. Notice that generally the full-potential code (TAWFIVE) predicts more expansion at 
the leading edge on the upper surface. Because of the grid density at the trailing edge, the CANTATA 
results seem to be more sensitive to the flaperon geometry. Also note the differences in the computational 
predictions at the wingtip. The TAWFIVE code predicts the upper-surface leading-edge peak quite well. 
However, the small-disturbance results compare better with the experiment aft of the leading edge. 

The comparison presented in figure 8 are at M = 0.9 and a = 6.9°. The comparisons on the 
outboard part of the wing are fairly good, although the experimental results show trailing-edge (rj = 0.7) 
and shock-induced (17 = 0.49) separations. On the inboard portion of the wing, the lower surface pressures 
are predicted fairly well; however, the upper surface pressures show some discrepancies. At rj = 0.31 and 
0.19 the experimental and computational results on the upper surface forward of the shock do not match 
well. While the experimental data are sparse in this region, it does not appear that the characteristics 
of the flow are predicted with either computational technique. In general, the two codes predict similar 
flow patterns over the wing span for the flow conditions observed. Comparisons for the fuselage-canard 
configuration also showed consistent predictions between the two codes. 

The next step in the calibration effort involved computations on the wing-fuselage-canard configura- 
tion. The first significant incongruence between the actual configuration and the geometric capabilities 
of the CANTATA code was uncovered. The grid structure for both the canard and wing grids is Carte- 
sian. This results in the boundary conditions for both surfaces being applied in parallel planes. Because 
of the dihedral in the wing, the locus of wing leading- edge points intersects the canard wake. This 
in conjunction with the closely coupled lifting surfaces results in a significant interaction between the 
wing and canard. The comparison presented in figure 9 shows the effect of relative wing-canard vertical 
position on the wing pressure distribution near mid span. The calculations were made with the wing 
3 inches (full scale) above and below the plane of the canard. When the wing is positioned above the 
canard plane, the loading over the inboard wing sections is collapsed. By comparisons with wind tunnel 
experimental data, it became obvious that positioning the wing slightly below the canard yielded the 
most reasonable flow simulation. 

Comparisons between wind tunnel data from reference 4 and computational predictions for the 
complete configuration are included in figures 10 and 11. The comparisons are strikingly similar for 
these M = 0.9, low lift (figure 10) and M = 0.95, moderate lift (figure 11) conditions. The inboard 
region on the uppers surface is predicted adequately. As one observes the comparisons outboard of the 
leading-edge break, the characteristic of the experimental pressures on the upper surface is not captured 
in the computations. This difference was observed in all the comparisons at M = 0.9 and above. The 
pressure distribution on the canard is predicted reasonably well except at the leading edge on the upper 
surface. An expansion observed in the experimental data is not evident at all in the computations. This 
difference could be the result of geometric compromises related to the inlet or to a vortex originating 
from the forebody or inlet lip. The accentuated expansions and compressions near the wing trailing 
edge axe the result of computationally modeling the discrete hinge-lines for the segmented trailing-edge 
flap system. Based on these comparisons at the higher Mach numbers, plus comparisons at subsonic 
conditions, it was felt that an adequate computational model had been generated. 

At this point the focus on the experimental data shifted from the data obtained in the controlled 
wind tunnel environment to data obtained in the dynamic flight environment where the uncertainties 
axe increased. Before proceeding with comparisons between the computations and the flight test data, 
it was of interest to determine how well the flight and wind tunnel data compared. Because of the large 
number of parameters which had to be matched (Mach, angle of attack, canard, wing flap, and strake 
flap deflections, etc.), only a few conditions were identified for comparison between the wind tunnel 
test and flight test. Comparisons at M = 0.8 and M = 0.9 are presented in figure 12. At M = 0.8 
and a = 2.2°, the comparisons between flight and wind tunnel data are excellent at both of the span 
stations shown. At higher Mach numbers and angles of attack the comparisons were not quite as good. 
It was interesting to note, as shown in figure 12b, that a very good comparison was obtained when the 
tunnel data were compared at an elevated angle of attack relative to the flight data. Several hypotheses 



were purposed to explain this anomaly including wind tunnel wall interference, Reynolds number effects, 
and aeroelastic deformation. Wing deflection data were available at some flight conditions. A flight test 
point of interest was identified and the configuration was analyzed in TAPS. The resulting predicted 
deflections are compared with the measured deflections in figure 13. The measured deflections from the 
flight test also show an error range with the data. TAPS overpredicted the deflections along the entire 
span. However, the inclusion of the computational (overpredicted) deflections in the solution did not 
yield differences between the rigid and aeroelastic solutions as large as those observed between the flight 
and tunnel experiments at comparable angles of attack. These anomalies accentuate the need for the 
experiment in the NTF. 

The final set of comparisons is between computational predictions and flight test data at M = 0.6, 
0.8, and 0.91 and moderate lift levels (figures 14, 15, and 16). The comparisons at M = 0.6 and 0.8 are 
quite good on both the wing and canard. At the outboard span location on the wing the upper-surface 
leading-expansion is slightly underpredicted at M = 0.6. At the inboard location on the wing at M = 0.8 
the leading edge expansion is overpredicted in the computations. The upper-surface leading-edge flow 
is predicted adequately at M = 0.91 for the inboard station yet the level of expansion from 20% to 
30% chord is not predicted accurately. At M = 0,91 the characteristic of the upper-surface flow is 
not predicted at the outboard span location. The characteristics of the experimental and computational 
pressure distributions are similar to those observed in the comparisons of wind tunnel and computational 
predictions in figures 10 and 11. The influence of the rather bulbous flap-track fairings is evident in 
comparisons on the lower surface at the outboard span location for M = 0.8 and 0.91. There was no 
attempt to model these fairings computationally. 

There is evidence of flow expansion at the inboard canard location leading edge for all three Mach 
numbers which is not predicted computationally. At M = 0.91 the expansion appears to terminate in a 
shock at approximately 10% chord. While not shown in these figures, this expansion migrates outboard 
along the canard with increasing angle of attack. Except for this anomaly the canard flow is predicted 
well. 

Overall the comparisons between the experimental data and computational predictions are en- 
couraging. However, there were differences observed for which the causes were difficult to isolate. In 
particular, at higher Mach numbers the flow on the outboard part of the wing, as well as the canard 
leading edge was not predicted well computationally. These observations identify the need for obtaining 
high-Reynold8-number data under controlled conditions. 


SUMMARY AND CONCLUSIONS 

This paper has presented an overview of a major cooperative effort between NASA Ames/Dryden 
Flight Research Facility and NASA Langley in support of the X-29A flight-test program. The effort 
involves flight testing, wind tunnel testing in the NTF, and computational support. Each phase in the 
effort has distinct advantages and disadvantages relative to the data which are obtained. For example, 
the flight data could be considered “gospel” without scaling or wall interference effects. However, as was 
previously discussed, it is difficult to isolate the influences of individual parameters. Wind tunnel testing 
allows the component build-up and parametric variation of independent variables. Hence, individual 
influences and interference effects can be isolated. However, the data can suffer from scaling effects 
and/or wall interference. In addition, although the range of variables tested can cover the airplane 
capability, testing all combinations of each variable quickly expands to a prohibitively large matrix. 
While physical geometry and fluid physics modeling comparisons are often made, CFD allows the gaps 
between tunnel and flight data to be filled. It also allows evaluation of flight conditions outside the 
cleared flight envelope and estimation of the effects of configuration modifications. Hence, it is easy to 
see from this simplified discussion that each phase of the effort complements the other phases. 

An adequate computational representation of the X-29 A configuration has been developed. Com- 
parisons between computational predictions and both flight and wind tunnel experimental data have 
been made a a range of transonic conditions. At freestream Mach numbers less than 0.9, the comparison 
of pressure distribution between the computations and experiment is quite good on both the wing 
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and the canard. Anomalies have been observed for moderate lift coefficients at M = 0.9 and above. 
The computations show no evidence of the flow expansion observed on the canard leading-edge. In 
addition, the experimental pressure distribution over the mid-span of the wing has a different character 
than that predicted computationally. Further computational and experimental investigations of these 
discrepancies are required to understand these anomalies. 

The flight testing of the X-29A aircraft is well under way and a significant amount of data is becoming 
available. The NTF wind tunnel test model is virtually complete and the NTF test will commence during 
the Fall 1988. Computational support is continuing with the focus on understanding the anomalies 
which have been identified. The data available from the three-phase effort will greatly enhance the 
understanding of the complex flow phenomena and aerodynamics of the X-29A aircraft. 
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Figure 5 - Comparison of NTF and X-29A flight 
envelopes 



Figure 6 - Wire frame representation of the computational model 


Figure 8 - Comparison of computations and wind tunnel results M=0.9 t a = 6.9 4 










a) Wing above canard b) Wing below canard 

Figure 9 - Effect of wing/canard position on wing pressure distribution - M=0.9, a = 4.7° 
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b) Canard 


Figure 10 - Comparison of computational results and wind tunnel data - M=0.9, a = 2.1° 
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b) Canard 

Figure 11 - Comparison of computational results and wind tunnel data - M=0.95, a = 5.9' 
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Figure 12 - Comparison of wind tunnel and flight test data 
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Figure 13 - Comparison of measured and predicted wing deflections - M=0.92, a = 7.5° 
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SUMMARY 

Results of correlative and design studies for transition location, laminar 
and turbulent boundary-layer parameters, and wake drag for forward swept and 
aft swept wings are presented in this paper. These studies were performed with 
the use of newly developed, improved integral -type boundary-layer and 
transition-prediction methods. Theoretical predictions were compared with 
flight measurements at subsonic and transonic flow conditions for the variable 
aft swept wing F-14 aircraft for which experimental pressure distributions, 
transition locations, and turbulent boundary-layer velocity profiles have been 
measured. Flight data were available at three spanwise stations for several 
values of sweep, freestream unit Reynolds number, Mach numbers, and lift 
coefficients. Theory /experiment correlations indicate excellent agreement for 
both transition location and turbulent boundary-layer parameters. The results 
of parametric studies carried out during the design of a laminar glove for the 
forward swept wing X-29 aircraft are also presented. These studies included 
the effects of a spanwise pressure gradient on transition location and wake 
drag for several values of freestream Reynolds numbers at a freestream Mach 
number of 0.9. 


INTRODUCTION 

The rise in jet fuel prices during the last decade, and uncertainties of 
the fuel costs for the future have been the cause of continued worldwide 
interest in improving military and commercial aircraft fuel efficiency. 
Additionally, improvements in the range of military aircraft are additional 
factors for achieving aircraft fuel efficiency and drag reductions. The most 
promising aerodynamic means of achieving significant gains in aircraft fuel 
efficiency is to design wings that provide extensive regions of laminar flow 
(refs. 1 and 2). Laminarization of the boundary layer on the surface of an 
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aircraft wing can be accomplished by the use of concepts such as Natural 
Laminar Flow (NLF), Laminar Flow Control (LFC), and a combination of NLF and 
LFC which is referred to as Hybrid Laminar Flow Control (HLFC). 

Extensive regions of laminar flow can be maintained on a wing surface with 
the NLF concept by suitably tailoring the spanwise and chordwise pressure 
gradients to limit the amplification of disturbances that trigger the 
transition of the laminar boundary layer to a turbulent one. Flight tests at 
the Dryden Flight Research Facility on a variable sweep TACT F-lll Fighter 
aircraft (ref. 3) with an NLF wing glove have demonstrated that natural laminar 
flow can be maintained over significant areas at transonic speeds. Recent 
flight tests of a laminar clean-up glove (a thin, smooth fiberglass coating) on 
an F-14 aircraft with a variable-sweep wing (ref. 4), equipped with 
instrumentation for the accurate determination of transition, also revealed 
similar phenomena. Correlative studies were performed during the present study 
for the F-14 aircraft laminar wing, clean-up glove using the methods of 
references 5 and 6 extended to account for spanwise pressure gradient 
effects. Results of these studies are presented in this paper in the form of 
theory /experiment correlations for transition location and turbulent boundary- 
layer parameters. In addition, comparisons of the hot-film experimental data 
on the F-14 aircraft's wing glove with theoretical predictions are described. 

The design of either forward- or aft-swept laminar flow wings at subsonic, 
transonic, or supersonic speeds, using NLF, HLFC, or LFC concepts requires the 
use of reliable theoretical computational methods. A number of finite- 
difference boundary-layer methods, stability methods, and full Navier-Stokes 
equation solvers are available; however, these methods are not well suited for 
the large number of calculations required for parametric wing design and 
optimization studies. Several new theoretical integral -boundary-layer methods 
(references 5, 6, 7, and 8) have small computer execution times and are very 
simple to use. These codes have the capability of predicting the laminar, 
transitional, and separating turbulent boundary layers for applications at 
subsonic, transonic, and supersonic Mach numbers. 

The design studies reported in this paper involve the design of an NLF 
glove for the forward swept X-29 wing. The methods of references 5, 6, 7, and 
8 have been modified for the present theoretical investigations to allow 
consideration of the effects of spanwise pressure gradient. Results are 
presented utilizing these modified codes for various spanwise pressure 
gradients and a prescribed streamwise pressure distribution to determine the 
extent of laminar flow on the corresponding X-29 wing glove. The wing 
sections' upper-surface geometrical changes required to obtain the prescribed 
pressure distribution are also given. Furthermore, the present design- 
optimization studies show qualitatively the procedure for optimizing values of 
spanwise gradients so as to minimize the wake drag with due consideration of 
the effect of imbedded shocks on the boundary layer. Finally a comparison of 
the extent of laminar flow on "equivalent" forward-swept and aft-swept wing 
configurations is made. 
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NOMENCLATURE 


AR 

ASW 


D i 


J wake 


C p 


sonic 


aspect ratio 
aft swept wing 

chord p 

r 

induced drag coefficient - - ^ — 

wake drag coefficient 

pressure coefficient or specific 


heat of air at constant 


pressure 

pressure coefficient corresponding to local Mach number 
of 1.0 


AC P or &C p 
K 1 c 


c 2 


r 


e 

FSW 

H 

L, M, N 

K, X 

L. E. 

M 

N.F. 


defined in figure 18 

indicator of either growth or decay of disturbance for U 
component of velocity in transformed plane 

downstream rate of propagation of disturbance for U 
component of velocity in transformed plane 

indicator of either growth or decay of disturbance for W 
component of velocity in transformed plane 

spanwise rate of propagation of disturbance of W component 
of velocity in the transformed plane 

wing efficiency factor for induced drag 

forward swept wing 
* 

form factor = |- 

groups of dimensionless parameters in transformed planes 
presenting the effects of suction, pressure gradient and 
curvature of velocity profiles. Each symbol defined in 
text. 

leading edge 

Mach number (used only with subscripts) 

transition or instability due to normal flow 


P 


static pressure 


Prandtl number 
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q 


Qx 

or 

Q z 


R 


R T,n 

R inst 

R tran 

s,Z,z 


S.F. 


T 


T.E. 


u,v,w 


U.V.M 


x,y 

x, Y 


dynamic head 

dimensionless parameter representing effects of suction and 
pressure gradient on the shape of velocity profile in 
transformed plane 


6 t c dU Q 

T,s e 
v_ dX 


Vt,s . _ 6 T,z 

v o ’ Qz ' % 


dW^ 

dZ" 


V .5 


s T,z 


freestream Reynolds number based on chord 


.( 


Ue^ 
or 
We/°T,n 


** 


m 


value of Ry n at neutral point 
value of Ry >n at transition 


U c 

oo 

V 

00 


curvilinear coordinates in physical plane; s is normal to 
local sweep lines on wing surface, £ normal to wing 
surface, and z normal to s-£ plane in spanwise direction 


transition or instability due to spanwise flow 
temperature, °R 


trailing edge 

velocity components in physical planes in s, 5 and z 
directions, respectively 


velocity components in transformed planes in X, Y, and Z 
directions, respectively ^ 

wing section coordinates in streamwise direction 

Stewartson's transformed coordinates 

5 T 


s a P 

x = / « ~ ds; 
0 0 0 


d p 

e K e 


■rr I 


~0 P 0 0 


*1 disturbance stream function for transformed plane X-Y 

^2 disturbance stream function for transformed plane Z-Y 

i'y disturbance amplitude function for X direction in 

(y) transformed plane 
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(y) 


u 

a 

v 

0 


*T, X 


di sturbance amplitude function for Z direction in 
transformed plane 

2Tr/(wavelength of disturbance in transformed plane in X- 
di recti on) 

2ir/(wavelength of disturbance in transformed plane in Z- 
di recti on) 

dynamic viscosity 

proportionality constant for viscosity -temperature 
relationship 

kinematic viscosity 

momentum thickness in physical s, 5 plane for 2-D or 3-D 

00 u 

normal direction 0 = / p |l- (1 - 

0 p e e e 

momentum thickness in transformed X, Y plane for 2-D or 3-D 

-/Tc 1 -^ dY 

0 e e 


’t.z 


momentum thickness in transformed Z, Y plane 

■ r 0 T 11- r> dY 

e e 


p,s 


/ P ’ S rir C 1 - fci* 


e e 




r ri - ^)d 5 

i p e w e ^ V 


0 


P»n 


/ p,n fr 11 - S- )d « 


e e 


6 greater of 6^ _ or 6 
p,n a p,s p,z 
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6 


boundary-layer thickness in physical plane 


* 

6 


■r 

0 


(1 - 


pu 

P u 
p e e 


)d£ 


boundary-layer thickness in transformed plane 


* 


6 


T,X 


* 


6 


T,Z 




T,X 


'k'k 


T,Z 


A 

Y 

n 

Subscripts 

o 


e,s 


e,z 

aw 

m 

N 

P 


displacement thickness in transformed X, Y plane for 
2-D or 3-D 

/ t 1 - -0- ) dv 

o e 

displacement thickness thickness in transformed Z, Y plane 

■ r a- /idv 

o w e 

/ T ’ n (1 - f ^ dY 
o e 

/ T,n (1- /) dY 
o e 

s T n greater of « T>X or s T _ z 
local sweep angle in degree 
ratio of specific heats for air 
= 2y/b 

stagnation condition 
edge of boundary layer in s-di recti on 
edge of boundary layer in z-direction 
adiabatic wall condition 

mean or average value over boundary-layer thickness 
normal section 
physical plane 
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P,S 


p,z 


p,n 

s 

T 

T,X 

T,Z 

w 

X 

Z 

L.E. 

T.E. 

6 

0 


S 


physical plane in s-di recti on 
physical plane in z-di recti on 
either p,s or p,z 

direction normal to local sweep lines 

transformed plane 

transformed plane X-di recti on 

transformed plane in Z-direction 

conditions at wall 

transformed X-di recti on 

transformed Z-direction 

conditions at leading edge 

conditions at trailing edge 

parameter based on boundary-layer thickness 

parameter based on momentum thickness 

freestream value 

measured in streamwise direction 


In addition to the above list of symbols and subscripts, symbols and 
subscripts are also defined at appropriate places in the text. 


DESCRIPTION OF THEORY 

A description of the theory for the integral methods used to compute the 
laminar boundary layer, transition, and the turbulent boundary layer on a swept 
wing is given in references 5, 7, and 8. These theoretical methods were 
developed for infinite-span swept wings at subsonic and transonic speeds, and 
for two-dimensional airfoils and bodies of revolution at supersonic speeds. 
However, they did not include the effects of spanwise pressure gradient and the 
resulting effects on the location of transition. The boundary-layer "wash-out" 
phenomena are present on finite aft swept wing configurations with a given 
taper and twist, whereas the boundary-layer "wash-in" phenomena are present on 
configurations with forward sweep. The above theoretical methods of references 
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5, 6, and 7 have been extended in the present paper to account for the effects 
of spanwise pressure gradients on laminar boundary-layer transition. At the 
present time, the effects of suction on the extended method are not considered. 

The theoretical programs, which have been developed and used in the 
present theoretical investigations of the wing glove of the F-14 aircraft and 
for the parametric design studies of the forward swept wing of the X-29 
aircraft with a laminar flow glove, are classified and enumerated as follows: 

1. Potential /Viscous Design and Analysis 

a. Perturbation Method of Characteristics for inverse design at transonic 
Mach numbers. 

b. Bauer-Garabedian-Korn Program for transonic analysis. 

c. NASA-Lockheed Multi -component Airfoil Program for subsonic analysis. 

2. Boundary-Layer Analysis 

a. Integral Compressible Laminar Boundary-Layer Method with arbitrary 
chordwise pressure distribution and spanwise pressure gradient, sweep 
and suction at subsonic through high supersonic speeds. 

b. Short bubble and reattachment criteria. 

c. Instability and transition prediction due to leading-edge-normal and 
spanwise flow including the effects of spanwise pressure gradient. 

d. Separating turbulent boundary-layer method from subsonic to transonic 
speeds. 

3. Method for computing C n for an infinitely swept wing from subsonic 

u wake 

through supersonic speeds. 

4. Modified strip method to account for the effects of taper, sweep, and twist 
on wing-section characteristics. 


The theoretical development of the laminar boundary-layer method for an 
infinitely swept wing will be briefly described. In addition, the methodology 
for determining the effects of spanwise pressure gradients on laminar boundary- 
layer growth and on transition location will be outlined. 


Governing Equations 


The usual governing equations" for a compressible hydrodynamic boundary 
layer on an infinitely swept wing for the coordinate axis system shown in 
figure 1 are the continuity, streamwise momentum in s-direction, and momentum 
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in the 5 and z directions. The coordinate axis of figure 1 is chosen for the 
purpose of facilitating the application of the present theory to the design of 
the wing sections for an arbitrary wing planform. Moreover, the determination 
of the effect of the spanwise pressure gradient on the extent of laminar flow 
on swept wings by the present methods is greatly facilitated by the use of the 
coordinate axis system selected. Thus, 


d( pu ) + 9 ( pV ) = 0 Continuity (la) 

3s 35 


3u . 
pU 3? + 


pv 


3U 


|4 + b (p If) * (s-direction 
9s 9 ^ 9 ^ momentum) 


3u ■ 


(lb) 


1L = 0 (Normal Momentum) 

35 


3W „ 3W _ 3 f 3W>* 

pU 8s + PV 35 ' + 3? 35 J 


(z-di recti on 
spanwise 
momentum) 


(lc) 

(l d) 


The above equations (1) contain the terms consisting of variable physical 
properties, such as density p and dynamic viscosity y. Variations of these 
properties across the boundary layer as well as along the flow direction are 
not negligible for laminar boundary-layer flow at transonic speeds and are 
accounted for. At supersonic speeds, variations in the physical properties are 
very large which means that there is a strong coupling between solutions of the 
hydrodynamic and thermal boundary-layer equations. 


In order to simplify the governing equations (1) for solution by integral 
techniques while maintaining realistic computational results for the 
hydrodynamic and thermal boundary layers, the use of the following Stewartson's 
transformation is made: 


X = 


s a P 

r G G j 

I a J~ F" ds 

0 a o 0 


Y = 


a p 
e e 


5 T, 


f l r dc 

^0 0 


(2) 
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where ot is the proportionally constant for the viscosity-temperature 
relationship, namely 



T awV /2 T o + 198 ' 6 
\) T aw + 19ti - 6 


(3a) 

(3b) 



(1 + /Pr ^ M 2 ) 

ITT^f 


(3c) 


Furthermore, if it is assumed that the viscous flow on the wing can be 
divided into a finite number of suitably oriented strips on its surface, and 
that each infinitesimal control volume of the strip is a portion of an 
infinite-span swept wing, then the governing equations of motion in the 
physical and transformed planes can be further simplified. The transformed 
boundary-layer equations using Stewartson 's transformations, for infinite-span 
swept wing conditions, which can be derived from equations (1) with adiabatic 
wall temperature, are written as, 


M + iX = o 
ax aY u 


au , v au _ .. dU e . 
u ax + v W u e dx“ + v c 


a 2 u 


aY 


n aw v aw _ a 2 w 
u ax + v aT * v o ^2 


(4a) 

(4b) 

(4c) 


The relationships between the velocities within the boundary layer and the 
pressure gradients in the physical and transformed planes can be written 
(Pr « 1) as 


U=-^u; V=-^v (for adiabatic wall ) ; W = w 

cl a _ 


dU , dM 
H5T = a a o di~ t 1 


+ ^ M 



(5a) 

(5b) 


C*3 
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while the applicable boundary conditions in the transformed planes are given 
by. 


and 


0Y = O + U = W = O, 



0 Y 


>T,X 


♦ U = U = a M , 
e o e 


a 2 U _ n 

TTv = o ~ U 


all 

aY 


7 ? 


(6a) 

(6b) 


@ Y = 5 t>z . W = W e = w e> 


aw a 2 w 

9Y "^ 


o 


(6c) 


The compatibility conditions at 
planes can be written as, 

0 Y = 0 -*■ - 


the surface of the wing in the transformed 


w _ ii e j_ r a 

v sWw " u e dx v o( „2' 


ft 

aY'- w 


(7a) 


0 Y = 0 + 


r a 2 u 


r a 3 iK 
V - jJ 
0 3Y J w 


(7b) 


and in the absence of suction for flows with spanwise pressure gradients, 


1 aP 
p az 


Y=o 



(7c) 


The integral equations in the transformed planes in the leading-edge normal 
direction X and in the spanwise direction Z can be derived by making use of 
equations (4), boundary conditions (6), and wall compatibility conditions 
equations (7), and by use of Leibnitz's rule. 

Momentum Integral Equation in X-Direction, Infinite Swept Wing Conditions 


U e dX (dlWdx) = 2 [ L ' k(h T,s + 2) " S e J 

v- A 


( 8 ) 
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Momentum Integral Equation in Z-Direction, Infinite Swept Wing Conditions 


3* [«, “ e / T ’ n I (1 - S-)«] 

0 e e 


■ 3W 


v o^3Y^Y=0 " V s W e (9) 


where 6j n is greater of 6 T z or <Sj x - 


The symbols L, K, and S Q appearing in the above equation (8) are defined 
in equation (20) for transformed X and Z coordinates. 

The velocity profile assumptions for the flow in the X- and Z-di recti ons 
in the transformed planes are given by the following equations, 


where 


U _ . » 2 . 3.ii 4 

U A l n X + A 2 n X A 3 n X A 4 n X 


n X = 6 


T S X 


and 


A^ = 2 + Q^/6; A z - -0.5Q^; 

A^ ~2 — A 2 , A^ 3 A | 






(10a) 


(10b) 


(ID 


where, n z = 5 


T,Z 


The numerical solutions of equations (8) and (9) are performed by the 
Euler method with repeated iterations. These solutions give the shape of 
velocity profiles jj- and in the transformed planes under the assumption of 

e e 

infinite-span swept-wing conditions. 
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EFFECTS OF BOUNDARY-LAYER MASH- IN/WASH-OUT DUE TO SPANWISE 
PRESSURE GRADIENTS ON TRANSITION ON FINITE, SWEPT, 
TAPERED, AND TWISTED WING 


The results of theory /experiment correlations for several wings using both 
wind tunnel and flight-test data have indicated that for sweeps greater than 
approximately 25°, and freestream Reynolds numbers greater than 20 million 
transition is triggered by spanwise flow in the z direction. Thus, for the 
purpose of designing laminar-flow wings and gloves at transonic and supersonic 
speeds at Reynolds numbers and sweeps of practical interest, it is necessary to 
consider the effects of boundary-layer wash-in/wash-out due to spanwise 
pressure gradients. This is especially important for wings of supersonic 
aircraft where the sweep of the leading edge is usually in excess of 45°. In 
order to account for the effects of boundary-layer wash-in/wash-out on finite 
swept wings the approach shown schematically in figure 2 is used. This figure 
shows the infinitesimal control volume composed of surfaces of a trapezoid. 

The coordinate system is the orthogonal s, £ and z axes where the s 
coordinate is normal to the local sweep line, z-coordinate is normal to s in 
the plane of wing, and the £-coordinate is normal to the s-z plane. 
and are constant s/c lines where c is the local chord in s-direction 

and A^Dj and are parallel to s-direction. 


The momentum theorem is applied in the z-di recti on, which states that the 
algebraic sum of the forces acting on the control volume ABCD - A^BiC^Di in 
the z-di recti on is equal to the rate of change of momentum in control volume 
plus the net flux of momentum in the z-di recti on across the control surfaces of 
the control volume. For steady mean flow conditions an equation can be derived 
for the boundary-layer wash-in/wash-out due to spanwise pressure gradient 


effects on velocity profile — in the z-di recti on. This equation in the 


w 


physical coordinate system for transformed boundary-layer quantities can be 


written as: 
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where 


(0y Z^N = Transformed momentum thickness in Z-direction 

* which is modified by the boundary-layer wash -in/ 

wash-out phenomena due to spanwise pressure 
gradients on finite-swept and tapered wing. 


(e 


T,Z 


o 


Transformed momentum thickness in Z-direction for 
infinitely swept wing assumption conditions. 


A 


X 


(d x ) 0 — 

T,X 0 

0 



Dimensionless pressure gradient parameter in 
X-direction in transformed plane for infinite-span, 
swept-wing assumption conditions. 


Transformed boundary-layer thickness in X-direction 
for infinite-span swept-wing assumption condition. 

T o + due Sutherland 

t + 198.6 theory of viscosity. 

aW 

The following approach is used for the purpose of calculating (<$ T _) and 
the velocity profiles in the transformed plane in the spanwise direction in 
order to account for the boundary-layer wash-in or wash-out phenomena. 

Define the quantities, 


and. 


2 1 dW e 

A 7 = (6 T 7 ) w ( — -r4 ) 
Z N T,Z N v Q d Z 

? i dW 

h = (0 T 7 >KJ (— ht^) 

Zfj T,z N v 0 dZ 


(13) 



where: 

( 6y Z ) N = Transformed boundary-layer thickness in Z direction which is 

’ modified by the boundary-layer wash-in/wash-out phenoena due 

to spanwise pressure gradients on finite-swept and tapered 
wing. 

If a fourth order velocity profile is assumed in the spanwise Z direction 
in the transformed plane, then the relation between A and a can be derived 

as the following polynomial: N N 


180 



where. 


A 7 = P. + P ? 8 7 + Po 6 2 7 + P,8y + Pj-By + PfiSy 
Z N 1 2 Z N 3 Z N 4 Z N 5 Z N 6 Z N 


Pj_ = 0.032015 
P 2 = 73.156 
P 3 = 10.587 
P 4 = -129.6 
P 5 = 4669.6 
P 6 = 43865.0 


thus, ( S T can be written as, (6 T 

I ,4 N 1,41 


Z ^ N ^ ^ 

, 4 n 4 n 


1/2 


W 

The velocity profile in the transformed plane can therefore be calculated 

e 

from the following assumed fourth order polynomial given by: 


W B l n Z + B 2 n Z + B 3 n Z + B 4 n Z 
e 


where. 


n, = 


(6 T,Z ) N 


and Bj = 2 + -g-; B 2 


7 A z. 


B 3 = -2 - B 2 ; B 4 = 3 


THERMAL BOUNDARY-LAYER SOLUTION 

Temperature variations within the boundary layer at subsonic speeds are 
small, hence the effects of temperature variations within the boundary layer on 
instability and transition at subsonic speeds are neglected. However, at high 
transonic and supersonic speeds, the heat generated by friction and adiabatic 
compression is quite significant. These phenomena give rise to large 
temperature variations across the boundary layer, and it has been found from 
computational experiments that these large temperature variations affect the 
following equations for computing temperature profiles; they can be derived 
using the procedures described earlier. 
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For determining the temperature profile at a given location on the wing 
use is made of either equations (14a) or (14b) depending upon whether 5-,- x is 

U W 

larger than 5 X , or vice versa. The velocity profiles jr- or y- in the 
T,Z u e w e 

transformed planes in X or Z direction, which are required in equations (14a) or 
(14b), are dependent on local pressure gradients in the leading edge normal -s 
or spanwise-z directions due to use of equations (9) or (11), respectively. 

Thus, the effect of pressure gradient and suction is implicitly accounted for 
in the determination of temperature profiles by the present method. 


Criteria for Determining Locations of Neutral Instabilities 
for Velocity Profiles in s and z Directions 

The Orr-Sommerfeld equations in transformed planes for X and Z 
directions can be written as follows: For the transformed coordinate in the X- 

di recti on, 

(U - Cj)^' - o 2 ^) - + aj^] (15) 

1 0 T,X 

and for the transformed coordinate in the Z-direction, 

( W - Cg ) ( 4>2 ' " a 2^2 ^ ~ ^yY^2 = a fT - — " ^ a 2^2. ' + ^ 

2 0 T,X 

The disturbance stream functions for the transformed X and Z directions 
are assumed, as for the physical plane, as the following types: 

ia, (X-C, t) ■ , jn: 

= <t>i( Y ) e 1 + X-direction (17a) 


ia 2 (Z-C 2 t) 

^2 = ^2(Y) e * Z-direction (17b) 
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where. 


$ im and <j)p , Y \ = disturbance amplitude function 
11 Y) for X and Z direction, respectively 

in transformed planes 

a, and a ? = 2rr / ( di sturbance wavelength in transformed 
1 1 planes in X- and Z-direction, respectively) 

X = transformed coordinate for s-di recti on 
Z = transformed coordinate for z-direction 

C 1 = C. + iC 1 
r i 

complex in general 

Cp = C? + 
c c r S' 

or =0 corresponds to neutral instability 

> 0 corresponds to amplification of disturbance 
< 0 corresponds to decay of disturbance 


The generalized solution between the dimensionless pressure gradient 
parameters and "equivalent" Reynolds number based on momentum thickness for the 
neutral stability of laminar boundary layers in either s or z direction is 
shown in figure 3. The curve shown in this figure is derived from solutions of 
the Orr-Sommerfel d equations, equations (15) and (16), in the transformed X-Y 
plane or Z-Y plane for appropriate fourth-order polynomial velocity profiles. 
The effects of pressure gradient in the s-di recti on or z-direction, local Mach 
number, temperature profile, and suction have been accounted for in the 
derivation of this curve. The abscissa of figure 3 is given as the zero 
suction ( S Q = 0) version of 


(0 dM N 
P.s e,s 
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where 


v m 


kl 


vd£ 



( 19 ) 


In equation (19), the term (t^— ) or ( v - T - ) 

e,s e,z 

substituted from equations (14a) or (14b). 


for temperature profile are 


It should be noted that because of the choice of the coordinate axis 
systems shown in figures 1 and 2, it is possible to use the same criteria shown 
in figure 3 for determining the neutral stability point for disturbances in 
both the normal and spanwise directions. 


Boundary -Layer Transition Criteria 


The criteria used in the present method for determining the transition 
location due to the amplification of disturbances either in the normal or 
spanwise directions are shown in figure 4. The parametric curves of figure 4 
are valid for arbitrary specified pressure gradients - both in the normal as 
well as spanwise directions. The independent arid dependent variables shown in 
this figure contain groups of several hydrodynamic and thermal boundary -layer 
parameters whicfTare significant during the transition process. In deriving 
the curve for Jow values of freestream turbulence (less than .05 percent), the 
use of experimental data obtained in flight or low- turbulence wind tunnels is 
made and this information is put in the transformed planes. The dimensionless 
parameters shown in figure 4 have been derived from dimensional analyses and by 
making use of the laws of dynamic similarity. These physical dimensionless 
parameters consist of the following: 
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In deriving the curves for several higher values of freestream turbulence 
greater than 0.05% shown in figure 4, use has been made of information 
presented in figure 5 for a zero pressure gradient and for values of Mach 
number from low subsonic to high supersonic (up to M ro = 4.0). It has been 
hypothesized in constructing the curves of figure 4 for values of freestream 
turbulence greater than 0.05 percent, that the effect of pressure gradient for 
freestream turbulence greater than 0.05 is similar to low values of freestream 
turbulence of 0.05. Curves for freestream turbulence higher than 0.05 percent, 
shown in figure 4, have been constructed by making use of information in figure 
5 and using the above hypothesis. 

It needs to be pointed out here that the curves of figure 5 can also be 
used for other types of disturbances, e.g., noise generated due to the 
turbulent boundary layer on the walls of the test section of supersonic wind 
tunnels. This objective can be accomplished after the application of 
appropriate calibration procedures. 


Turbulent Boundary Layer 

It has been known from experimental data on sections of NLF wings that if 
turbulent separation is present upstream of x/c rj 0.95, then the design 
pressure distribution necessary to maintain laminar flow is degraded. 
Additionally, the wake drag will be quite high due to presence of the separated 
flow region. In the present paper turbulent boundary-layer separation is 
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accounted for by the method developed in reference 8. This turbulent boundary- 
layer method has already been successfully applied in the design of NLF and 
HLFC wing sections and for the prediction of wake drag. 

To compute the turbulent boundary layer from the location of transition at 
various spanwise stations to the wing trailing edge, the method of reference 8 
was modified to account for effect of spanwise flow assuming infinite swept 
wing conditions. The important objectives for the computations of the 
turbulent boundary layer are: (1) to see if turbulent separation is present 

downstream of the transition location, and (2) to determine the turbulent 
boundary-layer quantities at the wing trailing edge. 

The theoretical method of reference 8 has been developed by analyzing the 
mean and fluctuating components of velocity profile data obtained by the use of 
a laser velocimeter for several airfoils with trailing-edge seperation. These 
experimental data indicated that as the turbulent boundary layer approached the 
separation location, the fluctuating velocity components become of the same 
order of magnitude as the mean velocity components within the turbulent 
boundary layer. This observation suggested the retention of the fluctuating- 
velocity components in the governing equations of motions as was done in 
reference 8. This method has been found to be extremely reliable in predicting 
turbulent boundary-layer quantities, including separation, on wing sections 
with thickness ratios in the range of 0.09 < t/c < 0.22 and Reynolds numbers 
from 2 to 30 million. 


Computations of Cn for Ming Sections 
"wake 

The derivation of the expression for C D for an infinitely swept wing 

W3K6 

is given in Appendix A. This expression, which is valid from subsonic to 
supersonic speeds, is given as 


wake 


c f^uuV.E. (p^T.E. ^ cos2a T.E. 

OO 00 

w 

* ( 8 uw)t.E. (r>T.E. 


where the exponent p is given by 

p = [H t>e _ + 5 + ( y -1)?]/2.0 (21) 

The C °wake ^° r 3 flnite ta P ered swept wing can be computed from equation 
(21) and by making use of modified strip theory. 
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APPLICATIONS OF THE PRESENT THEORETICAL JCTHODS 

The theoretical methods described in previous sections are utilized to 
perform two tasks in the following section, namely: 

(i) The predition of transition location and turbulent boundary-layer 
parameters on the "clean-up" glove of the F-14 aircraft and 
correlation with experimental data. 

(ii) Performing parametric studies for the design of a laminar glove for 
the forward swept wing of X-29 aircraft. 


Correlation with Experimental Data on the Glove of F-14 Aircraft 


The planform and instrumentation for the gloved F-14 wing with 20 sweep 
is shown in figure 6. With the "clean-up" glove installed, Q the sweep of the 
wing leading edge was varied from A L>£> = 20 to A L ^ = 35 , and the flaps an 

slats were locked in the up position! 'as shown in figure 6, the three rows of 
surface pressure orifices are located at spanwise stations n = 0.515, O.b/, an 
0.83 respectively, and were oriented in a direction parallel to freestream or 
the wing leading edge sweep A, E = 20°. In order to avoid interference 
between hot-film gauges, they were placed at an angle of 30 to rows of surface 
pressure orifices as shown in this figure. The boundary-layer rakes were 
installed at two spanwise locations at x/c = 0.55 in order to measure the 
turbulent boundary-layer velocity profiles. Transition and turbulent boundary- 
layer data were obtained for freestream Mach numbers from 0.6 to 0.8, wigg 
leading edge sweeps from 20° to 35°, and a Reynolds number range of 8x10 to 
25x10 . Pressure distributions for the aforementioned three spanwise 
locations were also measured for these Mach numbers and sweeps. These pressure 
distributions constituted the boundary conditions for computing both the laminar 
and turbulent boundary layers and the transition location on the F-14 glove. 


Figure 7(a) shows the pressure distributions on the upper surface of the 
qlove at three spanwise locations and hot film traces at station 2 
(n = 0.67) for the wing leading edge sweep of 20°, = 0.7 and a freestream 

Reynolds number based on chord of 13 million. In addition to using the 
pressure distributions as boundary conditions, they are used to evaluate the 
effect of spanwise pressure gradient on the development of the laminar 
boundary layer and transition location at station 2. Since the experimental 
pressure data only extended back to x/c = 0.6 theoretically computed pressures 
distributions were used from x/c = 0.6 to trailing edge in order to compute 
the growth of turbulent boundary layers from the transition location to the 
trailing edge. Figure 7(b) shows the computed results for chordwise 
transition and instability location due to normal flow (N.F.) and spanwise 
flow (S.F.) phenomena. The results are plotted in this figure as x/c of 
instability and transition versus freestream Reynolds number. These results 
indicate that at a Reynolds number of 13 million, the computed transition 
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location is due to N.F . phenomena and occurs at x/c = 0.33. The hot film 
trace shown in figure 7(a) shows that the experimental transition location at 
station 2 is between 0.3 < x/c < 0. 4. Thus, good agreement is obtained 
between the theoretically computed transition location and the hot-film- 
measured transition location. Figure 7(c) shows the results of computations 
for displacement thickness, S*/c, momentum thickness, e/c, and form factor H, 
plotted against x/c for station 2. Freestream Mach number, Reynolds number, 
and leading edge sweep for the results shown in figure 7(c) are the same as 
for figures 7(a) and 7(b). Experimental data for 6*/c, e/c, and H at x/c = 
0.55 for station 2 in the turbulent boundary-layer region are also plotted in 
figure 7(c). These results indicate that present theoretical computations for 
the turbulent boundary-layer parameters agree quite well with the 
experimentally measured quantities. Furthermore, these results indicate that 
the location of transition predicted by the present theory must be in close 
agreement with experimentally determined transition. It is also interesting 
to notice from figure 7(c) that the slope of e/c versus x/c in the laminar 
flow region is smaller than that in the turbulent region. This large slope of 
0/c versus x/c in the turbulent boundary-layer region will lead, of course, to 
large values of 0/c at the wing trailing edge. The wake drag coefficient, 
C °wake’ iS dlrectly P^portional to the momentum thickness, 0/c, at the 

trailing edge. 

Figures 8(a) and 8(b) show plots of experimental pressure distributions 
and hot film traces for = 0.801, a leading edge sweep of 20® and freestream 
Reynolds numbers of 11 million and 23 million, respectively . As can be seen 
from the comparison of pressure distributions in figures 8(a) and 8(b), there 
is not any noticeable difference in pressure distributions between Reynolds 
numbers of 11 and 23 million, and hence, the same pressure distributions at 
three spanwise stations were used for computations of transition and 
instability results for the range of Reynolds number shown in figure 8(c). 

These pressure distributions indicate that the shock strength is high, and for 
this reason, there seems to be some discrepancy in computed and measured 
turbulent boundary-layer parameters at x/c = 0.55, as shown in figure 8(e), 
even though the transition location is predicted quite accurately. The present 
shock boundary-layer model is being refined in order to provide better initial 
conditions for turbulent boundary-layer computations downstream of a strong 
shock. 


As seen from figure 8(c), the computed transition location for a 
freestream Reynolds number of 11 million is at x/c ~ 0.5 and caused by N.F. 
phenomena, however, at a freestream Reynolds number of 23 million the computed 
results indicate that transition is caused by S.F. phenomena and occurs at 
x/c = 0.41. Experimental data from the hot film sensor shown in figure 8(a) 
indicate that at a Reynolds number of 11 million, transition occurs at 
x/c ~ 0.5. Hot film traces at a Reynolds number of 23 million, which are 
shown in figure 8(b), indicate that transition is occurring, for this 
freestream condition, at x/c ** 0.4. Thus, as seen from comparisons of 
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computed results, shown in figure 8(c), with hot film traces of figures 8(a) 
and 8(b), the present theory reliably predicts both the location of transition 
as well as its mode, i.e., due to N.F. or S.F. phenomena. Figure 8(d) shows a 
plot of computed C n as a function of freestream Reynolds number. As seen 

from this figure, C n decreases as the freestream Reynolds number increases 

u wake 

from 8 to 20 million and above the Reynolds number of 20 million 

increase with Reynolds number. The Reynolds number of 20 million, where the 

computed Cn becomes minimum, corresponds to the initial triggering of 
u wake 

transition due to spanwise flow. This is an important phenomena and one that 
can be used to advantage as seen in the design of a proposed laminar glove for 
X-29 aircraft discussed subsequently in this paper. Finally, figure 8(e) 
shows the plots of computed values of 5*/c, e/c and H versus x/c for a 
Reynolds number of 12.4x10^. 

Figure 9(a) shows the experimental pressure distributions and hot film 
traces for a wing-glove leading-edge sweep of 25°, a freestream Mach number of 
0.706 and a Reynolds number of 11 million. By using the experimental pressure 
distributions as boundary conditions (see figure 9(a)), the transition and 
instability locations are computed and the results are shown plotted in figure 
9(b) as a function of freestream Reynolds number. Experimental data for 
transition obtained from hot film traces, shown in figure 9(a), indicate that 
transition for a wing glove leading edge sweep of 25.3°, Reynolds number of 11 
million, and M = 0.706 takes place for 0.4 < x/c < 0.5. Theoretical results 
for these same°°conditions shown in figure 9(b) also show that transition 
should occur at an x/c between 0.4 to 0.5. 

Now, it is generally believed that if the freestream Mach number is kept 
constant, that increasing wing sweep reduces the extent of laminar flow. 
However, experimental data and calculations for the wing-glove leading-edge 
sweep of 20° and M^ = 0.7, (figures 7(a) and 7(b)) when compared to results 
for a sweep of 25°T figures 9(a) and 9(b), show that transition is located 
closer to the leading edge for the lower sweep. Furthermore, comparison of 
the calculated results of figure 9(b) with those of figure 7(b) indicates the 
same phenomena. This is an unusual result and shows the ability of the 
present methods to optimize wing section geometries and planforms for 
extensive laminar flow regions at specified freestream conditions. 

Figure 9(c) shows the results of theoretical computations for the laminar 
and turbulent boundary-layer parameters 5*/c, 6/c and H plotted versus x/c for 
a leading-edge sweep of 25°, M of 0.704 and a Reynolds number of 
approximately 11 million. For°°this case, strong shocks are not present at any 
spanwise station, as can be seen from figure 9(a). The computed boundary- 
layer parameters agree quite well in the turbulent boundary-layer region at 
x/c = 0.55. This agreement of computations with experiment in the turbulent 
boundary-layer region is an indication that the computed transition location 
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is fairly close to the experimental transition location. Moreover, these 
results indicate the reliability of the present theoretical turbulent 
boundary-layer method. 

Figure 10(a) shows the experimental pressure distributions and hot film 
traces for a leading-edge sweep of 29.7°, M^ = 0.704, and a Reynolds number of 
9.5 million. The corresponding experimental data for a Reynolds number of 
19.2 million and leading edge sweep of 31.4° are shown in figure 10(b). 

Figure 10(c) shows theoretical results for instability and transition 
locations for a leading-edge sweep of 30° and = 0.704. These results show 
that for a Reynolds number of 9 million transitTon takes place at x/c = 0.34 
and is caused by S.F. instability phenomena. Experimental data obtained by 
hot-film sensors, shown in figure 10(a), for a Reynolds number of 9 million 
indicate similar results, namely transition takes place at an x/c between 0.3 
to 0.4. At a Reynolds number of 19.2 million (figure 10(c)), transition due 
to S.F. takes place at x/c ~ 0.175; experimental hot film trace data for the 
same freestream conditions, shown in figure 10(b), also indicate that 
transition occurs between x/c of 0.1 to 0.2, 


Design and Parametric Studies for Laminar Glove of X-29 Aircraft 

It has been observed from theoretical studies and experimental data that 
wing sweep have a pronounced effect on transition location when the freestream 
Mach number and Reynolds number are held constant. It is also believed by 
many, from a conceptual view point, that it is easier to achieve large 
laminar-boundary-layer-flow regions on a forward swept wing (FSW) than on an 
aft swept wing (ASW). This is so because for FSW the local values of sweep 
are lower near the leading edge and the local sweep angle increases from the 
leading edge to the trailing edge; whereas, for ASW, sweep angle is the 
largest near the leading edge and decreases progressively toward the trailing 
edge. Figure 11(a) shows schematically a comparison of planforms for 
"equivalent" FSW and ASW. The word "equivalent" implies that for both FSW and 
ASW, wing area, aspect ratio, taper ratio, and shock locations are 
identical. Figure 11(b) shows the plots of local sweep angle versus x/c for 
FSW and ASW of aspect ratio of 4 and taper ratio of 0.4. The planform and 
variation of sweep angle with x/c shown in figures 11(a) and 11(b) are a rough 
approximation of the X-29 aircraft wing. The point C in figure 11(b), which 
is at the intersection of local sweep versus x/c curves for the "equivalent" 
FSW and ASW, indicates that the local sweep angles for these "equivalent" 
wings become equal at x/c ~ 0.5. Thus the possibilities of maintaining 
laminar flow for x/c < 0.5 on an NLF glove for the FSW of X-29 aircraft are 
much greater than they would be for an "equivalent" ASW with a laminar glove. 

It is possible to realize an additional advantage from forward swept 
wings and this stems from the lower shock-wave drag. If the NLF design 
pressure distribution is such that the minimum C p on the upper surface occurs 
aft of x/c ~ 0.5, then the local Mach number at the shock location for FSW 
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will be lower than that ahead of the shock on 
value of the freestream Mach number. This is 
for FSW is higher than that for an equivalent 
(x/c ~ 0.5) as can be seen from figure 11(b). 
lower shock wave drag for FSW for the NLF desi 
moreover, the susceptibility of the NLF design 
is also diminished. This is because the ini ti 
thickness at the shock location would be lower 
equivalent ASW. 


an "equivalent" ASW for the same 
so because the local sweep angle 
ASW downstream of point C 
These phenomena would result in 
gn glove of X-29 aircraft, and 
for FSW to turbulent separation 
al turbulent boundary -layer 
for the FSW than for the 


The above discussion suggests that better prospects exist for overall 
lower total drag for an NLF concept with FSW than one with an equivalent 
ASW. The superior drag peformance due to application of NLF (or HLFC) 
concepts to a FSW arise due to the possibility of achieving larger chordwise 
runs of laminar boundary -layer flow, lower shock wave drag, and less 
susceptibility to turbulent separation drag. 


Figure 12 shows the planform of the X-29 aircraft. The canard and wing q 
are in the same plane. The wing has a leading-edge sweep of approximately 29 
and the trai ling-edge sweep is 45°. The chord of the first flap extends from 
x/c = 0.75 to x/c = 0.9 and that of the second flap extends from x/c - 0.9 to 
x/c = 1.0. The wing leading-edge break is at n = 0.25, and the tip of the 
canard is at n = 0.49. Experimental pressure distributions on the canard 
suggest that boundary -layer flow separation exists on the canard and the 
possibility exist that the separated wake of the canard may prohibit the 
achievement of laminar flow on the wing between the wing root and the midspan 
location. The forward swept wing has associated with it the boundary-layer 
wash-in phenomenon, and for this reason it is believed that a laminar glove 
can be designed so as to obtain extensive laminar flow on the outboard wing 
from n = 0.55 to n = 0.95. 


Simple-sweep methods or modified simple-sweep methods are frequently used 
in the design of airfoil sections for aft swept wings. Different forms of 
modified simple-sweep methods are employed but they give essentially the same 
results. The present modified simple sweep approach for ASW is illustrated in 
figure 13. The geometry of the normal section is arrived at from the geometry 
of the streamwise section by the following expression: 


where. 


_ ^streamwise section 

^ Z//c ^Normal section Cos A^ 

A^ = local value of sweep 
= a L£ + (a je - a le ) • (x/c) 


191 


The pressure distribution is computed for the normal section by airfoil 
programs, such as, Bauer, Garabedian, and Korn program (ref. 9) or Lockheed- 
NASA multi -Component airfoil program (ref. 10). The computations for the 
normal section Mach number and lift coefficient are made using 

M = M • CosA .. 

°° w » effective 

where, 

M = effective component of freesteram Mach number for normal section 
N 

M = streamwise freestream Mach number 

oo 

s 

and, a e ff ec tj ve = A @ x/c = 0.25 for subsonic and transonic freestream 
Mach numbers 


and 

C, _ C, /Cos*” A r f 

= L$ effective 

where. 


Ci = lift coefficient for normal section 
N 

Ci = lift coefficient for streamwise section 

o 


The thickness of the normal section for the aft wing section is thus 
greater than that of the streamwise section by the factor of 1/CosA, . The 
pressure distribution computed for this normal section (with L 

M c« = M « CosA ff and C, = C. / Cos^a ff ) is related to the pressure 
N S l n S e 

distribution for the streamwise section by, 

C p = C p • Cos 2 A, 

Streamwise Normal L 

The application of the above-mentioned modified sweep approach to the 
computation of the pressure distribution for the forward swept wing of the X- 
29 at n = 0.7, M^ = 0.91 and C^ = 0.351 is shown in figure 14. It is 
apparent from this figure that the theoretical pressure distribution computed 
by the above-mentioned modified simple sweep method for aft swept wings does 
not agree with the experimental data for the FSW of X-29 aircraft. 
Consequently a new modified simple sweep approach was developed for the 
computation of pressure distributions for forward swept wings. In this 
approach, the ordinates of the upper surface of the normal section are kept 
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the same as for the streamwise section but the thickness distribution for the 
normal section is increased by the cosine of the local sweep angle and is 
given by 


(t/c) Local = (t/c) Local / Cosa. 

Normal Streamwise 

Streamwise Section 

The remaining steps for the computation of the streamwise pressure 
distribution are the same as for the aft swept wing. Both procedures are 
illustrated in figure 13. The results of computations of the streamwise 
pressure distributions, using the new modified simple sweep method, for the X- 

29 aircraft for n = 0.7, = 0.91 and = 0.351 are shown in figure 14. As 

seen in this figure, the agreement between the theoretical upper surface 
pressure distribution computed by the new modified simple sweep method and 
experimental data is reasonably good. Hence, the parametric studies for the 
design of the outboard laminar-glove sections are conducted using this 
approach. 

In order to enlarge the extent of laminar flow for the X-29 wing, by the 
installation of a wing glove, it is first necessary to determine the extent of 
laminar flow on the existing X-29 wing. Theoretical computations for 
instability and transition locations on the existing X-29 wing at n = 0.7 are 
plotted in figure 15(a) as a function of freestream Reynolds number. 
Experimental pressure distributions are used for boundary conditions in this 
calculation. The computed wake drag coefficient for the existing wing section 
at n = 0.7 is plotted as a function of freestream Reynolds number in figure 
15(b). The results shown in figures 15(a) and 15(b) illustrate that the 
extent of laminar flow on the existing X-29 wing is quite limited and, 
consequently, the wake drag is high for Reynolds numbers in the range of 8 to 

30 mi 1 1 i on . 

A number of pressure distributions and corresponding glove geometries 
have been theoretically investigated in the present studies, under the 
assumption of infinite-span swept-wing conditions. It is necessary that the 
glove geometry be compatible with original X-29 wing such that the glove wraps 
the wing surface with continuous first and second derivatives, and that there 
be a smooth merging of the glove with the original wing surface. From these 
considerations, the glove with the geometrical shape designated as XTM2P3K, 
and shown in figure 16 for spanwise station n = 0.7 was developed. The 
geometry of the baseline X-29 for the spanwise station of n = 0.7 is also 
shown in figure 16 for comparison purposes. Figure 17(a) shows the computed 
streamwise pressure distribution for this glove at n = 0.7 for M = 0.905, 

C. = 0.435 and a leading-edge sweep of 30°* The computed pressure 
S 

distribution for C, = 0.351 for the baseline section is also shown in figure 

17(a) for comparison. It is important to note the shape of C p in figure 

K sonic 
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17(a). For the forward swept tapered wing, the Cp . line has an upward 

some 

slope, whereas for the "equivalent" ASW the C p line will have a downward 

r sonic 

slope. These phenomena suggest, for the specified freestream Mach number and 
value of Cp just upstream of the shock, that the FSW will have a smaller value 
of shock wave drag than the ASW. 

The streamwise pressure distribution computed for the glove with a L.E. 
sweep of 25°, and C L = 0.43 is shown in figure 17(b); since the sweep is 

lower, the design freestream Mach numb er was lowered to a value of 0.865 (as 
compared to an M^ of 0.905 for sweep of 30®). The reason that the Mach number 
was lowered was to keep the leading-edge normal Mach numbers constant. The 
theoretical transition characteristics for the X-29 glove for leading-edge 
sweeps of 30° and 25° and zero spanwise pressure gradients are shown in figure 
17(c). As seen from figure 17(c), the S.F. instability triggers the 
transition at a Reynolds number of 12 million for a leading edge sweep of 30°, 
and at a Reynolds number of 15 million for a leading edge sweep of 25°. 

Figure 17(d) shows a comparison of the wake drag characteristics for the 
X-29 glove for leading edge sweeps of 30° and 25°. Also plotted in figure 
17(d) is the wake drag of the baseline section for a leading-edge sweep of 
30°. These comparisons indicate the following phenomena, first, the wake drag 
of the glove is lower than the basic X-29 wing section by as much as 17 to 18 
counts at a Reynolds number of 10.5 million, and second, when the leading-edge 
sweep is reduced from 30° to 25°, the freestream Reynolds number at which 

Cn becomes minimum is increased from 10.5 million to 16 million. 
u wake 

However, it may not be possible to reduce the wing glove leading-edge sweep to 
25° because of stability and control considerations. 

As stated previously in the theoretical section, the methods of 
references 5, 6, and 7 have been extended in the present studies to take into 
account the effects of boundary-layer wash-in/wash-out and spanwise pressure 
gradients on the development of the laminar boundary layer as well as 
instability and transition locations. Thus, it was possible to conduct 
parametric studies on the effects of boundary-layer wash-in phenomena and 
spanwise pressure gradients on the triggering of transition due to S.F. 
phenomena. Proper use of spanwise gradients can delay transition to higher 
Reynolds numbers and result in a wake drag reduction at the cruise Reynolds 
number. Figures 18(a) and 18(b) show the two types of spanwise pressure 
gradients which have been utilized to study the effect of spanwise pressure 
gradients on transition. In figure 18(a), constant spanwise pressure 
gradients are assumed on the upper surface from 0 < x/c < 0.6; computations 
have been performed for spanwise gradients, AC p , of 0.05, 0.1 and 0.15. For 

c 

the linearly varying spanwise pressure gradients depicted in figure 18(b), the 
spanwise pressure gradient is a maximum near the leading edge and it decreases 
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linearly to zero at x/c = 0.6. Computational results are presented for 

ACp values of 0.05, 0.1, and 0.15. Instabilities and transition 
K 1 

characteristics due to N.F. and S.F. and wake drag characteristics for the 
various spanwise pressure gradients shown in figures 18(a) and 18(b) have been 
determined. A few representative results are shown in figures 19 and 20. 
Figure 19 shows the plots of N.F. and S.F. transition characteristics for a 
constant spanwise pressure gradient corresponding to AC p = 0.15. Also 

c 

plotted in figure 19 for comparison are the results for N.F. and S.F. 
transition characteristics on the laminar-flow glove for a spanwise pressure 
gradient of zero. It can be seen from these results that the Reynolds number 
at which S.F. triggers transition increases from 12 to 20 million as the 
spanwise pressure gradient is changed from zero to AC p = 0.15 (see figure 

c 

18(a)). Furthermore, notice the change in N.F. transition characteristics for 
these two cases. Calculations indicate that turbulent separation does not 
occur in either of the cases shown in figure 19. 

Figure 20 shows the plots of C n versus freestream Reynolds number for 

wake 

the various constant spanwise gradients, including zero, considered in the 
present study. Thus, the results shown in figure 20 illustrate the ability of 
the present theory for increasing the design Reynolds number for minimum wake 
drag by controlling spanwise pressure gradients. It is necessary to point out 
here that constant spanwise pressure gradient increases the shock strength of 
the inboard spanwise location, n = 0.49 and decreases it at the outboard 
station, n = 0.91. It is likely that these shock strengths (figure 18(a)) 
will yield different pressure recoveries than that assumed. 

Another way to obtain a spanwise gradient is to fix the shock strength 
and vary the pressure levels near the leading edge. Figure 21 shows the 
computational results for transition due to N.F. and S.F. phenomena for such a 
pressure distribution i.e., the linearly varying spanwise pressure gradient, 
shown in figure 18(b) (AC p =0.15) . Also shown in figure 21 are the 

theoretical transition characteristics for an X-29 glove when the spanwise 
pressure gradient is assumed to be zero (infinite-span swept wing 
conditions). Results of figure 21 illustrate that the triggering of 
transition due to S.F. is delayed from a Reynolds number of 12 million to 16 
million due to the effect of spanwise pressure gradient. In addition, the 
shock strength remains constant along the span when the linearly varying 
spanwise pressure gradient is assumed. Thus, when the linearly varying 
spanwise pressure gradient is specified, the shock wave drag remains the same 
as for a zero spanwise pressure gradient, and also the adverse effect of 
increased shock strength on the laminar glove design is eliminated. 
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versus freestream Reynolds 


Figure 22 shows the plots of computed C n 

u wake 

number for linearly varying spanwise pressure gradients for several values 
of aC p • Also shown are computations for a spanwise pressure gradient of 
K 1 

zero for comparison purposes. The results shown in figure 22 further 
illustrate how the present theory can be used for delaying transition Reynolds 
number due to triggering of S.F. instabilities by approximately 6 million 
above the zero spanwise-pressure gradient value by the use linearly varying 
spanwise pressure gradients. Neither induced drag, Cp., or shock wave drag 
are adversely affected. 1 


Comparison Between Equivalent Forward and Aft 
Swept Wing Concepts for Laminar Flow 

A comparison of planforms and local sweep variations for FSW and ASW with 
an aspect ratio of 4 and a taper ratio of 0.4 is shown in figure 11. The 
pressure distribution and spanwise pressure gradient selected for making a 
comparison of laminar flow regions on equivalent FSW and ASW is a linearly 
varying spanwise pressure gradients with AC p = 0.15. This spanwise pressure 

_ K 1 

gradient corresponds approximately to that provied by an elliptical spanwise 
load distribution for both FSW and equivalent ASW. The results of 
computations for transition locations as a function of the freestream Reynolds 
number for these "equivalent" wings are shown in figure 23(a) and calculations 

of Cr» versus freestream Reynolds number are shown in figure 23(b). Figure 
u wake 

23(a) shows that extent of laminar flow for FSW is 15 to 35 percent more than 
the equivalent ASW for Reynolds numbers from 8 to 30 million. The wake drag 
for FSW at freestream Reynolds number of 16 million, which corresponds to 
flight condition for the X-29 aircraft at = 0.905 and C|_ = 0.435, is lower 
by 21 counts than the equivalent ASW. 


The comparison shown in figures 23(c) and 23(d) assumes a favorable 
spanwise pressure gradient that provides a spanwise flow that is opposite to 
that norm al ly obtained on ASW; this provides a spanload that deviates from the 
elliptical and causes an increase in Cp. for ASW. Results of computations of 

the extent of laminar flow, shown in figure 23(c), indicate the extent of 
laminar flow for FSW is still larger by 7 to 20 percent over the "equivalent" 
ASW. A comparison of C n values under this assumption, as shown in figure 

23(d), indicates that for ASW the minimum value of C n occurs at a Reynolds 

u wake 

number of 11.5 million, whereas for FSW, C n becomes a minimum at a 

u wake 

freestream Reynolds number of 16.5 million. Moreover, values of C n for 

u wake 

the range of Reynolds number from 8 to 30 million are lower for FSW than the 
equivalent ASW. 
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It is necessary to point out that even though the FSW concept, in 
general, is better for a laminar-flow wing than the equivalent ASW, it is 
possible to determine other pressure distributions and spanwise pressure 
gradients for ASW that yield larger extents of laminar flow than shown in 
figure 23. The insight for obtaining the conditions for large laminar flow 
regions for highly swept ASW, especially for supersonic flow conditions, can 
be obtained by going through an optimizaiton process for chordwise pressure 
distribution and spanwise pressure gradients as has been done for FSW of X-29 
aircraft in the present studies. 


CONCLUDING REMARKS 

(1) The present integral methods for computing laminar boundary-layer 
properties and transition have been extended to take into account boundary- 
layer wash-in/wash-out phenomena and the effects of spanwise pressure gradient 
on transition for the purpose of delaying transition due to S.F. instabilities. 

(2) The above method has been applied to predict transition phenomena on 
the clean-up glove of the F-14 aircraft. Good agreement has been obtained 
between the predicted location of transition and experimental data as a 
function of Reynolds number, Mach number, and wing leading edge sweep. The 
results of the present studies for the F-14 aircraft wing glove, as well as 
those of references 5, 6, and 7, indicate that once transition occurs due to 
S.F. , the transition location moves upstream very fast with even slight 
increase in freestream Reynolds number. Thus, by the use of the present 
method the wing-section geometry can be designed and spanwise pressure 
gradient specified in such a manner that transition due to S.F. pheonomena can 
be substantially delayed. This facilitates obtaining larger laminar-flow 
regions and making the wing performance less susceptible to minor changes in 
Reynolds number and Mach number. 

(3) Comparisons of the computed displacement thicknesses, momentum 
thicknesses, and form factors in turbulent boundary-layer regions on the glove 
of a F-14 aircraft by the present turbulent boundary-layer method with 
experimental data for several freestream conditions and wing sweeps indicated 
good agreement. This implies that not only is the present turbulent boundary- 
layer theory quite reliable, but the transition location prediction is 
accurate as well. 

(4) The extent of laminar flow computed for the F-14 wing glove for a 
leading-edge sweep of 25° and M^ = 0.7 was larger than that computed for the 
smaller sweep value of 20°. The experimental data obtained from hot film 
traces also revealed this unusual pheonomena. This provides additional 
evidence of the ability of the present method to optimize wing sections as 
well as the planforms for obtaining large extents of laminar flow. 


197 


(5) The effective sweep procedure which is normally used for ASW, when 
used to compute the streamwise pressure distribution on the FSW of the X-29 
aircraft, did not yield good correlation with experimental results. The new 
effective sweep procedure established during the present study for FSW gives 
reasonable predictions of streamwise pressure distributions when compared to 
experimental data for the X-29 wing. 

(6) By the use of the method of characteristics, in conjunction with the 

effective sweep procedure established for FSW, section geometries were 
established for the outboard glove for the wing of the X-29 so as to obtain 
the desired chordwise pressure distributions for leading-edge sweeps of 30° 
and 25°. The computed transition and wake drag characteristics were also 
computed and compared for these sweeps. The methods were also used to compute 
transition and wake drag for the basic outboard wing section. This was done 
for assumed "infinite-span swept-wing" conditions. The wake drag for the 
investigated glove section was quite superior to the basic section, and the 
glove with a leading-edge sweep of 25° had a minimum value of C D at 

freestream Reynolds number of 16 million as compared to 10.5 million for a 
sweep of 30°. 

(7) Several spanwise pressure gradients, both constant and linearly 
varying, for x/c = 0 to x/c = 0.6 on the upper surface gloves, were 
investigated for their ability to delay the transition due to S.F. phenomena 
at M = 0.905. It was concluded from these parametric studies that both 
linearly varying and constant-increment spanwise pressure gradients provided 
beneficial effects on skin friction drag at the flight condition considered. 

(8) Computational experiments were performed to determine the extent of 
laminar flow on FSW and "equivalent" ASW for identical chordwise pressure 
distributions and spanwise pressure gradients. This comparison indicated that 
the FSW concept is superior for achieving large regions of laminar flow and 
lower wake drag than an "equivalent" ASW concept. 
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APPENDIX A 


Expression for Make Drag for Infinitely Swept Wing 


C c 




Figure Al. - Coordinate system used for wake drag derivations. 


Consider control volume ABCD-EFGH as shown in figure Al above for a wing 
with sweep angle A. The x-coordinate is perpendicular to the leading edge, 
the z-coordinate is parallel to leading edge, the y-coordinate is 
perpendicular to the wing surface, S is equal to wing span as shown, and c 
is the chord of wing in the direction of flow. 

Net mass balance ~ « 

across the control = / p U S cos A dy - / puS dy = 0 (Al) 

volume -oo -oo 


Drag = Net momentum flux across control volume 


- / plr S cos A dy - / puS(u cos A + w sin A)dy (A2) 


Multiply equation (Al) by and subtract from (A2) 
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00 o o 0 

Drag = / pU u$(cos 2 a + sin 2 A)dy - J 


puS(u cos A + w sin A)dy 


oo 

= s cos A P u 2 / (1 - y|— ) dy 

00 00 J pU ^ U' 

— 00 CO CO oo 

+ S sin A p (1 - £-)dy (A3) 


or 


1 — 7J 


Drag 


wake o- p U S*c cos A 
^ 00 00 


Drag 


12 2 
7 p u sec A*S»c cos A 

^ oo oo 


where 


0 0 
uu 0 uw 0 

_ f oo / oo . / i 

■ 2{— — cos A + — — Sin a} 


o = / _py_ (i . iU dy 
UU J p u ^ u ; J 

OO -oo 00 CO oo 


0 = J _pu_ (1 _ w_) dy 

UW J p u ^ w ; J 
00 —00 00 00 00 


(A4) 


For swept wings as well as bodies of revolution the boundary -layer 
parameters at the trailing edge can be calculated by the theoretical methods 
described in the main body of this paper. Equation (A4) gives Cp in 
terms of conditions downstream at infinity. The following derivatYSkf show 


how boundary-layer parameters at the trailing edge are related to viscous flow 
quantities at infinity. 

The momentum integral equation in the wake in the streamwise direction, 
either for swept wing or for body of revolution can be written as 


where 


dx 




1 d p e 

+ — 9 = 0 

p dx x 
e 


0 =0,,, of equation (A4) 

A U 14 


(A5) 
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By integrating equation (A5) in the wake from the trailing edge to downstream 
infinity, and after algebraic manipulations the following equation is obtained 


u (H +2) °° » 

[-le x p e (^) (] T _ E _ = T / E _ 


(A6) 


The hypothesis is made that in the wake the following relationship is valid. 


An 



An 



U 




F 





(A7) 


The following equation can be derived by making use of Stewartson's trans- 
formation for a Prandtl number equal to one. 


H 


x 



(1 + 


III 

2 



+ 


III 

2 


(A8) 


H x = streamwise form factor in physical plane 

H = streamwise form factor in transformed plane 
X T 

Substitution of equations (A7) and (A8) into (A5) and simplification results 
in the following equation 


0 


uu 


(oo) 


0 


UU 


T.E. 


(1 


Cl - O 


1 M 2 

e,s 

2n -2 .5 


•>-2.5 

^T.E. 



(A9) 


where. 


p = {H 


T.E. 


+ 5 + 


( y - DM ?}/2 


/Nr Root mean square Mach number in the wake from trailing 
e edge to infinity downstream 


Equation (A4) for C n for a swept wing contains the crossflow momentum 

u wake 

thickness term, 0 . In order to evaluate this crossflow component of 

uw 

oo 

momentum thickness, use is made of the momentum equation in the spanwise 
direction for an infinite-span swept wing. Thus, 
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(A10) 


pH 


+ v — + ow — = - i£ + i_ 
ax p ay pw az az ay 


r 3W-1 


The infinite-span swept-wing assumption implies 
a 


az 


( ) = o 


(All) 


and the continuity equation under the same assumption yields. 


p v = - / g^- (pu)dy 


0 


Along the locus of minimum velocity in the wake fy = 0. Substituting (All) and 

(A12) into (A10) and integrating with respect to y in the wake from 0 to 
6 z> the following form of momentum integral equation in the wake can be 
derived. 


(A12 ) 


siw,/ ^d-wrW-o 


e e 


W [°e u e w e e J * 0 


p e u e W P 

n = r_T.E.^ T.E., r e T.E.,^ ^ 

9 uw P«, u w ^ G uwH.E. 

oo 00 00 00 


(A13) 


Thus, C D can be written from equa^ons (A4), (A9), and (A13) as: 


c ^ 0 uu^T.E. (p H.E. ^u ) cos A T.E. 


wake 


KJt.E. (^T.E. (iT^ s1n **T.E.] 


(A14) 


Equation (A14) contain terms evaluated at the trailing edge, such as 
^ 0 uw^T.E.’ ^ 0 uu^T.E.’ ^ p e^T.E. 5 and ^ u e^T.E.' These values are computed in 
the framework of computer programs for boundary -layer methods, (refs. 5, 6, 7, 
and 8) and pressure-di stribution prediction methods such as those of references 
9 and 10. 
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M co sin A 



u(s,^ )— Velocity profile in 

streamwise direction 

v(s, £ )- Velocity profile in 

direction normal to wall 

w(s,^ )-Velocity profile in 
spanwise direction 


Figure 1. - Definition of coordinate axis system. 



AjBjC^D - Edge of boundary layer 

F z - Algebraic sum of forces in z directions 


Figure 2. - Schematic of control volume used to account for boundary-layer 
wash-in/wash-out phenomena. 



Figure 3. - Criteria for determining laminar boundary-layer neutral 
instability. 



x 


Figure 4. - Criteria for determining laminar boundary-layer transition 




Free-stream turbulence, % 

Figure 5. - Effect of free-stream turbulence on the difference of transition 
and instability Reynolds numbers. 



• Three rows of flush static pressure orifices 

• Three rows of surface pitot tubes 

• Five hot-film gauges 

• Two boundary-layer rakes 


Figure 6. - Planform of an instrumentation on the gloved F- 14 wing. 
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(M m = 0.700, A L £ = 20.0°, R«= 13 X 10 6 ) 


n = 0.52 


11 = 0.67 


n = 0.84 




HOT FILM TRACES 

Station 2 


1.0 


(a) Experimental pressure distributions and hot-film traces. 


(M „ =0.697, A le = 20°) 



(b) Computed transition locations as a function of free-stream chord Reynolds 
number. 

Figure 7. - Results for the F-14 "clean-up" gloved wing for a leading- edge 
sweep of 20° and of 0.7. 
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( M 00 = 0.697, A LE = 20° , Rco = 13.6 x 10 6 ) 



X/C 

(c) Comparison between computed and flight boundary-layer parameters for 
R*, = 13 x 10°. 


Figure 7. - Concluded. 


r| - 0.52 


ti = 0.67 


T1 = 0.84 



HOT FILM TRACES 

Station 3 


(a) Experimental pressure distributions and hot-film traces for 
R = 11 x 10 6 . 

CO 


Tl = 0.52 n = 0.67 T] = 0.84 



c 0.1 0.2 0.3 0.4 0.5 


HOT FILM TRACES 

Station 2 

(b) Experimental pressure distributions and hot-film traces for 
R = 23 x 10 6 . 

oo 

Figure 8. - Results for the F - 14 "clean-up" gloved wing for a leading-edge 
sweep of 20° and of 0.8. 
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NF (Normal flow) 


Transition 


SF (Spanwise flow) 


Flight 8 , Station 3 


Flight 18, Station 2 


Ins tabili ty. 


Rep/10 6 

(c) Computed transition and instability locations versus free-stream Reynolds 
number. 


(A le = 20° > Station 3) 


(d) Computed C n versus free-stream chord Reynolds number 
u wake 

Figure 8. - Continued. 





(e) Comparison between computed and flight boundary- layer parameters for 
R = 12.4 x 10 6 . 

CO 

Figure 8. - Concluded. 
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Tl = 0.52 


H = 0.67 


T\ = 0.84 




HOT FILM TRACES 

Station 3 


J I 

.8 1.0 


(a) Experimental pressure distributions and hot-film traces. 



(b) Computed transition locations as a function of free-stream chord Reynolds 
number. 

Figure 9. - Results for the F- 14 "clean-up" gloved wing for a leading-edge 
sweep of 25° and of 0.7. 
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(c) Comparison 

R = 10.8 x 
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Flight 5, Station 3 
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Figure 9. - Concluded. 


11 = 0.52 


n = 0.67 


ri = 0.84 



c 0.1 0.2 0.3 0.4 0.5 

HOT FILM TRACES 


(a) Experimental pressure distributions and hot-film traces for 
= 9.5 x 10 6 . 


i 

i 


n = 0.52 Ti = 0.67 Ti = 0.84 
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HOT FILM TRACES 


(b) Experimental pressure distributions and hot-film traces for 
R w = 19.2 x 10 6 . 

Figure 10. - Results for the F - 1 4 "clean-up" gloved wing for a leading-edge 
sweep of 30° and of 0.7. 
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(c) Computed transition and instability locations versus free-stream Reynolds 
number. 

Figure 10. - Concluded. 
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Figure 11. - Comparison between "equivalent" forward and aft swept wings. 
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x/c 

Figure 14. - Computational results for the pressure distributions on the 
X-29 wing section situated on both forward and aft winqs 
(M, = 0.91, r, = 0.70, a le = 30°) y ' 
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(a) Comparison between baseline and NLF-Glove for - 0.905 and leading-edge 
sweep of 30°. 



(b) Pressure distribution for = 0.865 and leading-edge sweep of 25 . 

Figure 17. - Chordwise pressure distributions for X-29 wing section for 
„ = 0.7 and transition location computations. 
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(c) Transition characteristics for leading-edge sweeps of 30° and 25 



Figure 17. - Concluded. 




-.8 


Cp, upper "* 4 



(b) Linearly varying spanwise pressure gradient. 

Figure 18. - Schematics showing types of spanwise pressure gradients studied. 
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NF (Normal flow), AC p =0 
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SF, ACp = .15 — \ 
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SF (Spanwise flow), ACp =0 

.1 - Kc 


Rco /1q6 


Figure 19, - Comparison of the effect of a constant spanwise pressure 

gradient of aC p = 0.15 on transition with that for zero 
c 

spanwise gradient, i.e., with aC p = 0. 

c 


(M „ = 0.905, A le =300) 


C-D wake 22 


AC = 0 




Figure 20 


versus with ACp as parameter for X-29 glove 











NF (Normal flow), ASW 


NF,FSW 


SF (Spanwise flow), FSW 


X/C .3 h SF, ASW 


ASW - Aft Swept Wing A le = 45 ° 

■1 “ FSW - Forward Swept Wing A le = 30 ° 


R„/10 6 

(c) Transition location versus for AC p ^ = -0.15. 


Forward Swept Wing A|_e ? 30° 
44 T Aft Swept Wing A |_£ = 45° — y 


versus R m for aC d 
ke 

Figure 23. - Concluded 
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